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A concept - quantum order - is introduced to describe a new kind of orders that generally 
appear in quantum states at zero temperature. Quantum orders that characterize universality 
classes of quantum states (described by complex ground state wave-functions) is much richer then 
classical orders that characterize universality classes of finite temperature classical states (described 
by positive probability distribution functions). The Landau's theory for orders and phase transitions 
does not apply to quantum orders since they cannot be described by broken symmetries and the 
associated order parameters. We introduced a mathematical object - projective symmetry group - 
to characterize quantum orders. With the help of quantum orders and projective symmetry groups, 
we construct hundreds of symmetric spin liquids, which have SU{2), (7(1) or Z2 gauge structures 
at low energies. We found that various spin liquids can be divided into four classes: (a) Rigid 
spin liquid - spinous (and all other excitations) are fully gaped and may have bosonic, fermionic, 
or fractional statistics, (b) Fermi spin liquid - spinous are gapless and are described by a Fermi 
liquid theory, (c) Algebraic spin liquid - spinons are gapless, but they are not described by free 
fermionic/bosonic quasiparticles. (d) Bose spin liquid - low lying gapless excitations are described 
by a free boson theory. The stability of those spin liquids are discussed in details. We find that stable 
2D spin liquids exist in the first three classes (a-c). Those stable spin liquids occupy a finite region 
in phase space and represent quantum phases. Remarkably, some of the stable quantum phases 
support gapless excitations even without any spontaneous symmetry breaking. In particular, the 
gapless excitations in algebraic spin liquids interact down to zero energy and the interaction does 
not open any energy gap. We propose that it is the quantum orders (instead of symmetries) that 
protect the gapless excitations and make algebraic spin liquids and Fermi spin liquids stable. Since 
high Tc superconductors are likely to be described by a gapless spin liquid, the quantum orders and 
their projective symmetry group descriptions lay the foundation for spin liquid approach to high Tc 
superconductors. 
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I. INTRODUCTION 

Due to its long length, we would like to first outline 
the structure of the paper so readers can choose to read 
the parts of interests. The section M summarize the main 
results of the paper, which also serves as a guide of the 
whole paper. The concept of quantum order is introduced 
in section I A. A concrete mathematical description of 
quantum order is described in section IVA and section 



IV B. Readers who are interested in the background and 
motivation of quantum orders may choose to read sec- 
tion I A. Readers who are familiar with the slave-boson 
approach and just want a quick int roduc tion to qu antum 
orders may choose to read sections IV A and IV B , Read- 



ers who are not familiar with the slave-boson approach 
may find the review sections II and III useful. Read- 



ers who do not care about the slave-boson approach but 
are interested in application to high Tc superconductors 
and experimental measurements of quantum orders may 
choose to read sections I A, IB, VII and Fig. 0- Fig. |l^, 
to gain some intuitive picture of spinon dispersion and 
neutron scattering behavior of various spin liquids. 



A. Topological orders and quantum orders 

Matter can have many different states, such as gas, 
liquid, and solid. Understanding states of matter is the 
first step in understanding matter. Physicists find matter 
can have much more different states than just gas, liquid, 
and solid. Even solids and liquids can appear in many 
different forms and states. With so many different states 
of matter, a general theory is needed to gain a deeper 
understanding of states of matter. 

All the states of matter are distinguished by their in- 
ternal structures or orders. The key step in developing 
the general theory for states of matter is the realization 
that all the orders are associated with symmetries (or 
rather, the breaking of symmetries). Based on the rela- 
tion between orders and symmetries. Landau developed a 
general theory of orders and the transitions between dif- 
ferent orders. ||l|, || Landau's theory is so successful and 
one starts to have a feeling that we understand, at in 
principle, all kinds of orders that matter can have. 

However, nature never stops to surprise us. In 1982, 
Tsui, Stormer, and GossardH discovered a new kind of 
state - Fractional Quantum Hall (FQH) liquid. Q Quan- 
tum Hall liquids have many amazing properties. A quan- 
tum Hall liquid is more "rigid" than a solid (a crystal), 
in the sense that a quantum Hall liquid cannot be com- 
pressed. Thus a quantum Hall liquid has a fixed and well- 
defined density. When we measure the electron density 
in terms of filling factor v, we found that all discovered 
quantum Hall states have such densities that the filling 
factors are exactly given by some rational numbers, such 
as i^ = 1,1/3,2/3,2/5,.... Knowing that FQH liquids 
exist only at certain magical filling factors, one cannot 
help to guess that FQH liquids should have some inter- 
nal orders or "patterns". Different magical filling fac- 
tors should be due to those different internal "patterns" . 
However, the hypothesis of internal "patterns" appears 
to have one difficulty - FQH states are liquids, and how 
can liquids have any internal "patterns" ? 

In 1989, it was realized that the internal orders in 
FQH liquids (as well as the internal orders in chiral spin 
liquids ra, o) are different from any other known orders 
and cannot be observed and characterized in any con- 
ventional ways.JTl B What is really new (and strange) 
about the orders in chiral spin liquids and FQH liquids 
is that they are not associated with any symmetries (or 
the breaking of symmetries) , and cannot be described by 
Landau's theory using physical order parameters. H] This 
kind of order is called topological order. Topological order 
is a new concept and a whole new theory was developed 
to describe it.p, [lO| 

Knowing FQH liquids contain a new kind of order - 
topological order, we would like to ask why FQH liquids 
are so special. What is missed in Landau's theory for 
states of matter so that the theory fails to capture the 
topological order in FQH liquids? 

When we talk about orders in FQH liquids, we are 
really talking about the internal structure of FQH liq- 



uids at zero temperature. In other words, we are talking 
about the internal structure of the quantum ground state 
of FQH systems. So the topological order is a property 
of ground state wave-function. The Landau's theory is 
developed for system at finite temperatures where quan- 
tum effects can be ignored. Thus one should not be sur- 
prised that the Landau's theory does not apply to states 
at zero temperature where quantum effects are impor- 
tant. The very existence of topological orders suggests 
that finite-temperature orders and zero-temperature or- 
ders are different, and zero-temperature orders contain 
richer structures. We see that what is missed by Landau's 
theory is simply the quantum effect. Thus FQH liquids 
are not that special. The Landau's theory and symmetry 
characterization can fail for any quantum states at zero 
temperature. As a consequence, new kind of orders with 
no broken symmetries and local order parameters (such 
as topological orders) can exist for any quantum states at 
zero temperature. Because the orders in quantum states 
at zero temperature and the orders in classical states at 
finite temperatures are very different, here we would like 
to introduce two concepts to stress their differences :pl| 

(A) Quantum orders:\92\ which describe the universal- 
ity classes of quantum ground states {ie the universality 
classes of complex ground state wave- functions with in- 
finity variables); 

(B) Classical orders: which describe the universality 
classes of classical statistical states {ie the universality 
classes of positive probability distribution functions with 
infinity variables). 

From the above definition, it is clear that the quantum 
orders associated with complex functions are richer than 
the classical orders associated with positive functions. 
The Landau's theory is a theory for classical orders, 
which suggests that classical orders may be characterized 
by broken symmetries and local order parameters. [p3| 
The existence of topological order indicates that quan- 
tum orders cannot be completely characterized by bro- 
ken symmetries and order parameters. Thus we need to 
develop a new theory to describe quantum orders. 

In a sense, the classical world described by positive 
probabilities is a world with only "black and white" . The 
Landau's theory and the symmetry principle for classical 
orders are color blind which can only describe different 
"shades of grey" in the classical world. The quantum 
world described by complex wave functions is a "colorful" 
world. We need to use new theories, such as the theory of 
topological order and the theory developed in this paper, 
to describe the rich "color" of quantum world. 

The quantum orders in FQH liquids have a special 
property that all excitations above ground state have fi- 
nite energy gaps. This kind of quantum orders are called 
topological orders. In general, a topological order is de- 
fined as a quantum order where all the excitations above 
ground state have finite energy gapes. 

Topological orders and quantum orders are general 
properties of any states at zero temperature. Non trivial 
topological orders not only appear in FQH liquids, they 



also appear in spin liquids at zero temperature. In fact, 
the concept of topological order was first introduced in a 
study of spin liquids, [pi FQH liquid is not even the first 
experimentally observed state with non trivial topolog- 
ical orders. That honor goes to superconducting state 
discovered in I9II.|Q In contrast to a common point 
of view, a superconducting state cannot be characterized 
by broken symmetries. It contains non trivial topological 
orders, 1 13 and is fundamentally different from a super- 
fluid state. 

After a long introduction, now we can state the main 
subject of this paper. In this paper, we will study a new 
class of quantum orders where the excitations above the 
ground state are gapless. We believe that the gapless 
quantum orders are important in understanding high Tc 
superconductors. To connect to high Tc superconduc- 
tors, we will study quantum orders in quantum spin liq- 
uids on a 2D square lattice. We will concentrate on how 
to characterize and classify quantum spin liquids with 
different quantum orders. We introduce projective sym- 
metry groups to help us to achieve this goal. The projec- 
tive symmetry group can be viewed as a generalization 
of symmetry group that characterize different classical 
orders. 



B. Spin-liquid approach to high Tc superconductors 

There are many different approaches to the high Tc 
superconductors. Different people have different points 
of view on what are the key experimental facts for the 
high Tc superconductors. The different choice of the key 
experimental facts lead to many different approaches and 
theories. The spin liquid approach is based on a point of 
view that the high Tc superconductors are doped Mott 
insulators. pj, nS, n& (Here by Mott insulator we mean a 
insulator with an odd number of electron per unit cell.) 
We believe that the most important properties of the 
high Tc superconductors is that the materials are in- 
sulators when the conduction band is half filled. The 
charge gap obtained by the optical conductance experi- 
ments is about 2eV, which is much larger than the anti- 
ferromagnetic (AF) transition temperature Taf ^ 250K, 
the superconducting transition temperature Tc ^ f OOK, 
and the spin pseudo-gap scale A ~ 40mcV. p7|, [iq , [l9| ] 
The insulating property is completely due to the strong 
correlations present in the high Tc materials. Thus the 
strong correlations are expect to play very important role 
in understanding high Tc superconductors. Many impor- 
tant properties of high Tc superconductors can be directly 
linked to the Mott insulator at half filling, such as (a) the 
low charge density [g^ and superfiuid density, pi|| (b) Tc 
being proportional to doping Tc oc x,^, p3| |24| (c) the 
positive charge carried by the charge carrier, pO| etc . 

In the spin liquid approach, the strategy is to try to 
understand the properties of the high Tc superconduc- 
tors from the low doping limit. We first study the spin 
liquid state at half filling and try to understand the par- 



ent Mott insulator. (In this paper, by spin liquid, we 
mean a spin state with translation and spin rotation 
symmetry.) At half filling, the charge excitations can 
be ignored due to the huge charge gap. Thus we can 
use a pure spin model to describe the half filled system. 
After understand the spin liquid, we try to understand 
the dynamics of a few doped holes in the spin liquid 
states and to obtain the properties of the high Tc su- 
perconductors at low doping. One advantage of the spin 
liquid approach is that experiments (such as angle re- 
solved photo-emission, 0, |l|, |5[ |26) NMR,[pi, neutron 



scattering, |2q, 29 



30| etc ) suggest that underdoped cu- 
perates have many striking and qualitatively new prop- 
erties which are very different from the well known Fermi 
liquids. It is thus easier to approve or disapprove a new 
theory in the underdoped regime by studying those qual- 
itatively new properties. 

Since the properties of the doped holes (such as their 
statistics, spin, effective mass, etc ) are completely de- 
termined by the spin correlation in the parent spin liq- 
uids, thus in the spin liquid approach, each possible spin 
liquid leads to a possible theory for high Tc supercon- 
ductors. Using the concept of quantum orders, we can 
say that possible theories for high Tc superconductors in 
the low doping limits are classified by possible quantum 
orders in spin liquids on 2D square lattice. Thus one 
way to study high Tc superconductors is to construct all 
the possible spin liquids that have the same symmetries 
as those observed in high Tc superconductors. Then an- 
alyze the physical properties of those spin liquids with 
dopings to see which one actually describes the high Tc 
superconductor. Although we cannot say that we have 
constructed all the symmetric spin liquids, in this paper 
we have found a way to construct a large class of symmet- 
ric spin liquids. (Here by symmetric spin liquids we mean 
spin liquids with all the lattice symmetries: translation, 
rotation, parity, and the time reversal symmetries.) We 
also find a way to characterize the quantum orders in 
those spin liquids via projective symmetry groups. This 
gives us a global picture of possible high Tc theories. 
We would like to mention that a particular spin liquid 
- the staggered-flux/d-wave state] 31, |2| - may be im- 
portant for high Tc superconductors. Such a state can 
explain p^, pj| the highly unusual pseudo-gap metallic 
state found in underdoped cuperates,|17, [iq , 2^, |2^ as 
well as the d-wave superconducting state p2r 

The spin liquids constructed in this paper can be di- 
vided into four class: (a) Rigid spin liquid - spinons are 
fully gaped and may have bosonic, fermionic, or frac- 
tional statistics, (b) Fermi spin liquid - spinons are gap- 
less and are described by a Fermi liquid theory, (c) Alge- 
braic spin liquid - spinons are gapless, but they are not 
described by free fermionic/bosonic quasiparticles. (d) 
Bose spin liquid - low lying gapless excitations are de- 
scribed by a free boson theory. We find some of the con- 
structed spin liquids are stable and represent stable quan- 
tum phases, while others are unstable at low energies due 
to long range interactions caused by gauge fluctuations. 




The algebraic spin liquids and Fermi spin liquids are in- 
teresting since they can be stable despite their gapless 
excitations. Those gapless excitations are not protected 
by symmetries. This is particularly striking for algebraic 
spin liquids since their gapless excitations interact down 
to zero energy and the states are still stable. We propose 
that it is the quantum orders that protect the gapless 
excitations and ensure the stability of the algebraic spin 
liquids and Fermi spin liquids. 

We would like to point out that both stable and unsta- 
ble spin liquids may be important for understanding high 
Tc superconductors. Although at zero temperature high 
Tc superconductors are always described stable quantum 
states, some important states of high Tc superconductors, 
such as the pseudo-gap metallic state for underdoped 
samples, are observed only at finite temperatures. Such 
finite temperature states may correspond to (doped) un- 
stable spin liquids, such as staggered flux state. Thus 
even unstable spin liquids can be useful in understand- 
ing finite temperature metallic states. 

There are many different approach to spin liquids. In 
addition to the slave-boson approach, [S 0^ 0^, ISll Isl, 

I—— I I I I I I I I I I I I L T— 1 I I I I I I I I 

|33| , p5| , |3q , |37| , p8| , |39| , [40| | spin liquids has been stud- 
ied using slave-fermion/cr-model approach 
ES, M quantum dimer model, p7 
various numerical methods. |52[ p3| 
the numerical results and recent experimental results |56| 
strongly support the existence of quantum spin liquids in 
some frustrated systems. A 3D quantum orbital liquid 
was also proposed to exist in XgTiOs .IpTII 

However, I must point out that there is no generally 
accepted numerical results yet that prove the existence 
of spin liquids with odd number of electron per unit cell 
for spin-1/2 systems, despite intensive search in last ten 
years. But it is my faith that spin liquids (with odd 
number of electron per unit cell) exist in spin-1/2 sys- 
tems. For more general systems, spin liquids do ex- 
ist. Read and SachdevEj found stable spin liquids in 
a Sp{N) model in large A^ limit. The spin-1/2 model 
studied in this paper can be easily generalized to SU{N) 
model with N/2 fermions per site.[^ Q In the large 
N limit, one can easily construct various Hamiltonians 
whose ground states realize various U{1) and Z2 spin 
liquids constructed in this paper. |53] The quantum or- 
ders in those large- A^ spin liquids can be described by 
the methods introduced in this paper. Thus, despite the 
uncertainty about the existence of spin-1/2 spin liquids, 
the methods and the results presented in this paper are 
not about (possibly) non-existing "ghost states". Those 
methods and results apply, at least, to certain large- Af 
systems. In short, non-trivial quantum orders exist in 
theory. We just need to find them in nature. (In fact, 
our vacuum is likely to be a state with a non-trivial quan- 
tum order, due to the fact that light exists. |5^) Knowing 
the existence of spin liquids in large- A^ systems, it is not 
such a big leap to go one step further to speculate that 
spin liquids exist for spin-1/2 systems. 



C. Spin-charge separation in (doped) spin liquids 

Spin-charge separation and the associated gauge the- 
ory in spin liquids (and in doped spin Hquids) are very 
important concepts in our attempt to understand the 
properties of high Tc superconductors. |14|, |l5t [l6 , 39, 6^] 
However, the exact meaning of spin-charge separation is 
different for different researchers. The term "spin-charge 
separation" has at lease in two different interpretations, 
fn the first interpretation, the term means that it is better 
to introduce separate spinous (a neutral spin-f /2 excita- 
tion) and holons (a spinless excitation with unit charge) 
to understand the dynamical properties of high T^ su- 
perconductors, instead of using the original electrons. 
However, there may be long range interaction (possibly, 
even confining interactions at long distance) between the 
spinous and holons, and the spinous and holons may not 
be well defined quasiparticles. We will call this inter- 
pretation pseudo spin-charge separation. The algebraic 
spin liquids have the pseudo spin-charge separation. The 
essence of the pseudo spin-charge separation is not that 
spin and charge separate. The pseudo spin-charge sepa- 
ration is simply another way to say that the gapless exci- 
tations cannot he described by free fermions or bosons. In 
the second interpretation, the term "spin-charge separa- 
tion" means that there are only short ranged interactions 
between the spinous and holons. The spinous and holons 
are well defined quasiparticles at least in the dilute limit 
or at low energies. We will call the second interpretation 
the true spin-charge separation. The rigid spin liquids 
and the Fermi spin liquids have the true spin-charge sep- 
aration. 

Electron operator is not a good starting point to de- 
scribe states with pseudo spin-charge separation or true 
spin-charge separation. To study those states, we usu- 
ally rewrite the electron operator as a product of sev- 
eral other operators. Those operators are called parton 
operators. (The spinon operator and the holon opera- 
tor are examples of parton operators). We then con- 
struct mean-field state in the enlarged Hilbert space of 
partons. The gauge structure can be determined as the 
most general transformations between the partons that 
leave the electron operator unchanged. |Q After identi- 
fying the gauge structure, we can project the mean-field 
state onto the physical {ie the gauge invariant) Hilbert 
space and obtain a strongly correlated electron state. 
This procedure in its general form is called projective 
construction. It is a generalization of the slave-boson 
approach. |l|, |l|, ||, |6[ |7[ H, |§ The general pro- 
jective construction and the related gauge structure has 
been discussed in detail for quantum Hall states. Q Now 
we see a third (but technical) meaning of spin-charge 
separation: to construct a strongly correlated electron 
state, we need to use partons and projective construc- 
tion. The resulting effective theory naturally contains a 
gauge structure. 

Although, it is not clear which interpretation of spin- 
charge separation actually applies to high Tc supercon- 



ductors, the possibility of true spin-charge separation in 
an electron system is very interesting. The first con- 
crete example of true spin-charge separation in 2D is 
given by the chiral spin liquid state, ||5|, |6|] where the 
gauge interaction between the spinous and holons be- 
comes short-ranged due to a Chern-Simons term. The 
Chern-Simons term breaks time reversal symmetry and 
gives the spinous and holons a fractional statistics. Later 
in 1991, it was realized that there is another way to 
make the gauge interaction short-ranged through the 
Anderson-Higgs mechanism. [08l H This led to a mean- 



field theory 1 38, ^ of the short-ranged Resonating Va- 
lence Bound (RVB) state [^, Q conjectured earher. We 
will call such a state Z2 spin liquid state, to stress the 
unconfined Z2 gauge field that appears in the low energy 
effective theory of those spin liquids. (See remarks at the 
end of this section. We also note that the Z2 spin liquids 
studied in Ref. [Ej all break the 90° rotation symmetry 
and arc different from the short-ranged RVB state stud- 
ied Ref. @ |0[ ^, |||.) Since the Z2 gauge fluctuations 
are weak and are not confining, the spinous and holons 
have only short ranged interactions in the Z2 spin liquid 
state. The Z2 spin liquid state also contains a Z2 vortex- 
like excitation, p^, p2| The spinous and holons can be 
bosons or fermions depending on if they are bound with 
the Z2 vortex. 

Recently, the true spin-charge separation, the Z2 gauge 
structure and the Z2 vortex excitations were also pro- 
posed in a study of quantum disordered superconduct- 
ing state in a continuum model p3| and in a Z2 slave- 
boson approach [p4|. The resulting liquid state (which 
was named nodal liquid) has all the novel properties of 
Z2 spin liquid state such as the Z2 gauge structure and 
the Z2 vortex excitation (which was named vison) . From 
the point of view of universality class, the nodal liquid is 
one kind of Z2 spin liquids. However, the particular Z2 
spin liquid studied in Ref. [^, ^ and the nodal liquid 
are two different Z2 spin liquids, despite they have the 
same symmetry. The spinous in the first Z2 spin liquid 
have a finite energy gap while the spinons in the nodal 
liquid are gapless and have a Dirac-like dispersion. In this 
paper, we will use the projective construction to obtain 
more general spin liquids. We find that one can construct 
hundreds of different Z2 spin liquids. Some Z2 spin liq- 
uids have finite energy gaps, while others are gapless. 
Among those gapless Z2 spin liquids, some have finite 
Fermi surfaces while others have only Fermi points. The 
spinons near the Fermi points can have linear E{k) oc |fc| 
or quadratic E{k) c>c k dispersions. We find there are 
more than one Z2 spin liquids whose spinons have a mass- 
less Dirac-like dispersion. Those Z2 spin liquids have 
the same symmetry but different quantum orders. Their 
ansatz are give by Eq. (||), Eq. (|^), Eq. (H), etc . 

Both chiral spin liquid and Z2 spin liquid states are 
Mott insulators with one electron per unit cell if not 
doped. Their internal structures are characterized by a 
new kind of order - topological order, if they are gapped 
or if the gapless sector decouples. Topological order is 



not related to any symmetries and has no (local) or- 
der parameters. Thus, the topological order is robust 
against all perturbations that can break any symmetries 
(including random perturbations that break translation 
symmetry) . [pi nOl (This point was also emphasized in 
Rcf. 1621 recently.) Even though there are no order pa- 
rameters to characterize them, the topological orders can 
be characterized by other measurable quantum numbers, 
such as ground state degeneracy in compact space as pro- 
posed in Ref. |g|, ^^. Recently, Ref. [^^ introduced a 
very clever experiment to test the ground state degen- 
eracy associated with the non-trivial topological orders. 
In addition to ground state degeneracy, there are other 
practical ways to detect topological orders. For example, 
the excitations on top of a topologically ordered state 
can be defects of the under lying topological order, which 
usually leads to unusual statistics for those excitations. 
Measuring the statistics of those excitations also allow us 
to measure topological orders. 

The concept of topological order and quantum order 
are very important in understanding quantum spin liq- 
uids (or any other strongly correlated quantum liquids). 
In this paper we are going to construct hundreds of dif- 
ferent spin liquids. Those spin liquids all have the same 
symmetry. To understand those spin liquids, we need to 
first learn how to characterize those spin liquids. Those 
states break no symmetries and hence have no order pa- 
rameters. One would get into a wrong track if trying to 
find an order parameter to characterize the spin liquids. 
We need to use a completely new way, such as topological 
orders and quantum orders, to characterize those states. 

In addition to the above Z2 spin liquids, in this paper 
we will also study many other spin liquids with different 
low energy gauge structures, such as U{1) and SU{2) 
gauge structures. We will use the terms Z2 spin liq- 
uids, U{1) spin liquids, and SU{2) spin hquids to describe 
them. We would like to stress that Z2, U{1), and SU{2) 
here are gauge groups that appear in the low energy ef- 
fective theories of those spin liquids. They should not be 
confused with the Z2, U{1), and SU{2) gauge group in 
slave-boson approach or other theories of the projective 
construction. The latter are high energy gauge groups. 
The high energy gauge groups have nothing to do with 
the low energy gauge groups. A high energy Z2 gauge 
theory (or a Z2 slave-boson approach) can have a low 
energy effective theory that contains SU{2), t/(l) or Z2 
gauge fluctuations. Even the t-J model, which has no 
gauge structure at lattice scale, can have a low energy 
effective theory that contains SU{2), U{1) or Z2 gauge 
fluctuations. The spin liquids studied in this paper all 
contain some kind of low energy gauge fluctuations. De- 
spite their different low energy gauge groups, all those 
spin liquids can be constructed from any one of SU{2), 
U{1), or Z2 slave-boson approaches. After all, all those 
slave-boson approaches describe the same t-J model and 
are equivalent to each other. In short, the high energy 
gauge group is related to the way in which we write down 
the Hamiltonian, while the low energy gauge group is a 



property of ground state. Thus we should not regard Z2 
spin liquids as the spin liquids constructed using Z2 slave- 
boson approach. A Z2 spin liquid can be constructed 
from the C/(l) or SU{2) slave-boson approaches as well. 
A precise mathematical definit ion of the low energy gauge 
group will be given in section [V A . 



D. Organization 

In this paper we will use the method outlined in 
Ref. IgS], |o| to study gauge structures in various spin 
liquid states. In section II we review SU{2) mean-field 
theory of spin liquids. In section P|, we construct sim- 
ple symmetric spin liquids using translationally invari- 
ant ansatz. In section IV, projective symmetry group is 
introduced to characterize quantum orders in spin liq- 
uids. In section M, we study the transition between dif- 
ferent symmetric spin liquids, using the results obtained 
in appendix B, where we find a way to construct all the 
symmetric spin liquids in the neighborhood of some well 
known spin liquids. We also study the spinon spectrum 
to gain some intuitive understanding on the properties of 
the spin liquids. Using the relation between two-spinon 
spectrum and quantum order, we propose, in section 



VII, a practical way to use neutron scattering to mea- 
sure quantum orders. We study the stability of Fermi 



spin liquids and algebraic spin liquids in section VIII 



We find that both Fermi spin liquids and algebraic spin 
liquids can exist as zero temperature phases. This is 
particularly striking for algebraic spin liquids since their 
gapless excitations interacts even at lowest energies and 
there are no free fermionic/bosonic quasiparticle excita- 
tions at low energies. We show how quantum order can 
protect gapless excitations. Appendix A contains a more 
detailed discussion on projective symmetry group, and a 
classification of Z2, C/(l) and SU{2) spin liquids using 
the projective symmetry group. Section H summarizes 
the main results of the paper. 



II. PROJECTIVE CONSTRUCTION OF 2D 

SPIN LIQUIDS - A REVIEW OF SU(2) 

SLAVE-BOSON APPROACH 



In this section, we are going to use projective con- 
struction to construct 2D spin liquids. We are going 
to review a particular projective construction, nameh 
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the SU{2) slave-boson approach. ||15|, |16| , 

E^ The gauge structure discovered by Baskaran and 
Anderson pq| in the slave-boson approach plays a cru- 
cial role in our understanding of strongly correlated spin 
liquids. 

We will concentrate on the spin liquid states of a pure 
spin-1/2 model on a 2D square lattice 



H. 



/ , JijJi ■ '-'j 



spin — / ^ ^ij 



(1) 



where the summation is over different Hnks (ze (ij) and 
{ji) are regarded as the same) and ... represents possi- 
ble terms which contain three or more spin operators. 
Those terms are needed in order for many exotic spin 
hquid states introduced in this paper to become the 
ground state. To obtain the mean-field ground state of 
the spin liquids, we introduce fermionic parton operators 
fia, a = 1, 2 which carries spin 1/2 and no charge. The 
spin operator Si is represented as 



Si = 



r,Jia^(^PJi. 



P 



(2) 



In terms of the fermion operators the Hamiltonian Eq. (|lj) 
can be rewritten as 



H 



2^ 2 "^^^ 

(ij) 



ha-ljajjpjip + ^haJio'JjpJiP 



(3) 



Here we have used crap ■ cTa'p' = '^^ap'Sa'p — SapSa'ff'. 
We also added proper constant terms J2i Hafia and 
J2{ij) fLfiafjpfjP to get the above form. Notice that 
the Hilbert space of Eq. dsh is generated by the parton 
operators /„ and is larger than that of Eq. (||). The 
equivalence between Eq. (|]) and Eq. (||) is valid only in 
the subspace where there is exactly one fermion per site. 
Therefore to use Eq. (H) to describe the spin state we 
need to impose the constraint [nsl [l6| 



JiaJia -L , Jiajip^al3 U 



(4) 



The second constraint is actually a consequence of the 
first one. 

A mean-field ground state at "zeroth" order is obtained 
by making the following approximations. First we replace 
constraint Eq. (0) by its ground-state average 



ifJafia) = Ij {fiafip^af}) = 



(5) 



Such a constraint can be enforced by including a site 
dependent and time independent Lagrangian multiplier: 

'^o(*)(/ia/i" ~^)t^ — 1;2,3, in the Hamiltonian. At the 
zeroth order we ignore the fluctuations [ie the time de- 



pendence) of ag. If we included the fluctuations of ag 



the constraint Eq. ( 
straint Eq. (B).i5| 




would become the original con- 
ISTII Second we replace the 



operators fl^foP ^'^'^ fiafifd by their ground-state ex- 



pectations value 

r]ij<^al3 = - 2(fia fjp), 






(6) 



Xij /Vj2 

again ignoring their fluctuations. In this way we obtain 
the zeroth order mean-field Hamiltonian: 

O r- 

(ij) 

-\Xij\'^\Vij\'] (7) 

+ X] ["-oifLfic^ - 1) + [(«() + iaDfiafifitali "f h.C 



Xij and rjij in Eq. (^ must satisfy the self consistency 
condition Eq. (H) and the site dependent fields ag(i) are 
chosen such that Eq. (g) is satisfied by the mean-field 
ground state. Such Xij, '^ij and Og give us a mean-field 
solution. The fiuctuations in Xiji Vij and ag(i) describe 
the collective excitations above the mean-field ground 
state. 

The Hamiltonian Eq. (0) and the constraints Eq. 



have a local 51/(2) symmetry. [36, M The local SU{2) 



symmetry becomes explicit if we introduce doublet 



and matrix 



^ 



Uij^ 



\.ij 



Hj 



Vij 

-Xij, 



-u}i 



(8) 



(9) 



Using Eq. (g) and Eq. 
Eq. i) as: 



we can rewrite Eq. (^ and 



4^'i^i 



= 



(10) 



Hr, 



X.s'^^J 



(ij) 



It'^uI Uij) ~ iijlUij^, + h.C.) 



Y^ a\)ipy^i 



(11) 



where r' , / = 1, 2, 3 are the Pauli matrices. From Eq. ( |Tl| ) 
we can see clearly that the Hamiltonian is invariant under 
a local SU{2) transformation W{i): 



iji -^ W{i) i>i 



Ui. 



wii) Uij wHj) 



(12) 



The SU{2) gauge structure is originated from Eq. (||). 
The SU{2) is the most general transformation between 
the partons that leave the physical spin operator un- 
changed. Thus once we write down the parton expres- 
sion of the spin operator Eq. (S) , the gauge structure of 
the theory is determined. ||6l[] (The SU{2) gauge structure 
discussed here is a high energy gauge structure.) 

We note that both components of if) carry spin-up. 
Thus the spin-rotation symmetry is not explicit in our 
formalism and it is hard to tell if Eq. (11) describes a 
spin-rotation invariant state or not. In fact, for a general 
Uij satisfying Uij — t/J^, Eq. ( [ll| ) may not describe a 
spin-rotation invariant state. However, if Uij has a form 

Uij = ipijWij, 



Pij 
Wii 



real number, 

5[/(2), 



(13) 



then Eq. (|ll|) will describe a spin-rotation invariant state. 
This is because the above Uij can be rewritten in a form 
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Eq. @. In this case Eq. (11) can be rewritten as Eq. (|^) 
where the spin-rotation invariance is expHcit. 

To obtain the mean-field theory, we have enlarged the 
Hilbert space. Because of this, the mean-field theory 

is not even qualitatively correct. Let l^mean) be the 
ground state of the Hamiltonian Eq. ( [ll| ) with energy 
E{Uij,a\T^). It is clear that the mean-field ground state 
is not even a valid wave-function for the spin system 
since it may not have one fermion per site. Thus it is 
very important to include fluctuations of Oq to enforce 
one-fermion-per-site constraint. With this understand- 
ing, we may obtain a valid wave-function of the spin sys- 
tem ^ spin{{ai\) by projecting the mean-field state to the 
subspace of one-fermion-per-sitc: 



^spin{{ai}) 



\\{h 



I* 



(U.: 



(14) 



Now the local SU{2) transformation Eq. (n3) can have 

a very physical meaning: I'i'mian) and |*,'Sn^''^ ^^^') 
give rise to the same spin wave-function after projection: 



n 



JiOi I ^mean/ — 



n 






(15) 



Thus Uij and C/^ = W {i)UijW'' {j) are just two differ- 
ent labels which label the same physical state. Within the 
mean-field theory, a local SU{2) transformation changes 

a mean-field state \'^mean) to a different mean-field state 

l^mean). If the two mcan-ficld states always have the 
same physical properties, the system has a local SU{2) 
symmetry. However, after projection, the physical spin 
quantum state described by wave-function ^spind^i}) is 
invariant under the local SU{2) transformation. A local 
SU{2) transformation just transforms one label, Uij, of a 
physical spin state to another label, f/^ • , which labels the 
exactly the same physical state. Thus after projection, 
local SU(2) transformations become gauge transforma- 
tions. The fact that Uij and U^- label the same physical 
spin state creates a interesting situation when we con- 
sider the fiuctuations of Uij around a mean-field solution 
- some fluctuations of Uij do not change the physical 
state and are unphysical. Those fluctuations arc called 
the pure gauge fluctuations. 

The above discussion also indicates that in order for 
the mean-field theory to make any sense, we must at least 
include the SU{2) gauge (or other gauge) fluctuations 
described by Uq and Wij in Eq. (p^), so that the SU{2) 
gauge structure of the mean-field theory is revealed and 
the physical spin state is obtained. We will include the 
gauge fluctuations to the zeroth-order mean-field theory. 
The new theory will be called the first order mean-field 
theory. It is this first order mean-field theory that rep- 
resents a proper low energy effective theory of the spin 
liquid. 

Here, we would like make a remark about "gauge sym- 
metry" and "gauge symmetry breaking". We see that 
two ansatz Uij and U'^- — W{i)UijW'^{j) have the same 



physical properties. This property is usually called the 
"gauge symmetry" . However, from the above discussion, 
we see that the "gauge symmetry" is not a symmetry. A 
symmetry is about two different states having the same 
properties. Uij and t/^'- are just two labels that label 
the same state, and the same state always have the same 
properties. We do not usually call the same state hav- 
ing the same properties a symmetry. Because the same 
state alway have the same properties, the "gauge symme- 
try" can never be broken. It is very misleading to call the 
Anderson-Higgs mechanism "gauge symmetry breaking" . 
With this understanding, we see that a superconductor is 
fundamentally different from a superfluid. A supcrfiuid 
is characterized by [/(I) symmetry breaking, while a su- 
perconductor has no symmetry breaking once we include 
the dynamical electromagnetic gauge fiuctuations. A su- 
perconductor is actually the first topologically ordered 
state observed in experiments, |T3] which has no symme- 
try breaking, no long range order, and no (local) order 
parameter. However, when the speed of light c = oo, 
a superconductor becomes similar to a superfiuid and is 
characterized by U{1) symmetry breaking. 

The relation between the mean-field state and the 
physical spin wave function Eq. (|l^) allows us to con- 
struct transformation of the physical spin wave-function 
from the mean-field ansatz. For example the mean-field 

state \^mean) with U[j = Ui-ij-i produces a phys- 
ical spin wave-function which is translated by a dis- 
tance I from the physical spin wave-function produced 

by l^mean). Thc physical state is translationally sym- 
metric if and only if the translated ansatz f/^' • and the 
original ansatz Uij are gauge equivalent (it does not re- 
quire C/j'- = Uij). We see that the gauge structure can 
complicates our analysis of symmetries, since the phys- 
ical spin wave-function ^spindai}) may has more sym- 
metries than the mean-field state |^mean) before projec- 
tion. 

Let us discuss time reversal symmetry in more detail. 
A quantum system described by 



ihdt^it) = ff*(i) 



(16) 



has a time reversal symmetry if ^{t) satisfying the equa- 
tion of motion implies that \E'*(— t) also satisfying the 
equation of motion. This requires that H = H* . We 
see that, for time reversal symmetric system, if Vt is an 
eigenstate, then ^* will be an eigenstate with the same 
energy. 

For our system, the time reversal symmetry means that 

if thc mean-field wave function "^mean * is a mean-field 
ground state wave function for ansatz {Uij,a\T^), then 

'^mean'^^ ) will be the mean- field ground state wave 

function for ansatz {Ulj,a\{T^)*). That is 



^ meaji — ^ meaji 



\ ^ mean j 



(17) 



For a system with time reversal symmetry, the mean- field 
energy E{Uij,a\T^) satisfies 



E{Uij,a^T')^E{U:„a^{Tr) 



(18) 



Thus if an ansatz {Uij, a^r'') is a mean- field solution, then 
(t/j*-, a\{T'')*) is also a mean-field solution with the same 
mean-field energy. 

From the above discussion, we see that under the time 
reversal transformation, the ansatz transforms as 



Uij 



i-^T^)ia\T'n^T^) ^ -a[rK (19) 



a^' 



Note here we have included an additional SU{2) gauge 
transformation Wi — —ir^. We also note that under the 
time reversal transformation, the loop operator trans- 



forms as Pc = e*^+*^ 



(-zr2)P5(zr2) 



— „-ie+ie 



We see that the t/(l) flux changes the sign while the 
SU{2) fiux is not changed. 

Before ending this review section, we would like to 
point out that the mean-field ansatz of the spin liquids 
Uij can be divided into two classes: unfrustrated ansatz 
where Uij only link an even lattice site to an odd lattice 
site and frustrated ansatz where Uij are nonzero between 
two even sites and/or two odd sites. An unfrustrated 
ansatz has only pure SU{2) flux through each plaquette, 
while an frustrated ansatz has U{1) flux of multiple of 
tt/2 through some plaquettes in addition to the SU{2) 
flux. 



symmetric spin liquid. First let us introduce Uij: 
3 

■rJijUij = Uij 



(21) 



For translationally invariant ansatz, we can introduce a 
short-hand notation: 



"ij = "-i+j^'' = "- 



^+] 



(22) 



where u^ '' ' are real, u^ is imaginary, t'^ is the identity 
matrix and r^'^''^ are the Pauli matrices. The fermion 
spectrum is determined by Hamiltonian 

H = -J2 {y^Uj-i^J + ^-c-) + J2 ^l«o^Vi (23) 



In fc-space we have 



where /i = 0, 1,2,3, 



'(fc)-E<^''^'' 



(24) 



(25) 



ag = 0, and N is the total number of site. The fermion 
spectrum has two branches and is given by 



E±{k) ^u°{k)±Ea{k) 



^o(fc)=./E(^'(^)-«oF 



(26) 



III. SPIN LIQUIDS FROM TRANSLATIONALLY 
INVARIANT ANSATZ 



In this section, we will study many simple examples 
of spin liquids and their ansatz. Through those simple 
examples, we gain some understandings on what kind of 
spin liquids are possible. Those understandings help us 
to develop the characterization and classification of spin 
liquids using projective symmetry group. 

Using the above SU{2) projective construction, one 
can construct many spin liquid states. To limit ourselves, 
we will concentrate on spin liquids with translation and 
90° rotation symmetries. Although a mean-field ansatz 
with translation and rotation invariance always generate 
a spin liquid with translation and rotation symmetries, 
a mean-field ansatz without those invariances can also 
generate a spin liquid with those symmetries. p4| Because 
of this, it is quite difficult to construct all the translation 
and rotation symmetric spin liquids. In this section we 
will consider a simpler problem. We will limit ourselves 
to spin liquids generated from translationally invariant 
ansatz: 



Ui- 



-1,3+1 



= U, 



iji 



alii) = 



(20) 



In this case, we only need to find the conditions under 
which the above ansatz can give rise to a rotationally 



The constraints can be obtained from 
have a form 



dE„ 



da[. 



N{4,y4,i 



u\k) — Oq 

Eo{k) 



E 

k,E+{k)<0 



u\k) — ag 
^o(fc) 



= and 



= (27) 



which allow us to determine Oq, / = 1,2,3. It is inter- 
esting to see that if u^ = and the ansatz is unfrus- 
trated, then we can simply choose Oq = to satisfy the 
mean-field constraints (since u^{k) ~ —u^{k + (tt, tt)) for 
unfrustrated ansatz). Such ansatz always have time re- 
versal symmetry. This is because Uij and —Uij are gauge 
equivalent for unfrustrated ansatz. 

Now let us study some simple examples. First let us 
assume that only the nearest neighbor coupling u^ and 
Uy are non-zero. In order for the ansatz to describe a 
rotationally symmetric state, the rotated ansatz must be 
gauge equivalent to the original ansatz. One can easily 
check that the following ansatz has the rotation symme- 
try 



= 



XT +f]T 

3 1 

XT -r]T 



(28) 
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since the 90° rotation followed by a gauge transforma- 
tion Wi = ir^ leave the ansatz unchanged. The above 
ansatz also has the time reversal symmetry, since time 



reversal transformation 



Hj 



^ij 



followed by a gauge 



transformation Wi — it leave the ansatz unchanged. 

To understand the gauge fluctuations around the above 
mean-field state, we note that the mean-field ansatz 
may generate non-trivial SU{2) flux through plaquettes. 
Those flux may break SU{2) gauge structure down to 
C/(l) or Z2 gauge structures as discussed in Ref. [pSl Eoj. 
In particular, the dynamics of the gauge fluctuations in 
the break down from SU{2) to Z2 has been discussed in 
detail in Ref. |§. According to Ref. [||, |§, the SU{2) 
flux plays a role of Higgs fields. A non-trivial SU{2) flux 
correspond to a condensation of Higgs fields which can 
break the gauge structure and give SU{2) and/or C/(l) 
gauge boson a mass. Thus to understand the dynamics 
of the gauge fluctuations, we need to find the SU{2) flux. 

The SU{2) flux is defined for loops with a base point. 
The loop starts and ends at the base point. For example, 
we can consider the following two loops Ci^2 with the 
same base point i: Ci = i — > i-\-x — > i+x-\-y -^ i+ij -^ i 
and C2 is the 90° rotation of Ci: C2 = i ^ i + y ^ 
i — x + y-^i — x^i. The SU{2) flux for the two loops 
is defined as 



the projective-symmetry-group classification, the SU{2)- 
gapless ansatz Eq. ( |3(]| ) is labeled by SU2AnO and the 
5t/(2)-linear ansatz Eq. 1^) by SU2BnO (see Eq. (|2|)). 
When X y^ V ¥^ 0, The flux Pc is non trivial. How- 
ever, Pc commute with Pc as long as the two loops 
C and C" have the same base point. In this case the 
SU{2) gauge structure is broken down to a f/(l) gauge 
structure. Pq, [4Q| The gapless spinon still only appear 
at isolated k points. We will call such a state U(\)- 
linear state. (This state was called staggered flux state 
and/or d-wave pairing state in literature.) After a proper 
gauge transformation, the f7(l)-linear state can also be 
described by the ansatz 



UiA+v =«X+ {-YVT^ 



(32) 



where the U{1) gauge structure is explicit. Under the 
projective-symmetry-group classification, such a state is 
labeled by UlCnOln (see Eq. (g) and |VC|). The low 
energy effective theory is described by massless Dirac 
fermions (the spinous) coupled to a f/(l) gauge field. 

The above results are all known before. In the follow- 
ing we are going to study a new class of translation and 
rotation symmetric ansatz, which has a form 



^c\ 



:,U 



t,i-\-y ^z-\-y.z-\-x-\-y ^i-\-x-\-y.t-\-x 



Ui 



Pc 



C2 



""?, i. — o'. ^i. — :7! . 7, — n 



,Ua 



x.z — x+y ^z — x+y.z+y ^t-\-y 



Ui 



ululuyus, 

Uxululuy 



(29) 



As discussed in Ref. K, 40 , if the SU{2) flux Pc for aU 



loops are trivial: Pc oc t°, then the SU{2) gauge struc- 
ture is unbroken. This is the case when x = 77 or when 
77 = in the above ansatz Eq. (ESl). The spinon in the 
spin liquid described by r/ = has a large Fermi surface. 
We will call this state S'C/(2)-gapless state (This state 
was called uniform RVB state in literature). The state 
with X = 77 has gapless spinons only at isolated k points. 
We will call such a state 5'C/(2)-linear state to stress the 
linear dispersion E oc |fc| near the Fermi points. (Such a 
state was called the 7r-flux state in literature). The low 
energy effective theory for the 5'C/(2)-linear state is de- 
scribed by massless Dirac fermions (the spinons) coupled 
to a SU{2) gauge field. 

After proper gauge transformations, the S'[/(2)-gapless 
ansatz can be rewritten as 



Ux = «X 
Uy = ix 



(30) 



and the 5f7(2)-linear ansatz as 



IX 



■.,i+y 



= ii-f^X 



(31) 



In these form, the SU{2) gauge structure is explicit since 
Uij oc iT^ . Here we would also like to mention that under 



=0 



Ux =iT]T" - x{t^ - T^) 

Uy =iT]T" - xir^ + T^) 



(33) 



with X s-nd 77 non-zero. The above ansatz describes the 
5t/(2)-gapless spin Hquid if x = 0, and the 5C/(2)-linear 
spin liquid if 77 = 0. 

After a 90° rotation Rgo, the above ansatz becomes 



i = -^77r"-x(T^+Tl) 

y = iTJT° - x{t^ - T^) 



(34) 



The rotated ansatz is gauge equivalent to the origi- 
nal ansatz under the gauge transformation Gjig„{i) = 
(—)*="(! — ^r^)/^/2. After a parity x -^ —x transforma- 
tion Px , Eq. (p3|) becomes 



- = -ir]T° - x{t^ - T^) 
Uy ^ iT]T° ~ xir^ + T^) 



(35) 



which is gauge equivalent to the original ansatz under the 
gauge transformation Gp^ (i) = {—)'^'=i{T^+T^)/^/2. Un- 
der time reversal transformation T, Eq. (p3) is changed 
to 



Ux ^ -iT]T° + x(r^ - r^) 
Uy = -iT]T° + x{t^ + t'^) 



(36) 



which is again gauge equivalent to the original ansatz un- 
der the gauge transformation Gr(i) = (— )*. (In fact any 
ansatz which only has links between two non-overlapping 
sublattices {ie the unfrustrated ansatz) is time reversal 
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symmetric if a[) = .) To summarize the ansatz Eq. (^) 
is invariant under the rotation Rqq, parity P^, and time 
reversal transformation T, followed by the following gauge 
transformations 

G'pJi)-(-^i(T^+Ti)/V2 
Grii) H-r (37) 

Thus the ansatz Eq. ( p3| ) describes a spin hquid which 
translation, rotation, parity and time reversal symme- 
tries. 

Using the time reversal symmetry we can show that 
the vanishing Oq in our ansatz Eq. (|33| ) indeed satisfy the 



changed to 



constraint Eq. (|27|). This is because c 
time reversal transformation. Thus 



dE„ 



a'o 



da'o 



3 under the 
— when 



for any time reversal symmetric ansatz, including 
the ansatz Eq. (p^). 

The spinon spectrum is given by (see Fig. pk) 



E± = 2?7(sin(fc^) + sin(fcy)) ± 2|x|y^2cos2(fc^) 



2cos^(fcj,) 
(38) 



The spinous have two Fermi points and two small Fermi 
pockets (for small ry). The SU{2) flux is non-trivial. Fur- 
ther more Pc^ and Pc2 do not commute. Thus the SU{2) 
gauge structure is broken down to a Z2 gauge structure 
by the 5f/(2) flux Pc, and Pchj^ @ We will caU the 
spin liquid described by Eq. (|3^) Z2-gapless spin liquid. 
The low energy effective theory is described by massless 
Dirac fermions and fermions with small Fermi surfaces, 
coupled to a Z2 gauge field. Since the Z2 gauge inter- 
action is irrelevant at low energies, the spinous are free 
fermions at low energies and we have a true spin-charge 
separation in the Z2-gapless spin liquid. The Z2-gapless 
spin liquid is one of the Z2 spin liquids classified in ap- 
pendix A. Its projective symmetry group is labeled by 
Z2 A tJl^tJ ^t^t^ or emiivalently by Z2Ax2(12)n (see sec- 
tion |V^ and Eq. @)). 

Now let us include longer links. First we still limit 
ourselves to unfrustrated ansatz. An interesting ansatz 
is given by 



Ux 



XT 



Uy = 


= XT - 


U2x+y = 


= At2 


-x+2y = 


= -Xt' 


'a2x-y - 


= Xt^ 


Ux+2y = 


= -Xt^ 



rjT 
rjT 



(39) 



By definition, the ansatz is invariant under translation 
and parity x — > —x. After a 90° rotation, the ansatz is 



Ux 

Uy 

U2x+y 

U-x+2y 

U2x-y 

'ax+2y 



= -XT 
= -XT^ 

= -Xt' 
= +Xr^ 
= -Ar2 

= +Xt' 



rjT 
rjT 



(40) 



which is gauge equivalent to Eq. (39) under the gauge 
transformation GR^gii) = ir^. Thus the ansatz describe 
a spin liquid with translation, rotation, parity and the 
time reversal symmetries. The spinon spectrum is given 
by (see Fig. ||c) 



E± 
ei 

£2 

63 



±v/ei(fc)2+e2(fe)2-fe3(/c)2 

-2x(cos(fc:E) -I- cos(A;j;)) 

— 277(cos(fca;) — cos{ky)) 

— 2A[cos(2A:a; 4- ky) + cos(2A;2; — ky) 

— cos{kx — 2ky) — cos{kx + 2ky)] (41) 



Thus the spinous are gapless only at four k points 
(±7r/2,±7r/2). We also find that Pc^ and Pc^ do not 
commute, where the loops C3 = i^i + x^i + 2x^ 
i + 2x + y -^ i and C4 = i^i + y^i + 2y^ 
i + 2y-x^i. Thus the SU{2) flux Pc^ and P^ break 
the SU{2) gauge structure down to a Z2 gauge struc- 
ture. The spin liquid described by Eq. ( p9|) will be called 
the Z2-linear spin liquid. The low energy effective theory 
is described by massless Dirac fermions coupled to a Z2 
gauge field. Again the Z2 coupling is irrelevant and the 
spinous are free fermions at low energies. We have a true 
spin-charge separation. Accor ding to the classification 
scheme summarized in section [VB| , the above Z2-linear 
spin liquid is labeled by Z2A003n. 

Next let us discuss frustrated ansatz. A simple Z2 
spin liquid can be obtained from the following frustrated 
ansatz 



al^O, 







Ux ^XT^ + r/T^ 



--XT 



U-x+y =7T 



rjT 



(42) 



The ansatz has translation, rotation, parity, and the time 
reversal symmetries. When Oq 7^ 0, x 7^ ^V and XV 7^ Oi 
OqT does not commute with the loop operators. Thus the 
ansatz breaks the SU{2) gauge structure to a Z2 gauge 
structure. The spinon spectrum is given by (see Fig. ^a) 



€{k) =2x(cos(fca,) 
27(cos(fca, - 
A(fc) =2ry(cos(fc^) 



-A2(fe) 

f cos(fcy)) + al 
ky) +cos{kx - ky)) 
- cos{ky)) + al 



(43) 
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which is gapless only at four k points with a hnear dis- 
persion. Thus the spin hquid described by Eq. (|4^) is a 
Z2-hnear spin hquid, which has a true spin-charge separa- 
tion. The Z2-hnear spin hquid is described by the projec- 
tive symmet ry gr oup Z2A0032 or equivalently Z2A0013. 
(see section [VB| .) From the above two examples of Z2- 
linear spin liquids, we find that it is possible to obtain 
true spin-charge separation with massless Dirac points 
(or nodes) within a pure spin model without the charge 
fluctuations. We also find that there are more than one 
way to do it. 

A well known frustrated ansatz is the ansatz for the 
chiral spin liquid M 



XT 

Xr' 



Ux = 


= -XT 


Uy = 


= -XT 


x+y = 


. rr' 



-x+y 



-rjT 



= 



(44) 



The chiral spin liquid breaks the time reversal and parity 
symmetries. The SU{2) gauge structure is unbroken. [p8| 
The low energy effective theory is an SU{2) Chern- 
Simons theory (of level 1). The spinous are gaped and 
have a semionic statistics. R M The third interesting frus- 
trated ansatz is given in Ref. Bq, Eo] 



ua 



=u. 



V 



XT 
Ux+y =VT^ + ^T^ 

1 Ar2 



'^ — x+y 
,2.3 



a. 







=0, al ^ 



(45) 



This ansatz has translation, rotation, parity and the time 
reversal symmetries. The spinous are fully gaped and 
the SU{2) gauge structure is broken down to Z2 gauge 
structure. We may call such a state Z2-gapped spin liquid 
(it was called sRVB state in Ref. |3^, ^ ) . It is described 
by the projective symmetry group Z2Axa;0z. Both the 
chiral spin liquid and the Z2-gapped spin liquid have true 
spin-charge separation. 



IV. QUANTUM ORDERS IN SYMMETRIC 
SPIN LIQUIDS 

A. Quantum orders and projective symmetry 
groups 

We have seen that there can be many different spin 
liquids wi th th e same symmetries. The stability analysis 
in section 



VIII shows that many of those spin liquids oc- 



cupy a finite region in phase space and represent stable 
quantum phases. So here we are facing a similar situa- 
tion as in quantum Hall effect: there are many distinct 
quantum phases not separated by symmetries and order 
parameters. The quantum Hall liquids have finite energy 



gaps and are rigid states. The concept of topological 
order was introduced to describe the internal order of 
those rigid states. Here we can also use the topological 
order to describe the internal orders of rigid spin liquids. 
However, we also have many other stable quantum spin 
liquids that have gapless excitations. 

To describe internal orders in gapless quantum spin 
liquids (as well as gapped spin liquids), we have intro- 
duced a new concept - quantum order - that describes 
the internal orders in any quantum phases. The key point 
in introducing quantum orders is that quantum phases, 
in general, cannot be completely characterized by bro- 
ken symmetries and local order parameters. This point 
is illustrated by quantum Hall states and by the stable 
spin liquids constructed in this paper. However, to make 
the concept of quantum order useful, we need to find 
concrete mathematical characterizations the quantum or- 
ders. Since quantum orders are not described by symme- 
tries and order parameters, we need to find a completely 
new way to characterize them. Here we would like to 
propose to use Projective Symmetry Group to character- 
ize quantum (or topological) orders in quantum spin liq- 
uids. The projective symmetry group is motivated from 
the following observation. Although ansatz for different 
symmetric spin liquids all have the same symmetry, the 
ansatz are invariant under transformations followed by 
different gauge transformations. We can use those dif- 
ferent gauge transformations to distinguish different spin 
liquids with the same symmetry. In the following, we 
will introduce projective symmetry group in a general 
and formal setting. 

We know that to find quantum numbers that charac- 
terize a phase is to find the universal properties of the 
phase. For classical systems, we know that symmetry 
is a universal property of a phase and we can use sym- 
metry to characterize different classical phases. To find 
universal properties of quantum phases we need to find 
universal properties of many-body wave functions. This 
is too hard. Here we want to simplify the problem by lim- 
iting ourselves to a subclass of many-body wave functions 
which can be described by ansatz {uij, Aqt') via Eq. (M). 
Instead of looking for the universal properties of many- 
body wave functions, we try to find the universal prop- 
erties of ansatz {uij,aQT'-). Certainly, one may object 
that the universal properties of the ansatz (or the sub- 
class of wave functions) may not be the universal prop- 
erties of spin quantum phase. This is indeed the case for 
some ansatz. However, if the mean-field state described 
by ansatz (wjj, agT') is stable against fluctuations (ie the 
fluctuations around the mean-field state do not cause any 
infrared divergence), then the mean-field state faithfully 
describes a spin quantum state and the universal proper- 
ties of the ansatz will be the universal properties of the 
correspond spin quantum phase. This completes the link 
between the properties of ansatz and properties of phys- 
ical spin liquids. Motivated by the Landau's theory for 
classical orders, here we whould like to propose that the 
invariance group (or the "symmetry" group) of an ansatz 
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is a universal property of the ansatz. Such a group will 
be called the projective symmetry group (PSG). We will 
show that PSG can be used to characterize quantum or- 
ders in quantum spin liquids. 

Let us give a detailed definition of PSG. A PSG is a 
property of an ansatz. It is formed by all the transfor- 
mations that keep the ansatz unchanged. Each trans- 
formation (or each element in the PSG) can be written 
as a combination of a symmetry transformation U (such 
as translation) and a gauge transformation Gy. The in- 
variance of the ansatz under its PSG can be expressed 
as 

GuU{uij) =Uij 

U{uij) =uu(i)^u{j) 
Gu{uij) =Gu{i)uijG\j{j) 

Gu{i) £SU{2) (46) 

for each GuU e PSG. 

Every PSG contains a special subgroup, which will be 
called invariant gauge group (IGG). IGG (denoted by Q) 
for an ansatz is formed by all the gauge transformations 
that leave the ansatz unchanged: 

Q = {Wi\WiUijW] = Uij, Wi e 5[/(2)} (47) 

If we want to relate IGG to a symmetry transformation, 
then the associated transformation is simply an identity 
transformation. 

If IGG is non-trivial, then for a fixed symmetry trans- 
formation [/, there are can be many gauge transforma- 
tions Gjj that leave the ansatz unchanged. If GjjU is in 
the PSG of Uij, GGuU wiU also be in the PSG iff G e a. 
Thus for each symmetry transformation f7, the different 
choices of Gjj have a one to one correspondence with the 
elements in IGG. From the above definition, we see that 
the PSG, the IGG, and the symmetry group (SG) of an 
ansatz are related: 



SG = PSG I IGG 



(48) 



This relation tells us that a PSG is a projective repre- 
sentati on o r an extension of the symmetry group. pj5| (In 
section A 1 we will introduce a closely related but differ- 



ent definition of PSG. To distinguish the two definitions, 
we will call the PSG defined above invariant PSG and 



the PSG defined in section A 1 algebraic PSG.) 

Certainly the PSG's for two gauge equivalent ansatz 
Uij and W {i)uijW'^ [j) are related. From WGuU{uij) — 
W{uij), where W{uij) = W{i)uijW\j), we find 
WGuUW-^W{uij) = WGuW^^UW{uij) = W{uij), 
where Wu = UWU-^ is given by Wu{i) = W{U{i)). 
Thus if GjjU is in the PSG of ansatz Ujj, then 
{WGuWu)U is in the PSG of gauge transformed ansatz 
W {i)uijW'^ {j) . We see that the gauge transformation 
Gjj associated with the symmetry transformation U is 
changed in the following way 



after a gauge transformation W{i). 

Since PSG is a property of an ansatz, we can group 
all the ansatz sharing the same PSG together to form a 
class. We claim that such a class is formed by one or 
several universality classes that correspond to quantum 
phases. (A more detailed discussion of this important 

It is in this sense we say 



VIII E 



point IS given m section 

that quantum orders are characterized by PSG's. 

We know that a classical order can be described by 
its symmetry properties. Mathematically, we say that a 
classical order is characterized by its symmetry group. 
Using projective symmetry group to describe a quantum 
order, conceptually, is similar to using symmetry group 
to describe a classical order. The symmetry description 
of a classical order is very useful since it allows us to 
obtain many universal properties, such as the number of 
Nambu-Goldstone modes, without knowing the details of 
the system. Similarly, knowing the PSG of a quantum 
order also allows us to obtain low energy properties of 
a quantum system without knowing its details. As an 
example, we will discuss a particular kind of the low en- 
ergy fluctuations - the gauge fluctuations - in a quantum 
state. We will show that the low energy gauge fluctua- 
tions can be determined completely from the PSG. In fact 
the gauge group of the low energy gauge fluctuations is 
nothing but the IGG of the ansatz. 

To sec this, let us assume that, as an example, an 
IGG Q contains a U{1) subgroup which is formed by the 
following constant gauge transformations 



{Wi 



l^e [0,2^)} eg 



(50) 



Now we consider the following type of fluctuations 
around the mean- field solution Uij: Uij — Uije^""^^^ . 
Since Uij is invariant under the constant gauge trans- 



formation e*^"^ 
tion e*' 



a spatial dependent gauge transforma- 
transform the fluctuation a?- to a?- = 
a?o +di — 9i . This means that a?- and a?, label the same 
physical state and a|- correspond to gauge fluctuations. 
The energy of the fluctuations has a gauge invariance 
_E({a?-}) — -E({a?-}). We see that the mass term of the 
gauge field, (a?-)^: i^ not allowed and the U{1) gauge 
fluctuations described by a?- will appear at low energies. 
If the U{1) subgroup of Q is formed by spatial depen- 
dent gauge transformations 



{H^i = e^»"--|0e[O,27r),|ni|-l}cg, 



(51) 



Gu{i) -^ W{i)Gu{i)W\U{i)) 



(49) 



we can always use a SU{2) gauge transformation to ro- 
tate rii to the z direction on every site and reduce the 
problem to the one discussed above. Thus, regardless 
if the gauge transformations in IGG have spatial depen- 
dence or not, the gauge group for low energy gauge fluc- 
tuations is always given by G- 

We would like to remark that some times low energy 
gauge fluctuations not only appear near fc = 0, but 
also appear near some other k points. In this case, we 
will have several low energy gauge fields, one for each k 
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points. Examples of this phenomenon are given by some 
ansatz of SU{2) slave-boson theory discussed in section 
VI, which have an SU{2) x SU{2) gauge structures at 
low energies. We see that the low energy gauge struc- 
ture SU{2) X SU{2) can even be larger than the high en- 
ergy gauge structure SU{2). Even for this complicated 
case where low energy gauge fluctuations appear around 
different k points, IGG still correctly describes the low 
energy gauge structure of the corresponding ansatz. If 
IGG contains gauge transformations that are indepen- 
dent of spatial coordinates, then such transformations 
correspond to the gauge group for gapless gauge fluctua- 
tions near k = 0. If IGG contains gauge transformations 
that depend on spatial coordinates, then those transfor- 
mations correspond to the gauge group for gapless gauge 
fluctuations near non-zero k. Thus IGG gives us a unified 
treatment of all low energy gauge fluctuations, regardless 
their momenta. 

In this paper, we have used the terms Z2 spin liquids, 
[/(I) spin liquids, SU{2) spin liquids, and SU{2) x SU{2) 
spin liquids in many places. Now we can have a pre- 
cise definition of those low energy Z2, U{1), SU{2), and 
SU{2) X SU{2) gauge groups. Those low energy gauge 
groups are nothing but the IGG of the corresponding 
ansatz. They have nothing to do with the high energy 
gauge groups that appear in the SU{2), C/(l), or Z2 slave- 
boson approaches. We also used the terms Z2 gauge 
structure, [/(I) gauge structure, and SU{2) gauge struc- 
ture of a mean-field state. Their precise mathematical 
meaning is again the IGG of the corresponding ansatz. 
When we say a C/(l) gauge structure is broken down to 
a Z2 gauge structure, we mean that an ansatz is changed 
in such a way that its IGG is changed from U{1) to Z2 
group. 



B. Classification of symmetric Z2 spin liquids 

As an application of PSG characterization of quantum 
orders in spin liquids, we would like to classify the PSG's 
associated with translation transformations assuming the 
IGG Q = Z2. Such a classification leads to a classification 
of translation symmetric Z2 spin liquids. 

When Q = Z2, it contains two elements - gauge trans- 
formations Gi and 6*2: 



Gi(i) =rO, G2{i) - -r". 



(52) 



Let us assume that a Z2 spin liquid has a translation sym- 
metry. The PSG associated with the translation group 
is generated by four elements ±GxTx^ ±GyTx where 



-^xy^ij) — ^i—x.j — xj -^y\^ij) ^i — y,j' 



(53) 



Due to the translation symmetry of the ansatz, we can 
choose a gauge in which all the loop operators of the 
ansatz are translation invariant. That is Pci — Pc2 if 



the two loops Ci and C2 are related by a translation. 
We will call such a gauge uniform gauge. 

Under transformation GxT^, a loop operator Pc 
based at i transforms as Pc -^ Gx{i')PT^cG\.{i') — 
Gx{i')PcG\.{i') where i' — T^i is the base point of the 
translated loop Tx{C). We see that translation invariance 
of Pc in the uniform gauge requires 



Gxii)=±T", Gy{i) 



±rO. 



(54) 



since different loop operators based at the same base 
point do not commute for Z2 spin liquids. We note that 
the gauge transformations of form W{i) = ir*^ do not 
change the translation invariant property of the loop op- 
erators. Thus we can use such gauge transformations 
to further simplify Gx.y through Eq. (|49|). First we can 
choose a gauge to make 



Gy{i) 



T°. 



(55) 



We note that a gauge transformation satisfying W{i) ~ 
W{ix) does not change the condition Gy{i) = r°. We 
can use such kind of gauge transformations to make 

Gx{lx,iy=0) = T''. (56) 

Since the translations in x- and y-direction commute, 
Gx,y must satisfy (for any ansatz, Z2 or not Z2) 



GxTxGyTy(GxTx) (GyTy) — 

GxTxGyTyT^ Gx Ty Gy G G- 
That means 

Gx{i)Gy{i - x)G-Hi ~ y)Gy{i)-^ e Q 
For Z2 spin liquids, Eq. ( |5q ) reduces to 
G,(i)G-i(i-y)=+rO 
or 

G,(i)G-i(i-y)--rO 



(57) 



(58) 



(59) 



(60) 



When combined with Eq. (55) and Eq. (b6|), we find that 



there are only two gauge inequivalent extensions of the 
translation group when IGG \s, Q = Z2- The two PSG's 
are given by 



and 



Gx{i) 



Gx{i) ^{-r^T\ 



Gyii) =r° 



Gy{i) 



(61) 



(62) 



Thus, under PSG classification, there are only two types 
of Z2 spin liquids if they have only the translation sym- 
metry and no other symmetries. The ansatz that satisfy 
Eq. (|6l] ) have a form 



^t,x-\-m 



(63) 
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and the ones that satisfy Eq. (62) have a form 



={-7' 



(64) 



calc ulated through the algebraic relations listed in section 
A 1. The algebraic PSG's are different from the invariant 



Through the above example, we see that PSG is a very 
powerful tool. It can lead to a complete classification of 
(mean-field) spin liquids with prescribed symmetries and 
low energy gauge structures. 

In the above, we have studied Z2 spin liquids which 
have only the translation symmetry and no other sym- 
metries. We find there are only two types of such spin 
liquids. However, if spin liquids have more symmetries, 
then they can have much more types. In the appendix 
A, we will give a classification of symmetric Z2 spin liq- 
uids using PSG. Here we use the term symmetric spin 
liquid to refer to a spin liquid with the translation sym- 
metry Tx^yi the time reversal symmetry T: Uij — > —Uij, 
and the three parity symmetries Px'. (zrc, iy) -^ {—ix, iy), 

-fy. y^x^'^y) ^ \^a;5 ^y)i 9-110- ^xy- \^X7^y) ^ y^y^'^xj- 

The three parity symmetries also imply the 90° ro- 
tation symmetry. In the appendix A, we find that 
there are 272 different extensions of the symmetry group 
{Tx,y,Px,y,xy,T} if IGG g = Z2. Those PSG's are gen- 
erated by iGxTx,GyTy,GTT,Gp^Px,Gp^Py,Gp^^Pxy). 

The PSG's can be divided into two classes. The first 
class is given by 

Gx{i) =rO, 

GP. (i) ^VxpxVx'pygP. 

Gp^yii) ^gP^y 

and the second class by 

G.(i) H-r-r^ 

GP^ (i) ^VxpxVx'pygP. 

Gp.y{i)={-r-'^9P.y 

Here the three 77's can independently take two value s ±1 . 
g^s ha ve 17 different choices which are given by Eq. ( A39) 
- Eq. ( [A55| ) in the appendix A. Thus there are 2 x 17 x 
2^ = 272 different PSG's. They can potentially lead to 
272 different types of symmetric Z2 spin liquids on 2D 
square lattice. 

To label the 272 PSG's, we propose the following 
scheme: 

Z2A{gpx)r,^^S9py)ru„y9pxy{9t)r,t, (67) 

Z2B{gpx)rj^^^ {9py)v.py9pxy{gt)nt ■ (68) 



Gy{i) 


=t" 


Gpyii) 


~'lxpy'lxpxyPy 


Grii) 


^Vigr (65) 


Gyii) 


=t" 


Gpyii) 


~~ 'Ixpy 'ixpxljPy 


Grii) 


^Vlgr (66) 



The label Z2A... correspond to the case Eq. (|65|), and 
the label Z2B... correspond to the case Eq. (pq). A typi- 
cal label will looks like Z2Arj|: r^r^^r^. We will also use 
an abbreviated notation. An abbreviated notation is ob- 
tained by replacing {t^ , t^ , t^ , t^) or (rij^ , r^ , r^ , r^ ) by 
(0,1,2,3) and {t^ , t}_ , t^ , t^) by {n,x,y,z). For exam- 
ple, Z2At|t°t^^t^ can be abbreviated as Z2Aln(12)2. 
Those 272 different Z2 PSG's, strictly speaking, are the 
so called algebraic PSG's. The algebraic PSG's are de- 
fined as extensions of the symmetry group. They can be 



PSG's which are defined as a collection of all transfor- 
mations that leave an ansatz Uij invariant. Although an 
invariant PSG must be an algebraic PSG, an algebraic 
PSG may not be an invariant PSG. This is because cer- 
tain algebraic PSG's have the following properties: any 
ansatz Uij that is invariant under an algebraic PSG may 
actually be invariant under a larger PSG. In this case 
the original algebraic PSG cannot be an invariant PSG 
of the ansatz. The invariant PSG of the ansatz is re- 
ally given by the larger PSG. If we limit ourselves to the 
spin liquids constructed through the ansatz Uij , then we 
should drop the algebraic PSG's are not invariant PSG's. 
This is because those algebraic PSG's do not characterize 
mean-field spin liquids. 

We find that among the 272 algebraic Z2 PSG's, at 
least 76 of them are not invariant PSG's. Thus the 272 al- 
gebraic Z2 PSG's can at most lead to 196 possible Z2 spin 
liquids. Since some of the mean-field spin liquid states 
may not survive the quantum fluctuations, the number of 
physical Z2 spin liquids is even smaller. However, for the 
physical spin liquids that can be obtained through the 
mean-field states, the PSG's do offer a characterization 
of the quantum orders in those spin liquids. 



C. Classification of symmetric U{1) and SU{2) spin 
liquids 

In addition to the Z2 symmetric spin liquids studied 
above, there can be symmetric spin liquids whose low 
energy gauge structure is U(l) or SU{2). Such f/(l) 
and SU{2) symmetric spin liquids (at mean-field level) 
are classified by U{1) and SU{2) symmetric PSG's. The 
U{1) and SU{2) symmetric PSG's are calculated in the 
appendix A. In the following we just summarize the re- 
sults. 

We find that the PSG's that characterize mean-field 
symmetric U{1) spin liquids can be divided into four 
types: UlA, UIB, UlC and Ul™. There are 24 type 
UlA PSG's: 

Gx ==.93(^3;), Gy = gaiOy), 

GP^ =Vllxg3{0px), Gp^ = Vypx93{0py) 

Gp^y —gsidpxy), g'i{Spxy)iT 

Gt =?7?33(^t)U*=-i, 4g3{et)iT^ (69) 



and 



Gx —93(9 x), Gy — .93 (^y), 



Gp^ ^Vllxga 



pxHiK^pxl 



Gpy =- Vxixg^iOpy^T^ 



Jpxyl 



Gp^y —93i0pxy), 53 

Gt =vi93idt)\vt=-u 493{0t)iT^ 



where 



ga{e) 



(70) 



(71) 
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We will use VlAarf^^^bn^^^cdrn to label the 24 PSG's. 
a, b, c, d are associated with Gp^ , Gp^ , Gp^^ , Gt re- 
spectively. They are equal to r^ if the corresponding G 
contains a r^ and equal to t'^ otherwise. A typical nota- 
tion looks like UIAtIt^t^tI which can be abbreviated 
as UlAxlOx. 

There are also 24 type UIB PSG's: 

G,={-yyg3{e,), G,=53(e,), 

GP. =Vy'pcc93{0px), Gp^ = Vypx93{9py) 

i-y-'^Gp^^ ^g3{0p^y), gz{ep^y)iT^ 

Gt =vi93{et)W=-i, 493{0tyr' (72) 



and 



G. =(-^53(^x), Gy^g^iOy), 
Gp^ =Vxpx93{0px)iT^ , Gp^ ^ r]l^px93{0py)iT^ 
{-y-'^Gp^^ ^gsiOp^y), g3{9pxy)iT^ 

Gt ^493{0t)U=-i, 493{et)^r^ (73) 

We will use UlBa,,^^^6,,^p^crf,,j to label the 24 PSG's. 
The 60 type Ulc'pSG's are given by 

Gx =93{Sx)-iT^, Gy = g3{0y)iT^, 

GP^ ^VxpxVllx93{0px), Gp^ = VypxVllxgsiOpy) 



xpx lypxjo v"/^j. / ; -* y lypx ixpx-:^ 

7T 
Gp^y ^Vpxy93{vlxv^ + ^pxy), 

Gt =vi93{0t)\nt = -l, Vpxy93{0t)lT^ 



(74) 



Gx ^g3{0x)iT^, Gy = g3{9y)iT'^, 

Gp, =Vl^px93iOpx)iT\ Gp^ = Vx'pxVlxy93iOpy)iT^ 

Gp^y ^Vply93{llpxy^ +dpxy), 

Gt ^493{0t)K^-i, v;iyVi93{et)tr' (75) 



Gx =g3{0x)-iT'^ , Gy = g3{0y)iT'^, 

Gp^ =Vl'pxVy'px93{0px), Gp^ = VllxVxpx93{0py) 

Gp^y =g3{0pxy)iT^ 

Gt =4g3{et%,=-i (76) 



Gx =g3{0x)iT^, Gy = g3{9y)iT^, 



Gp^ ^rfxpxVllx93{ 



GPy ^ VypxVZx93{0py) 



Gp^y ^93{rfp^y^+0pxy)iT^ 
Gt =VpxyVi93{0tW 



Gx 
Gp. 



=93{0x)iT'^, Gy = g3{dy)iT'^ , 
' Gp 



^Vxpx93[t^px)lT 

7T 
Gp^y ^93{rilxy- + Opxy)iT^ 



Ixpx 'Ipxyy^x 



Gt =vlg3 



pxy ^ I ^px 
't)l>)t = -l: 



vWply93iOtW (78) 

which will be labeled by VlCar,^^^br^^p^Cr,p^^d,^^. 

The type Ul™ PSG's have not been classified. How- 
ever, we do know that for each rational number m/n e 
(0, 1), there exist at least one mean-field symmetric spin 
liquid, which is described by the ansatz 



^i^i+x 



XT 



U. 



i,i+y 



I '"■" • \ 3 

X93[ lx)T 



(79) 



It has nm/n flux per plaquette. Thus there are infinite 
many type Ul™ spin liquids. 

We would like to point out that the above 108 
U1[A,B,C] PSG's are algebraic PSG's. They are only a 
subset of all possible algebraic U{1) PSG's. However, 
they do contain all the invariant U{1) PSG's of type 
UlA, UIB and UlC. We find 46 of the 108 PSG's are 
also invariant PSG's. Thus there are 46 different mean- 
field U{1) spin liquids of type UlA, UIB an d UlC Thei r 
ansatz and labels are given by Eq. ( AllOJ ), Eq. (Alll), 
Eq. ( |A123| ), Eq. ( |AT2^ ), and Eq. ( |A152D - Eq. ( |A169D . 

To classify symmetric SU{2) spin liquids, we find 8 
different SU{2) PSG's which are given by 

Gx{i)^gx, Gy{i) = gy 

Gp^ (i) ^rfxpxVxpygp. : Gp^ (i) = VxpyVl^xgPy 

Gp.yii) ^gp.y, GTii) = i-ygT (80) 

and 

Gx{i) ={-)'" 9x, Gy{i)^gy 

GP^ (i) ^VxpxVx'pygP. , Gp^ (i) = VxpyVllxgPy 

Gp.yii)H-y-'^gp.y, GT{i) = {-ygT (81) 

where g's are in SU{2). We would like to use the following 
two notations 



SU2At" t" 

'Ixpx 'Ixpy 

SU2Bt° t° 

'Ixpx 'Ixpy 



(82) 



to denote the above 8 PSG's. SU2At^^^^ r^^^_^ is for 
Eq. ^) and SU2BTO^^^rO^^^ for Eq. (|8|). We find only 
4 of the 8 SU{2) PSG^s^ SU2A[nO, On] and SU2B[nO, On], 
leads to SU{2) symmetric spin liquids. The SU2AnO 
state is the uniform RVB state and the SU2BnO state 
is the TT-flux state. The other two SU{2) spin liquids are 
given by SU2A0n: 

'U'i^i+2x+y = + IXT 



^i.i — 2x-\-y 



IXT 



(77) 



Ui,i-x+2y = + iXT° 



(83) 
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and SU2B0n: 



^i,i-\-2x-\-y 
Ui,i-~2x-\-y 
^i,i-\-x-\-2y 
Ui,i—x-\-2y 



y = + iX^° 

= + iXT° 



(84) 



The above results give us a classification of symmet- 
ric C/(l) and 51/(2) spin liquids at mean-field level. If 
a mean-field state is stable against fluctuations, it will 
correspond to a physical U{1) or SU{2) symmetric spin 
liquids. In this way the C/(l) and the SU{2) PSG's also 
provide an description of some physical spin liquids. 



V. CONTINUOUS TRANSITIONS AND 

SPINON SPECTRA IN SYMMETRIC SPIN 

LIQUIDS 

Continuous phase transitions witliout symmetry 
breaking 



After classifying mean- field symmetric spin liquids, we 
would like to know how those symmetric spin liquids 
are related to each other. In particular, we would like 
to know which spin liquids can change into each other 
through a continuous phase transition. This problem 
is studied in detail in appendix B, where we study the 
symmetric spin liquids in the neighborhood of some im- 
portant symmetric spin liquids. After lengthy calcula- 
tions, we found all the mean-field symmetric spin liq- 
uids around the Z2-linear state Z2A001n in Eq. (39), 
the t/(l)-hnear state UlCnOln in Eq. (H), the SUjT)- 
gapless state SU2AnO in Eq. (|^, and the 5t/(2)-linear 
state SU2BnO in Eq. (|3l|). Those ansatz are given by 
Eq. ( ^ f or the Za-linear state, by Eq. (H), Eq. (|B24[ ), 
Eq. ( tB25| ), Eq. ( p2^, a nd E g. (|B28|) for the C/(l)- 
linear state, by Eq. (|B3lD , Eq. (|B44D - Eq. (p49| ), and 



Eq. ( p9^ - E q. ( |B10S|) fo r the 5 L/(2)-g apless stat e, and 
by Eq. (|B111| ), Eq. ( |B117D - Eq. (|B122| ), and Eq. ( |B139| ) 
- Eq. (B154) for the SU{2)-\mea,r state. According to 
the above results, we find that, at the mean-field level, 
the ?7(l)-linear spin hquid UlCnOln can continuously 
change into 8 different Z2 spin liquids, the 5'C/(2)-gapless 
spin liquid SU2AnO can continuously change into 12 U{1) 
spin liquids and 52 Z2 spin liquids, and the 5C/(2)-lincar 
spin liquid SU2BnO can continuously change into 12 U{1) 
spin liquids and 58 Z2 spin liquids. 

We would like to stress that the above results on the 
continuous transitions are valid only at mean-field level. 
Some of the mean-field results survive the quantum fiuc- 
tuations while others do not. One need to do a case by 
case study to see which mean-field results can be valid 
beyond the mean-field theory. In Ref. [^, a mean-field 
transition between a SU{2) x S'C/(2)-linear spin liquid 
and a Z2-gapped spin liquid was studied. In particular 
the effects of quantum fiuctuations were discussed. 



We would also like to point out that all the above 
spin liquids have the same symmetry. Thus the contin- 
uous transitions between them, if exist, represent a new 
class of continuous transitions which do not change any 
symmetries. [pGl 



B. Symmetric spin liquids around the U{1) 
spin Hquid UlCnOln 



■linear 



The S'[/(2)-hnear state SU2BnO (the 7r-fiux state), the 
J7(l)-linear state UlCriOln (the staggered-flux/d-wave 
state), and the S'{/(2)-gapless state SU2AnO (the uniform 
RVB state), are closely related to high Tc superconduc- 
tors. They reproduce the observed electron spectra func- 
tion for undoped, underdoped, and overdoped samples 
respectively. However, theoretically, those spin liquids 
are unstable at low energies due to the f/(l) or SU{2) 
gauge fluctuations. Those states may change into more 
stable spin liquids in their neighborhood. In the next a 
few subsections, we are going to study those more sta- 
ble spin liquids. Since there are still many different spin 
liquids involved, we will only present some simplified re- 
sults by limiting the length of non-zero links. Those spin 
liquids with short links should be more stable for simple 
spin Hamiltonians. The length of a link between i and j 
is defined as \ix — jx\-'r \iy — jy\- By studying the spinon 
dispersion in those mean-field states, we can understand 
some basic physical properties of those spin liquids, such 
as their stability against the gauge fluctuations and the 
qualitative behaviors of spin correlations which can be 
measured by neutron scattering. Those results allow us 
to identify them, if those spin liquids exist in certain sam- 
ples or appear in numerical calculations. We would like 
to point out that we will only study symmetric spin liq- 
uids here. The above three unstable spin liquids may also 
change into some other states that break certain symme- 
tries. Such symmetry breaking transitions actually have 
been observed in high Tc superconductors (such as the 
transitions to antiferromagnetic state, d-wave supercon- 
ducting state, and stripe state). 

First, let us consider the spin liquids around the 
J7(l)-linear state UlCnOln. In the neighborhood of the 
UlCnOln ansatz Eq. (p2|), there are 8 classes of symmet- 
ric ansatz Eq. ( p24| ), Eq. ( p25| ) Eq. ( p27| ), and Eq. ( p2^ ) 
that break the U(i) gauge structure down to a Z2 gauge 
structure. The first one is labeled by Z2A0013 and takes 
the following form 



Ui,i+y =XT^ + m"^ 



u 



t,i+x+y 



= + llT 



^i.i—x-\-y """^ ' Tl"^ 

Ui,i+2x =12T^ + A2T^ 
Ui,i+2y =72T^ - A2T^ 



4 7^ 0' %' = 



(85) 
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It has the same quantum order as that in the ansatz 
Eq. (H). The label Z2A0013 tells us the PSG that char- 
acterizes the spin liquid. 

The second ansatz is labeled by Z2Azzl3: 



,i+y =XT^ 



rjT 



u. 



'i,i+x+y 



^1,1— £c+y 



-7ir 



f^i,i+2x —Uia+2y 



= 



1,2,3 
^0 



= 



(86) 



The third one is labeled by Z2A001n (or equivalently 
Z2A003n): 



<=0 

^i,i-\-y X^ 
y'i.i-x+2ij = ~ At 



r]T 



,,i — x-\-2y 
^i,i-\-2x — % 



--\t' 



UiA+x+2y — — At 



(87) 



Such a spin liquid has the same quantum order as 
Eq. (pa). The fourth one is labeled by Z2Azzln: 



«o=0 



'XT + VT 



^i,i+y X^ 



rjT 



Ui,i+2x+y =Xl'^ + Vl'T + 



^i.i+2x — y Xl^ 






Ui,i+x+2y ^XlT^ - Tjit"^ - At^ (88) 

The above four ansatz have translation invariance. The 
next four Z2 ansatz do not have translation invariance. 
(But they still describe translation symmetric spin liquids 
after the projection.) Those Z2 spin liquids are Z2B0013: 

Ui^i+x =XT^ - VT^ 

Ui,i+y ={-)'■'' iXT^ +7]t'^) 

Ui,i+2x = - 72T^ + A2T^ 

1 \ 2 

"■" >2S = - 72T - A2T 



^z^z+zy 



al^O, 



Z2Bzzl3: 



Ui 



--XT 



2,3 



rjT 



= 0, 



(89) 



.1 



Ui4+2x+2y 



-2x+2y =1\T 



71 r 
1 



1,2,3 



0, 



(90) 



^^4#j' V 




FIG. 1: Contour plot of the spinon dispersion E+{k) as a 
function of {kx/2n,ky/2TT) for the Z2-linear spin liquids, (a) 
is for the Z2A0013 state in Eq. {^, (b) for the Z2A2zl3 state 
in Eq. (||), (c) for the Z2A001n state in Eq. (|^), (d) for the 
Z2Azzln state in Eq. 



Z2B001n: 







^iX+x 


1 , 2 
=XT +riT 






""iX+y 


= (- 


-y^ixr' 




"i 


i+2x+y 


= (- 


-Y-Xt' 




"i 


i-x+2y 


= - 


-Ar3 




"i 


i+2x—y 


= (- 


-)^-Ar3 




"i 


i+x+2y 


=0, 


-Ar^ 


and Z2Bzz 


In: 










Ux ^XT^ 


+ Vr' 



(91) 



UyH-T'^iXT^-Vr^) 



^2x+y 



'■-x+2y 



U 



2x-y 



U 



x+2y 



--XiT^ - m-r'^ + Ar^ 
= {-y-{XiT'+mr^-^r^) 
=Xir^ - mr^ - Ar^ 



a^=0. 



(92) 



The spinous are gapless at four isolated points with a 
linear dispersion for the first four Z2 spin liquids Eq. (|8q), 
Eq. (|^), Eq. 1^), and Eq. (g). (See Fig. ^ Therefore 
the four ansatz describe symmetric Z2-linear spin liquids. 
The single spinon dispersion for the second Z2 spin liq- 
uid Z2Azzl3 is quite interesting. It has the 90° rotation 
symmetry around k — (0, tt) and the parity symmetry 
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FIG. 2: Contour plot of the spinon dispersion 

niin(_Bi(fe), _B2(fe)) as a function of {kx/2Tr, ky/2Ti) for the Z2- 
linear states, (a) is for the Z2B0013 state in Eq. 1^, (b) for 
the Z2Bz2l3 state in Eq. @, (c) for the Z2B001n state in 
Eq. ^, (d) for the Z2Bz2ln state in Eq. (||). 



about k — (0,0). One very important thing to notice 
is that the spinon dispersions for the four Z2-hnear spin 
hquids, Eq. (||), Eq. ®, Eq. (|87|), and Eq. (|8|) have 
some quahtative differences between them. Those differ- 
ences can b e use d to physicahy measure quantum orders 
(see section VII). 



Next let us consider the ansatz Z2B0013 in Eq. (89) 



The spinon spectrum for ansatz Eq. (p9[) is determined 

by 



- 2xcos{ky)Ti + 277Cos(/i;j^)r3 + Ar4 
where kx G (0, tt), ky e (— 7r,7r) and 



(93) 



r2 =T^(g)T\ r3=T^(g>T\ 

ri=T^(g)T°. 



(94) 



assuming 71^2 = -^2 = 0. The four bands of spinon 
dispersion have a form zLEi{k), ztE2{k). We find the 
spinon spectrum vanishes at 8 isolated points near fc — 
(7r/2,±7r/2). (See Fig. |a.) Thus the state Z2B0013 is a 
Z2-linear spin liquid. 

Knowing the translation symmetry of the above Z2- 
linear spin liquid, it seems strange to find that the spinon 
spectrum is defined only on half of the lattice Brillouin 
zone. However, this is not inconsistent with translation 
symmetry since the single spinon excitation is not phys- 
ical. Only two-spinon excitations correspond to physical 
excitations and their spectrum should be defined on the 
full Brillouin zone. Now the problem is that how to ob- 
tain two-spinon spectrum defined on the full Brillouin 



zone from the single-spinon spectrum defined on half of 
the Brillouin zone. Let \k, 1) and \k, 2) be the two eigen- 
states of single spinon with positive energies Ei{k) and 
E2{k) (here k^ e (-7r/2,7r/2) and ky G (-7r,7r)). The 
translation by x (followed by a gauge transformation) 
change |fc, 1) and |fc, 2) to the other two eigenstates with 
the same energies: 



|fc,l) 
|fc,2) 



\k 
\k 



Ti-y,!) 
^y,2) 



(95) 



Now we see that momentum and the energy of two-spinon 
states |fci,Q;i)|fe2,a2) ± |fci + 7ry,ai)|fc2 + '^y,oi2} are 
given by 



E2spinon — -E'Qi(fei) + Ea2{k2) 

k =fci + fe2, ki + k2 + TTX 



(96) 

Eq. ( |9q ) allows us to construct two-spinon spectrum from 
single-spinon spectrum. 

Now let us consider the ansatz Z2Bzzl3 in Eq. (|9C 
The spinon spectrum for ansatz Eq. (p3) is determined 
by 



(97) 

2ky)T4 



H ^- 2xcos(A:a;)ro - 277COs(fca;)r2 

— 2xcos(A:j^)ri -I- 2r7COs(fcy)r3 

— 271 cos(2A:j; -I- 2fcj,)r4 + 271 cos(2fca; 
where k^ G (0,7r), ky e (— 7r,7r) and 

T2^T^(g>T\ T3=T^(g)T\ 

Ti=T^(g>T". (98) 

We find the spinon spectrum to vanish at 2 isolated points 
fc = (7r/2,±7r/2). (See Fig. |b.) The state Z2Bzzl3 is a 
Z2-linear spin liquid. 

The spinon spectrum for the ansatz Z2B001rt in 
Eq. (pi| ) is determined by 

H = - 2xcos{kx)To - 2-qcos{kx)T2 

— 2xcos{ky)Ti + 2ricos{ky)r3 

+ 2X{cos{kx + 2ky) + cos{-kx + 2ky))Vi 

- 2\{cos{2kx + ky) + cos(2A;;j; - ky))T5 (99) 
where k^ G (0,7r), ky G (— 7r,7r) and 



Fo 
F2 



Fi =r' . 

T3 =T^ ' 

Fs=r^ 



(100) 



The spinon spectrum vanishes at 2 isolated points k = 
(7r/2,±7r/2). (See Fig. |c.) The state Z2B001n is also a 
Z2-linear spin liquid. 

The spinon spectrum for the ansatz Z2Bzzln in 
Eq. (P^ ) can be obtained from 

H = - 2xcos(/i;j;)Fo - 277Cos(fc:r)F2 

— 2xcos(/i;y)Fi + 277COs(fcy)F3 

-I- 2A(cos(fc;j, + 2ky) — cos{—kx + 2ky))r4 

- 2A(cos(2fc^ + ky) - cos(2fc^ - ky))r5 (101) 
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where k^ € (0, tt), ky G (— tt, tt) and 



rn=ri 



T,=t' 






(102) 

We have also assumed that Xi = '?i = 0- The spinon 
spectrum vanishes at 2 isolated points k = (7r/2,±7r/2). 
(See Fig. Qd.) The state Z2Bzzln is again a Z2-linear 
spin liquid. 



C. Symmetric spin liquids around the 
5't/(2)-gapless spin liquid SU2A7iO 

There are many types of symmetric ansatz in the 
neighborhood of the 5'f7(2)-gapless state Eq. (pO|). Let us 
first consider the 12 classes of symmetr ic U ( 1 ) spin l iquid s 
around the 5'C/(2)-gapless state Eq. (p4^) - Eq. (|B49| ). 
Here we just present the simple cases where Uij are non- 
zero only for links with length < 2. Among the 12 classes 
of symmetric ansatz, We find that 5 classes actually give 
us the S't/(2)-gapless spin liquid when the link length is 
< 2. The other 7 symmetric U{1) spin liquids are given 
bellow. 



From Eq. (B50|) we get 



1,2,3 



=0 



rjT 



H.i+v ^XT 



rjT 



(jTx — Cry — T 



Gp^ =Gpy ~ T 
Gt =(-)V" 







(103) 

In the above, we have also listed the gauge transforma- 
tions Gx V, Gp. p p and Gt associated translation, 
parity and time reversal transformations. Those gauge 
transformations define the PSG that characterizes the 



C/(l) spin liquid. In section IV C, we have introduced a 
notation UlCnOlri to label the PSG and its associated 
ansatz. In the following, we will list ansatz together with 
their labels a nd th e associated gauge transformations. 
From Eq. (^5^) we get UlCnOOa; state 

Ui.i+y =XT^ 



^iA-\-x X^ 

3 
^i,i-]-x-\-y ^1^ 

Ui,i+2x =mT^ 

,3 



^i.i—x-\-y ^l'^ 



a, 







=%, 



Gx — Ctii — T 



Gp = 
UlCnOla; state 



= XT 



Ui^i+2x 
1,2,3 



'72T 



Ui,i+2y 


=?72T 


1,2 


=0 


GPx 


=Gpy 


Gt 


=zt2 


Uia+y = 


-Xr' 


H,i+2y = 


^V2T^ 



Gx Gy T 



Gp. 



-IT 



Gp^ =Gp^ =t", 
Gt =iT^ 



(105) 
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FIG. 3: Contour plot of the spinon dispersion E+{k) 
as a function of (kx/2 n, fc„ /27r) for (a) the [/(l)-linear 
state UlCnOOa: in Eq. (p^ ), and (b) the ;7(l)-quadratic 
state UlCilOa; in Eq. (106). In the f/(l)-quadratic state, 
the spinon energy vanishes as Afe^ near two points k = 
(7r,0),(0,7r). 



and UlCxlOx state 

^i.i-\-x X^ 
^i^i-\-x+y V^ 



1,2,3 



=0 



Gx — Cry — T 



Gf 



^i,i-\-y X^ 
3 
^i,i — x-\-y V^ 



Gp^^Gp^ ^it\ 
Gt =iT^ 



(106) 



From Eq. ( p4^ ) we get UlAOOOl state 



Ui,i+x =iXT'^ 
i,i+x-\-y ^l'^ 



Ui,i+2x ^mT 

1,2,3 



Ui,i+y =iXT 
3 
^i,i — x-\-y ^1^ 

Ui,i+2y =V2T^ 



a 



=0 



Gx =Gy = T°, {-T'^Gp^ ={-y^Gp^ = rO, 
Gpx, =T°, Gt =i{-yr' (107) 



and UlAOOll state 
Ui,i+x =iXT° 



UiA+2x 
1,2,3 



a 



=0 



Ui,i+2y =r]2T^ 



'0 

Ga; =Cry = T , 



Gp. 



=?r 



i-y^Gp^={-y-Gp^=T°, 

Gt =i(-)Vi (108) 



(104) state 



From Eq. ( p4^ ) we get, for Gp^^ = gsidpxy), UlAxlOx 



^i.i-\-x ^X^ 



1,2,3 



^i.i-\-y ^X^ 
3 



a. 







=0 



Gx =Gy = r", {-y^Gp^ =(-)'^Gp„ = iri, 



Gp „ =tO 



T , 



Gt =*(-)Vi 



(109) 
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FIG. 4: Contour plot of the spinon dispersion i5+(fc) as a 
function of {kx/2Ti,ky/2-K) f or th e (/(l)-gapless states, (a) is 
for the UlCnOla; state Eq. (p^), (b) for the UlAOOOl state 
Eq. ( |i07| ), (c) for the UlAOOll state Eq. (|io^), and (d) for 
the UlAa;10a; state Eq. (p^. 



Eq. (|10|) is the UlCnOln [/(l)-linear state (the stag- 
gered flux state) studied in the last section. After exam- 
ining the spinon dispersion, we find that the UlCnOOa; 
state in Eq. ( |104| ) can be a [/(l)-linear or a C/(l)-gapped 
state depending on the value of aj]. If it is a U{1)- 
linear state, it will have 8 isolated Fer mi points (see Fig. 
|a). The UlCnOlx state in Eq. (|l0|) is a C/(l)-gapless 
state (see Fig. |a). The UlCxlOx state in Eq. ( 10^ ) 
has two Fermi points at fci = (tt, 0) and fc2 = (0,7r). 
(see Fig. 0b). However, the spinon energy has a form 
E{k.) (X {k — fei,2)^ near fei and fc2- Thus we call the 



UlCxlOa; spin hquid Eq. (106) a C/(l)-quadratic state. 
The UlAOOOl state in Eq. (fo^), the UlAOOll state in 



Eq. ( [TOSD , and the UlAxlOx state in Eq. (|lO|) are C/(l)- 
gapless states (see Fig. |[). Again the spinon dispersions 
for the C/(l) spin liquids have some qualitative differences 
between each other, which can be used to detect different 
quantum orders in those C/(l) spin liquids. 

We next consider the 52 classes of symmetri c Z2 
spin liquids around the S'C/(2)-gapless state Eq. (B92) 
- Eq. (BIOS). Here we just present the simplest case 



where Uij are non-zero only for links with length < 1. 
We find that 48 out of 52 classes of ansatz describe U(l) 
or SU{2) spin liquids when the link's length is < 1. In 
the following we discuss the 4 remaining Z2 ansatz. 

W e ob tain one Z2 spin liquid Z2Aa:2(12) n from 
Eq. ( pQ^ . It is described by Eq. (||). From Eq. ( |B100[ ) 



we obtain a Z2 spin liquid Z2A0 013. It is described by 
Eq. (H) or Eq. (||). From Eq. ( |B102| ), we obtain a Z2 
spin liquid Z2Bj/l(12)n (note Z2Byl(12)n is gauge equiv- 
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FIG. 5: Contour plot of the spinon dispersion E+{k) as a 
function of (kx/2-K,ky/2Ti) for the Z2 spin liquids, (a) is for 
^2-gapless state Z2Aa;2(12)n in Eg. (p3[ ), and (b) is for Z2- 
quadratic state Z2Ba;2(12)n in Eq. (110). Despite the lack of 
rotation and parity symmetries in the single spinon dispersion 
in (a), the two-spinon spectrum does have those symmetries. 



alent to Z2Ba;2(12)n ): 



^,^+V 
1,2,3 



=0 



(110) 



From Eq. (|B108|), we obtain a Z2 spin hquid Z2B0013, 
which is described by Eq. (89). 

The ansatz Z2Ba::2(12)n in Eq. (|llC| ) is a n ew Z 2 spin 
liquid. The spinon spectrum for ansatz Eq. ( |110D is de- 
termined by 



H = - 2xsin(fc^)ro + 277Cos(A:^)r2 
- 2xsin(fcy)ri -I- 277C0s(/cj^)r3 

where kx € (— 7r/2,7r/2), ky £ (— tt, tt) and 



To 
El 



Tz 
Tg 



(111) 



(112) 



The spinon spectrum can be calculated exactly and 
its four branches take a form zLEi{k) and ±E2{k). 
The spinon energy vanishes at two isolated points fc = 
(0, 0), (0, tt). Near fc = the low energy spectrum is given 
by (see Fig. ^d) 



E: 



±^^ V(X2 + V')Hkl - fc2)2 + 4^4fc2fc2 (113) 



It is interesting to see that the energy does not vanish 
linearly as fc ^ 0, instead it vanishes like fc . 

We find that the loop operators for the following loops 
i^i + x^i + x + y^i + y^'i and i ^ i + y ^ 
i — x + y^i — x^i do not commute as long as both 
X and 77 are non-zero. Thus the spin liquid described 
by Eq. ( |110| ) indeed has a Z2 gauge structure. We will 
call such a state Z2-quadratic spin liquid to stress the 
E oc k dispersion. Such a state cannot be constructed 
from translation invariant ansatz, and it is the reason 
why we missed this state in the last section. The two- 
spinon spectrum is still related to the one-spinon spec- 
trum through Eq. WW. 
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D. Symmetric spin liquids around the SU{2)-linear 
spin liquid SU2B71O 

Last, we consider symmetric states in the neighbor- 
hood of the SU {2)-\mea.r state Eq. (|l|). The PSG's for 
those sym metric states can be obtained through the map- 
ping Eq. ( |A112 ) from the PSG's of symmetric spin hq- 
uids around the S'C/(2)-gapless spin hquid. Here we will 
only consider the 12 classes of symmetric U{1) spin liq- 
uids aroun d the 5t/(2)-linear state given by Eq. (B117) - 
Eq. ( B122| ). We will just present the simple cases where 



^ZJ 



are non-zero only for links with length < 2. We find 
that 7 of 12 classes of ansatz actually give us SU{2)- 
gapless spin liquids when the link length is < 2. Thus we 
only obtain the following 5 symmetric f/(l) spin liquids. 



From Eq. (B131) we get UlCnOln ansatz 
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FIG. 6: Contour plot of the spinon dispersion E+{k) as 
a function of (k x/2T r, ky/2Ti) for (a) the L''(l)-linear state 
UlCnOa;l in Eq. (pl5|) and (b) the (/(l)-hnear state Eq. (pl3|). 
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which has the same quantum order as in the C/(l)-linear 
state Eq. (|3|) (the staggered-flux state) . 

From Eq. (B132|) we get UlCnOxl ansatz 



FIG. 7: Contour plot of the spinon dispersion 

rmn{Ei{k) , E2{k)) as a function of {kx/2Ti,ky/2-K) for the 
C/(l) spin liquid state s, (a) is for the C/(l)-gapless state 
UlBOOOl in Eq. (|ll7| ) and (b) is for the (7(l)-linear state 
UlBOOll in Eq. (FW 
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From Eq. (B121) we get UlBOOOl ansatz 
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(115) 



(116) 



(117) 



and UlBOOll ansatz 



Ui^i+2x = — TjT 



Ui,i+y =?■(-)'" XT° 
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^0 
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i-yyGx ^Gy 
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■..H^i^ 


Gt ^i{-yr^ 



(118) 



Now let us discuss spinon dispersions in the above U{1 ) 
spin liquids. The spinon in the UlCnOxl state Eq. (115) 
has 4 linear nodes at (±7r/2, ±7r/2). Thus UlCnOxl state 
is a C/(l)-linear spin liquid. The UlCnOnl state Eq. ( |116| ) 
has fully gapped spinous and is a C/(l)-gapped spin liq- 
uid. 



The four spinon bands in the UlBOOOl state Eq. (|117| ) 
are given by (see Fig. 0a) 



±2x\l s\v? [kx) +sin^(fcj^)± {2rjCos{2kx 



2ricos{2ky) + i]i) 
(119) 



We find that the UlBOOOl state is a t/(l)-gapless spin 
liquifl The four spinon bands in the UlBOOll state 
Eq. (118) are given by (see Fig. 0b) 



±2x\/ sin^ (kx) + sm^{ky) ± 2?7(cos(2fca;) — cos(2A;j,)) 

(120) 
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Hence, the UlBOOll state is a C/(l)-linear spin liquid. 

To summarize we list all the spin liquids discussed so 
far in the following table: 



phase (E) is described by an ansatz 



Z2-gapped 


Z2Aa;a:0z 


Z2-linear 


Z2A0013, Z2Azzl3, Z2A00l7i 
Z2Azzln, Z2B0013, Z2Bzzl3 
Z2B001n, Z2Bzzln 


Z2-quadratic 


Z2Ba;2(12)n 


Z2-gapless 


Z2A2;2(12)n 


L/(l)-gapped 


UlCnOOx 


C/(l)-linear 


UlBOOll, UlCnOOx, UlCnOln 
UlCnOxl 


C/(l)-quadratic 


UlCxlOx 


t/(l)-gapless 


UlAOOOl, UlAOOll, UlAxlOx 
UlBOOOl, UlCnOls 


SU{2)-\mc&T 


SU2BnO 


5f7(2)-gapless 


SU2AnO 



VI. MEAN-FIELD PHASE DIAGRAM OF J1-J2 
MODEL 



To see which of the Z2, U{1), and SU{2) spin liquids 
discussed in the last section have low ground energies 
and may appear in real high Tc superconductors, we cal- 
culate the mean-field energy of a large class of transla- 
tion invariant ansatz. In Fig. ^ we present the result- 
ing mean-field phase diagram for a J1-J2 spin system. 
Here Ji is the nearest-neighbor spin coupling and J2 is 
the next-nearest-neighbor spin coupling. We have fixed 
Ji -\- J2 = 1. The 2/-axis is the mean-field energy per 
site (multiplied by a factor 8/3). The phase (A) is the 
TT-flux state (the SU2BnO S'C/(2)-hnear state) Eq. ^. 
The phase (B) is a state with two independent uniform 
RVB states on the diagonal links. It has SU{i) x SU{2) 
gauge fiuctuations at low energies and will be called an 
SU{2) X S'[/(2)-gapless state. Its ansatz is given by 



i,i+x+y X''' 
i.i-\-x — y X^ 

4=0 



(121) 



The phase (C) is a state with two independent tt-Aux 
states on the diagonal hnks. It has SU{2) x SU{2) gauge 
fiuctuations at low energies and will be called an SU{2) x 
S'J7(2)-linear state. Its ansatz is given by 



^z^x+x+y 



^i,i-\-x—y 
J 



-Xir'-r') 



ao-0 



(122) 



^i,i+x+y 


=XlT 


'U-i,i+x-y 


=XlT 


Ui,i+y 


=rjT^ 


Oo 


=0 



X2T 
X2T^ 



(123) 



The phase (D) is the chiral spin state Eq. (B3). The 



which break the 90° rotation symmetry and is a U{1)- 
linear state (see Fig. rab). Th e ph ase (F) is described by 
the UlCnOOx ansatz in Eq. ( p4|) . The UlCnOOx state 
can be a C7(l)-linear or a t/(l)-gapped state. The state 
for phase (F) turns out to be a t/(l)-gapped state. The 
phase (G) is described by the Z2Azzl3 ansatz in Eq. (8^) 
which is a Z2-linear state. The phase (H) is described 
by the Z2A0013 ansatz in Eq. ( [8^ ) and is also a Z2- 
linear state. The phase (I) is the uniform RVB state (the 
S'L/(2)-gapless state SU2AnO Eq. (|o|)). 

From Fig. Q, we see continuous phase transitions (at 
mean-field level) between the following pairs of phases: 
(A,D), (A,G), (B,G), (C,E), and (B,H). The three contin- 
uous transitions (B,G), (B,H) and (A,G) do not change 
any symmetries. We also note that the SU{2) gauge 
structure in the phase (A) breaks down to Z2 in the 
continuous transition from the phase (A) to the phase 
(G). The SU{2) x SU{2) gauge structure in the phase 
(B) breaks down to Z2 in the two transitions (B,G) and 
(B,H). 



VII. PHYSICAL MEASUREMENTS OF 
QUANTUM ORDERS 

After characterizing the quantum orders using PSG 
mathematically, we would like to ask how to measure 
quantum orders in experiments. The quantum orders 
in gapped states are related to the topological orders. 
The measurement of topological orders are discussed in 
Rcf. [^ ^, ^. The quantum order in a state with gap- 
less excitations can be measured, in general, by the dy- 
namical properties of gapless excitation. However, not 
all dynamical properties are universal. Thus we need to 
identify the universal properties of gapless excitations, 
before using them to characterize and measure quantum 
orders. The PSG characterization of quantum orders al- 
lows us obtain those universal properties. We simply 
need to identify the common properties of gapless exci- 
tations that are shared by all the ansatz with the same 
PSG. 

To demonstrate the above idea, we would like to study 
the spectrum of two-spinon excitations. We note that 
spinous can only be created in pairs. Thus the one- 
spinon spectrum is not physical. We also note that the 
two-spinon spectrum include spin-1 excitations which can 
be measured in experiments. At a given momentum, 
the two-spinon spectrum is distributed in one or several 
ranges of energy. Let E2s{k) be the lower edge of the 
two-spinon spectrum at momentum k. In the mean-field 
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FIG. 8: The mean-field energies for various phases in a Ji- 
J2 spin system. (A) the 7r-flux state (the Sf/(2)-hnear state 
SU2BnO). (B) the 517(2) x 5'C/(2)-gapless state in Eq. ( |l2l|) . 
(C) the SU{2) X 5'C/(2)-linear state in Eq. ( pj^ ) . (D) the chiral 
spin stat e (an St/(2)-gapped state). (E) the L'^(l)-linear state 
Eq. ( p^ which breaks 90° rotation symmetry. (F) the [/(!)- 
gapped state UlCnOOx in Eq. (104). (G) the Z2-linear state 
Z2A22I3 in Eq. (||). (H) the Z2-linear state Z2A0013 in 
Eq. (H). (I) the uniform RVB state (the S[/(2)-gapless state 
SU2AnO). 




FIG. 9: Contour plot of the dispersion for spin-1 excitation, 
E2s{k), as a function of (k^/2n. ky/ 2TT) for (a) the 5(7(2)- 
linear spin liquid SU2BnO in Eq. (^l|) (the 7r-flux ph ase) and 
(b) the (7(l)-quadratic spin liquid UlCilOa:: in Eq. (106). 



theory, the two-spinon spectrum can be constructed from 
the one-spinon dispersion 

£^2-spinon(fc) — -Bi_spinon (?) + £^l-spinon(fe — Q) (124) 

In Fig. |9|- ^we present mean- field E2S for some simple 
spin liquids. If the mean-field state is stable against the 
gauge fluctuations, we expect the mean-field i?2s should 
qualitatively agrees with the real E2S ■ 

Among our examples, there are eight Z2-linear spin 
liquids (see Fig. |l^ and Fig. ^. We see that some of 
those eight different Z2-linear spin liquids (or eight dif- 
ferent quantum orders) have different number of gapless 
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FIG. 10: Contour plot of E2s{k) as a function of 

{kx/2T7,ky/2n) for the Z2-linear spin liquids, (a) is for the 
Z2A0013 state in Eq. (||), (b) for the Z2Azzl3 state in 
Eq. (H), (c) for the Z2A001n state in Eq. (js^), (d) for the 
Z2A2:zln state in Eq. 
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FIG. 11: Contour plot of E2s{k) as a function of 

{kx/2n,ky/2TT) for the ^2 -linear spin liquids, (a) is for the 
Z2B0013 state in Eq. (||), (b) for the Z2Bz2l3 state in 
Eq. (|^), (c) for the Z2B001n state in Eq. (|^), (d) for the 
Z2B2:2ln state in Eq. 



points. The gapless points of some spin liquids are pinned 
at position k = (tt, it) and/or k — (tt, 0), (0, tt). By mea- 
suring the low energy spin excitations (say using neutron 
scattering), we can distinguish those Z2 spin liquids. We 
note that all the two-spinon spectra have rotation and 
parity symmetries around fe = 0. This is expected. Since 
the two-spinon spectra are physical, they should have all 
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FIG. 12: Contour plot ol E2s{k) as a function of 

{kx/2Tv,ky/2TT) for (a) the Z2-gapless state Z2Ax2(12)n in 
Eq. (g), and (b) the Z2-quadratic state Z2Ba;2(12)n in 
Eq. dnO). 
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FIG. 15: Contour plot of the two-spinon dispersion E2a(k) 
as a function of {kx/2-n, ky/2-K) for the (7(1) spin liqu id states, 
(a) is for the f/(l)-gapless state UlBOOOl in Eg. (pTTJ) and (b) 



is for the [/(l)-hnear state UlBOOll in Eq. (118) 
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FIG. 13: Contour plot of E2s(k) as a function of 

{kx/2-K,ky/2-n) for two L^(l)-linear spin liquids, (a) is for the 
UlCnOln state Eq. (P2|) (t he st aggered flux phase), and (b) 
for the UlCnOOa:: state Eq. (104) in the gapless phase. 



the symmetries the spin liquids have. 

We also have four f/(l)-linear spin liquids. Some of 
them can be distinguished by their different numbers of 
gapless points. It is interesting to note that all the U{1) 
spin liquids discussed here have a gapless point in the 
two-spinon spectrum pinned at position k = (7r,7r). The 
L''(l)-linear spin liquids are also different from the Z2- 
linear spin liquids in that the spin-spin correlations have 
different decay exponents once the C/(l) gauge fluctua- 
tions are included. We also see that i?2s has a quadratic 
form E2S oc k for the C/(l)-quadratic spin liquid. E^s 
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FIG. 14: Contour plot of the two-spinon dispersion E2s{k) 
as a function of {kx/2TT, ky/2n ) for (a) the [/(l)-linear spin 
liquid state Ul CwOx l in Eq. (115) and (b) the L''(l)-linear 
spin hquid Eq. (p^). 



vanishes in two finite regions in fe-space for the Z2-gapless 
spin liquids. 

Neutron scattering experiments probe the two-spinon 
sector. Thus low energy neutron scattering allows us to 
measure quantum orders in high Tc superconductors. 

Let us discuss the U{1) linear state UlCnOln (the 
staggered-flux state) in more detail. The UlCriOln 
state is proposed to describe the pseudo-gap metallic 
state in underdoped high T^ superconductors. [p3| p3 The 
UlCnOln state naturally explains the spin pseudo-gap 
in the underdoped metallic state. As an algebraic spin 
liquid, the UlCnOln state also explain the Luttinger-like 
electron spectral functionjpJ| and the enhancement of the 
(tt, tt) spin fluctuations [|67[ in the pseudo-gap state. From 
Fig. [13^, we see that gapless points of the spin-1 exci- 
tations in the UlCnOln state are always at fc = {tt,tt), 
(0, 0), (tt, 0) and (0, tt). The equal energy contour for the 
edge of the spin-1 continuum has a shape of two over- 
lapped ellipses at all the four k points. Also the energy 
contours are not perpendicular to the zone boundary. All 
those are the universal properties of the UlCnOln state. 
Measuring those properties in neutron scattering experi- 
ments will allow us to determine if the pseudo-gap metal- 
lic state is described by the UlCnOln (the staggered-flux) 
state or not. 

We have seen that at low energies, the UlCnOln 
state is unstable due to the instanton effect. Thus the 
UlCnOln state has to change into some other states, such 
as the 8 Z2 spin liquids discussed in section M or some 
other states not discussed in this paper. From Fig. 10 1, 
we see that the transition from the UlCnOln state to the 
Z2-linear state Z2A0013 can be detected by neutron scat- 
tering if one observe the splitting of the node at (tt, tt) 
into four nodes at (tt ± J, tt ± S) and the splitting of the 
nodes at (tt, O) and (0, tt) into two nodes at (tt ± 6, 0) 
and (0,7r ± S). From Fig. OQb, we see that, for the 
transition from the UlCnOln state to the Z2-linear state 
Z2Azzl3, the node at {tt,7t) still splits into four nodes 
at (tt ± (5, TT ± 5) . However, the nodes at (tt, 0) and (0, tt) 
split differently into two nodes at (tt, ±5) and {±5, n). We 
can also study the transition from the UlCnOln state to 
other 6 Z2 spin liquids. We find the spectrum of spin- 
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1 excitations all change in certain characteristic ways. 
Thus by measuring the spin-1 excitation spectrum and 
its evolution, we not only can detect a quantum transi- 
tion that do not change any symmetries, we can also tell 
which transition is happening. 

The neutron scattering on high T^ superconductor in- 
deed showed a splitting of the scattering peak at (tt, tt) 
into four peaks at (tt ± (5, tt) , (tt, tt ± (5) 1 30 , 68, p3, uOl u\ 
\!% |7|, [74| o r into two peaks at (tt, tt) ^> (tt+iJ, 'k—5)^{'k— 
5, 7r-|-(5) |28|, [75| as we lower the energy. This is consistent 
with our belief that the UlCnOln state is unstable at low 
energies. However, it is still unclear if we can identify the 
position of the neutron scattering peak as the position of 
the node in the spin-1 spectrum. If we do identify the 
scattering peak as the node, then non of the 8 Z2 spin 
liquids in the neighborhood of the UlCnOln state can 
explain the splitting pattern (tt ± 5, tt), (tt, tt ± 5). This 
will imply that the UlCnOln state change into another 
state not studied in this paper. This example illustrates 
that detailed neutron scattering experiments are power- 
ful tools in detecting quantum orders and studying new 
transitions between quantum orders that may not change 
any symmetries. 



VIII. FOUR CLASSES OF SPIN LIQUIDS AND 
THEIR STABILITY 



to Chern-Simons terms or Anderson-Higgs mechanism. 
By definition, the rigid spin liquids are locally stable and 
self consistent. The rigid spin liquids are characterized 
by topological orders and they have the true spin-charge 
separation. The low energy effective theories for rigid 
spin liquids are topological field theories. The Z2-gapped 
spin liquid and chiral spin liquid are examples of rigid 
spin liquids. 



B. Bose spin liquid 

The C/(l)-gapped spin liquid discussed in the last sec- 
tion is not a rigid spin liquid. It is a Bose spin liquid. 
Although the spinon excitations are gapped, the U{\) 
gauge fluctuations are gapless in the C/(l)-gapped spin 
liquid. The dynamics of the gapless U[\) gauge fluctua- 
tions are described by low energy effective theory 
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(125) 



where f^^ is the field strength of the U{1) gauge field. 
However, in l-|-2 dimension and after including the in- 
stanton effect, the U{\) gauge fluctuations will gain an 
energy gap.|Q The properties of the resulting quantum 
state remain to be an open problem. 



We have concentrated on the mean-field states of 
spin liquids and presented many examples of mean-field 
ansatz for symmetric spin liquids. In order for those 
mean-field states to represent real physical spin liquids, 
we need to include the gauge fluctuations. We also 
need to show that the inclusion of the gauge fluctuations 
does not destabilize the mean-field states at low ener- 
gies. This requires that (a) the gauge interaction is not 
too strong and (b) the gauge interaction is not a rele- 
vant perturbation. (The gauge interaction, however, can 
be a marginal perturbation.) The requirement (a) can 
be satisfied through large N limit and/or adjustment of 
short-range spin couplings in the spin Hamiltonian, if 
necessary. Here we will mainly consider the requirement 
(b). We find that, at least in certain large N limits, 
many (but not all) mean-field states do correspond to 
real quantum spin liquids which are stable at low ener- 
gies. In this case, the characterization of the mean-field 
states by PSG's correspond to the characterization of real 
quantum spin liquids. 

All spin liquids (with odd number of electron per unit 
cell) studied so far can be divided into four classes. In 
the following we will study each classes in turn. 



A. Rigid spin liquid 

In rigid spin liquids, by definition, the spinous and all 
other excitations are fully gaped. The gapped gauge field 
only induces short range interaction between spinous due 



C. Fermi spin liquid 

The Fermi spin liquids have gapless excitations that are 
described by spin 1/2 fermions. Those gapless excitations 
have only short range interactions between them. The 
Z2-linear, Z2-quadratic and the .Z2-gapless spin liquid 
discussed above are examples of the Fermi spin liquids. 

The spinous have a massless Dirac dispersion in Z2- 
linear spin liquids. Thus Z2-linear spin liquids are lo- 
cally stable since short range interactions between mass- 
less Dirac fermions are irrelevant at 1-1-2 dimensions. We 
would like to point out that the massless Dirac dispersion 
of the Z2-linear spin liquids are protected by the PSG (or 
the quantum order). That is any perturbations around, 
for example, the Z2-linear ansatz Eq. (p9) cannot destroy 
the massless Dirac dispersion as long as the PSG are not 
changed by the perturbations. To understand this result, 
we start with the most general form of symmetric pertur- 
bations Eq. ( |B^ ) around the Z2-linear ansatz Eq. (|3S|). 
We find that such perturbations vanish in the momen- 
tum space at fc = (±7r/2, ±7r/2). The translation, parity, 
and the time reversal symmetries do not allow any mass 
terms or chemical potential terms. Thus the Z2-linear 
spin liquid is a phase that occupy a finite region in the 
phase space (at T — Q). One does not need any fine tun- 
ing of coupling constants and Uij to get massless Dirac 
spectrum. 

Now let us con sider the stability of the Z2-quadratic 
spin liquid Eq. (110). The spinous have a gapless 
quadratic dispersion in the 2'2-quadratic spin liquid. The 
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gapless quadratic dispersion of the Z^2-quadratic spin 
liquid is also protected by the symmetries. The most 
general form of symmetr ic p erturbations arou nd the 
Z2-quadratic ansatz Eq. (IIC) is given by Eq. (B102) 
(Z2Ba;2(12)n). In the momentum space, the most gen- 
eral symmetric Z2-quadratic ansatz give rise to the fol- 
lowing Hamiltonian (after considering the 90° rotation 
symmetry) 

H = - 2 ^ Xm«[sin(nfc^ - mky)To + sin(mfc^ + nky)Ti] 

4- 2 \ 7y,„„[cos(7ifcj. — mky)T2 + cos[mkj. + nky)T^] 

+ 2 2^ ^mn[cos{nkx — mky)Ti -\- cos{mkx + nky)T5] 

(126) 

where 
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and the summation is over m =even, n =odd. We find 
that the spinon dispersion still vanish at fc = (0, 0), (0, tt) 
and the energy still satisfy E oc k . The translation, par- 
ity, and the time reversal symmetric perturbations do not 
change the qualitative behavior of the low energy spinon 
dispersion. Thus, at mean-field level, the Z2-quadratic 
spin liquid is a phase that occupy a finite region in the 
phase space (at T = 0). One does not need any fine tun- 
ing of coupling constants to get gapless quadratic disper- 
sion of the spinous. However, unlike the Z2-linear spin 
liquid, the short range four-fermion interactions between 
the gapless spinous in the .Z2-quadratic state are marginal 
at 1-1-2 dimensions. Further studies are needed to under- 
stand the dynamical stability of the Z2-quadratic spin 
liquid beyond the mean-field level. 

The Z2-gapless spin liquid is as stable as Fermi liquid in 
1-1-2 dimensions. Again we expect Z2-gapless spin liquid 
to be a phase that occupy a finite region in the phase 
space, at least at mean-field level. 



D. Algebraic spin liquid 

C/(l)-lincar spin liquids are examples of algebraic spin 
liquids. Their low lying excitations are described by 
massless Dirac fermions coupled to U{1) gauge field. 
Although the massless Dirac fermions are protected by 
quantum orders, the gauge couplings remain large at low 
energies. Thus the low lying excitations in the U{1)- 
linear spin liquids are not described by free fermions. 
This makes the discussion on the stability of those states 
much more difficult. 

Here we would like to concentrate on the t/(l)-linear 
spin liquid UlCnOln in Eq. (|3^). The spinous have a 
massless Dirac dispersion in the C/(l)-linear spin liquid. 
First we would like to know if the massless Dirac disper- 
sion is generic property of the C/(l)-linear spin liquid, ie if 



the massless Dirac dispersion is a property shared by all 
the spin liquids that have the same quantum order as that 
in Eq. (p2). The most general perturbations arou nd the 
C/(l)-linear ansatz Eq. ( p2|) are given by Eq. (B6), if the 
perturbations respect translation, parity, and the time re- 
versal symmetries, and if the perturbation do not break 
the C/(l) gauge structure. Since Suf^ = (5u^ — for 
m — even, their contributions in the momentum space 
vanish at fe = (0,0) and k = (0,7r). The spinon en- 
ergy also vanish at those points for the ansatz Eq. ([3^). 
Thus the massless Dirac dispersion is protected by the 
symmetries and the U{1) gauge structure in the U{1)- 
linear spin liquid Eq. (p2|). In other words, the massless 
Dirac dispersion is protected by the quantum order in 
the C/(l)-linear spin liquid. 

Next we consider if the symmetries and the U{1) gauge 
structure in the C/(l)-linear spin liquid can be broken 
spontaneously due to interactions/fluctuations at low en- 
ergy. The low energy effective theory is described by La- 
grangian (in imaginary time) 



C = 



yal°K,and^+^a^Ma (128) 
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where /i = 0,1,2, a — 1,2, 7^ are 4x4 7-matrices, 
fg = 1, and {vi,a,V2.a) are velocities for a^ fermion in 
X and y directions. We make a large N generalization 
of the above effective theory and allow a = 1,2,..,N. 
Our first concern is about whether the self energy from 
the gauge interaction can generate any mass/chemical- 
potential term, due to infrared divergence. It turns out 
that, in the 1/N expansion, the gauge fluctuations repre- 
sent an exact marginal perturbation that does not gener- 
ate any mass/chemical=potential term.|Q Instead the 
gauge interaction changes the quantum fixed point de- 
scribed by free massless Dirac fermions to a new quantum 
fixed point which has no free fermionic excitations at low 
energies. 1^, |7^ The new quantum fixed point has gapless 
excitations and correlation functions all have algebraic 
decay. Such a quantum fixed point was called algebraic 
spin liquid. [Q Actually, it is easy to understand why the 
gauge fluctuations represent an exact marginal perturba- 
tion. This is because the conserved current that couple to 
the gauge potential cannot have any anomalous dimen- 
sions. Thus if the gauge interaction is marginal at first 
order, then it is marginal at all orders. Gauge interaction 
as an exact marginal perturbation is also supported by 
the following results. The gauge invariant Green's func- 
tion of ^ is found to be gapless after coupling to gauge 
field, to all orders in 1/iV expansion.]?^ Recently it was 
argued that the C/(l) gauge interaction do not generate 
any mass perturbatively even when A'^ is as small as 2.|7q| 

Now let us discuss other possible instabilities. First we 
would like consider a possible instability that change the 
[/(l)-linear state to the Z2-linear state. To study such an 
instability we add a charge-2 Higgs field to our effective 
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theory 

C^^j]j°[Y{d^ + ia^)]xj;a (129) 

+ 1(^0 - 2mo)0l' + v\d, - 2m,)0p + V{(j)) 

where V{(t)) has its minimum at (j) = and we have as- 
sumed vi^a = V2,a = 1 for simphcity. (Note that (j) corre- 
sponds to A in Eq. ( ^ . It is a non-zero A that break the 
C/(l) gauge structure down to the Z2 gauge structure.) If 
after integrating out ^ and a^, the resulting effective po- 
tential V'e//(|0|) has its minimum at a non zero (/>, then 
the C/(l)-linear state has a instability towards the Z2- 
linear state. 

To calculate T4//(|(?!)|), we first integrate out ip and get 
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Now the effective potential Ve//(|0|) can be obtained by 
integrating out a^ (in the ao = gauge) and the phase 9 
of the field, = pe*^:i§ 
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We find that K// = V^ — Ci|0|^ In |(/)| where C'l is a con- 
stant. Now it is clear that the gapless gauge fluctuations 
cannot shift the minimum of V from = and the U{1)- 
linear state is stable against spontaneously changing into 
the .Z2-linear state. 

So far we only considered the effects of perturbative 
fluctuations. The non-perturbative instanton effects can 
also cause instability of the algebraic spin liquid. The 
instanton effects have been discussed in Ref. [Q for the 
case w^ = w^. It was found that the instanton effects 
represent a relevant perturbation which can destabilize 
the algebraic spin liquid when N < 24. In the following, 
we will generalize the analysis of Ref. [M to v^ ^ v^ 
case. First we rewrite 

d^k 1 

S = I j;^i^a^,{-k)TT^,ya^{k) 

d^k 1 

--f^(-k)K^,f,(k) (134) 
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When 7r„ 



/p = e^^xduO-x (135) 

fc2(5^^ - kfj^kiy, we find Kf,„ = 5^^- When 
TT/jiy = {k^Sfji, — kfj_kiy)/vk^, we may assume Kf^^, = 
Sf_iu/vk^- When vi^a 7^ V2,a we have 

(/v) = y: k^' + viA + vij^i)-'" X 
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(136) 



The instanton field /^ minimize the action Eq. (134) and 
satisfies 

i^p./. = c(fc)fc^ (137) 

where c(fc) is chosen such that fc^/^^ — 2iTr We find that 

8ccij 



/o 
/i = 

/2 = 



8cki 

Ea(^'+<afc?+4a^l)-^/VaM,a 

8ck2 

T.a{^'+vlakl+vlakl)-'"^l,alv2.,a 



(138) 



and 
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Using the above solution, we can calculate the action for 
a single instanton, which has a form 



TV 



a{v2/vi)\n{L) 



(140) 



where L is the size of the system and we have assumed 
that N/2 fermions have velocity {vx,Vy) = (wi,W2) and 
the other N/2 fermions have velocity {vxTVy) = (u2,i'i). 
We find a(l) = 1/4 + 0{l/N) and a(0.003) = 3.-1- 
0{1/N). When ^a{v2/vi) > 3, the instanton effect is 
irrelevant. We see that even for the case N = 2, the in- 
stanton effect can be irrelevant for small enough W2/i'i. 
Therefore, the algebraic spin liquid exists and can be sta- 
ble. 

It has been proposed that the pseudo-gap metallic 
state in underdoped high T^ superconductors is described 
by the (doped) staggered flux state (the C/(l)-linear 
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state UlCriOln) which contains a long range U{1) gauge 
interaction. |33[ |34| From the above result, we see that, 
for reahstic f2/i'i '--^ 0.1 in high Tc superconductors, the 
UlCnOl?! spin hquid is unstable at low energies. How- 
ever, this does not mean that we cannot not use the al- 
gebraic spin liquid UlCnOln to describe the pseudo-gap 
metallic state. It simply means that, at low temper- 
atures, the algebraic spin liquid will change into other 
stable quantum states, such as superconducting state or 
antiferromagnetic state |]79[ as observed in experiments. 

The unstable algebraic spin liquid can be viewed as 
an unstable quantum fixed point. Thus the algebraic- 
spin-liquid approach to the pseudo-gap metallic state 
in underdoped samples looks similar to the quantum- 
critical-point approach po|, pl[ . However, there is an im- 
portant distinction between the two approaches. The 
quantum-critical-point approach assumes a nearby con- 
tinuous phase transition that changes symmetries and 
strong fluctuations of local order parameters that cause 
the criticality. The algebraic-spin-liquid approach does 
not require any nearby symmetry breaking state and 
there is no local order parameter to fluctuate. 



E. Quantum order and the stability of spin liquids 

After introducing quantum orders and PSG, we can 
have a deeper discussion on the stability of mean-field 
states. The existence of the algebraic spin liquid is a 
very striking phenomenon, since gapless excitations in- 
teract down to zero energy and cannot be described by 
free fermions or free bosons. According to a conventional 
wisdom, if bosons/fermions interact at low energies, the 
interaction will open an energy gap for those low lying ex- 
citations. This implies that a system can either has free 
bosonic/fermionic excitations at low energies or has no 
low energy excitations at all. According to the discussion 
in section VIII, such a conventional wisdom is incorrect. 
But it nevertheless rises an important question: what 



protects gapless excitations (in particular when they in- 
teract at all energy scales). There should be a "reason" 
or "principle" for the existence of the gapless excitations. 
Here we would like to propose that it is the quantum or- 
der that protects the gapless excitations. We would like 
to stress that gapless excitations in the Fermi spin liq- 
uids and in the algebraic spin liquids exist even without 
any spontaneous symmetry breaking and they are not 
protected by symmetries. The existence of gapless ex- 
citations without symmetry breaking is a truly remark- 
able feature of quantum ordered states. In addition to 
the gapless Nambu-Goldstone modes from spontaneous 
continuous symmetry breaking, quantum orders offer an- 
other origin for gapless excitations. 

We have seen from several examples discussed in sec- 
tion that the quantum order (or the PSG) not only 
protect the zero energy gap, it also protects certain qual- 
itative properties of the low energy excitations. Those 
properties include the linear, quadratic, or gapless disper- 



sions, the k locations where the 2-spinon energy i?2s(fe) 
vanishes, etc . 

Since quantum order is a generic property for any 
quantum state at zero temperature, we expect that 
the existence of interacting gapless excitations is also a 
generic property of quantum state. We see that algebraic 
state is a norm. It is the Fermi liquid state that is special. 

In the following, we would like to argue that the PSG 
can be a stable (or universal) property of a quantum 
state. It is robust against perturbative fluctuations. 
Thus, the PSG, as a universal property, can be used to 
characterize a quantum phase. From the examples dis- 
cussed in sections VIII C and VIII D, we see that PSG 



protects gapless excitations. Thus, the stability of PSG 
also imply the stability of gapless excitations. 

We know that a mean-field spin liquid state is char- 
acterized by Uij — (ipiip-). If we include perturbative 
fluctuations around the mean-field state, we expect Uij 
to receive perturbative corrections SUij. Here we would 
like to argue that the perturbative fluctuations can only 
change Uij in such a way that Uij and Uij + 6Uij have 
the same PSG. 

First we would like to note the following well know 
facts: the perturbative fluctuations cannot change the 
symmetries and the gauge structures. For example, if 
Uij and the Hamiltonian have a symmetry, then SUij 
generated by perturbative fluctuations will have the same 
symmetry. Similarly, the perturbative fluctuations can- 
not generate SUij that, for example, break a U{1) gauge 
structure down to a Z2 gauge structure. 

Since both the gauge structure (described by the IGG) 
and the symmetry are part of the PSG, it is reasonable 
to generalize the above observation by saying that not 
only the IGG and the symmetry in the PSG cannot be 
changed, the whole PSG cannot be changed by the per- 
turbative fluctuations. In fact, the mean-field Hamilto- 
nian and the mean-field ground state are invariant under 
the transformations in the PSG. Thus in a perturbative 
calculation around a mean-field state, the transforma- 
tions in the PSG behave just like symmetry transforma- 
tions. Therefore, the perturbative fluctuations can only 
generate SUij that are invariant under the transforma- 
tions in the PSG. 

Since the perturbative fluctuations (by definition) do 
not change the phase, Uij and Uij + SUij describe the 
same phase. In other words, we can group Uij into classes 
(which are called universality classes) such that Uij in 
each class are connected by the the perturbative fluctu- 
ations and describe the same phase. We see that if the 
above argument is true then the universality classes are 
classified by the PSG's (or quantum orders). 

We would like to point out that we have assumed the 
perturbative fluctuations to have no infrared divergence 
in the above discussion. The infrared divergence implies 
the perturbative fluctuations to be relevant perturba- 
tions, which cause phase transitions. 
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IX. RELATION TO PREVIOUSLY 
CONSTRUCTED SPIN LIQUIDS 

Since the discovery of high T^ superconductor in 1987, 
many spin hquids were constructed. After classifying and 
constructing a large class of spin liquids, we would like 
to understand the relation between the previously con- 
structed spin liquids and spin liquids constructed in this 
paper. 

Anderson, Baskaran, and ZouJMl 11^, 11^ first used the 
slave boson approach to construct uniform RVB state. 
The uniform RVB state is a symmetric spin liquid which 
has all the symmetries of the lattice. It is a S'C/(2)-gapless 
state characterized by the PSG SU2AnO. Later two more 
spin liquids were constructed using the same [/(I) slave 
boson approach. One is the 7r-flux phase and the other 
is the staggered-flux/d-wave state. [^l| ^2], ^ The tt-Aux 
phase is a S'C/(2)-linear symmetric spin liquid character- 
ized by PSG SU2BnO. The staggered-flux/d-wave state 
is a J7(l)-linear symmetric spin liquid characterized by 
PSG UlCnOln. The UlCnOln state is found to be the 
mean-field ground for underdoped samples. Upon doping 
the UlCnOln state becomes a metal with a pseudo-gap 
at high temperatures and a d-wave superconductor at low 
temperatures. 

It is amazing to see that the slave boson approach, 
which is regarded as a very unreliable approach, pre- 
dicted the d-wave superconducting state 5 years before 
its experimental confirmation. ||2l|, ^ |8J, ^ Maybe pre- 
dicting the d-wave superconductor is not a big deal. Af- 
ter all, the d-wave superconductor is a commonly known 
state and the paramagnon approach |p6[ p7[ predicted d- 
wave superconductor before the slave boson approach. 
However, what is really a big deal is that the slave bo- 
son approach also predicted the pseudo-gap metal which 
is a completely new state. It is very rare in condensed 
matter physics to predict a new state of matter before 
experiments. It is also interesting to see that not many 
people believe in the slave boson approach despite such 
a success. 

The above C/(l) and SU{2) spin liquids are likely to 
be unstable at low energies and may not appear as the 
ground states of spin systems. The first known stable 
spin liquid is the chiral spin liquid, g, g]. It has a true 
spin-charge separation. The spinous and holons carry 
fractional statistics. Such a state breaks the time rever- 
sal and parity symmetries and is a S'f/(2)-gapped state. 
The SU{2) gauge fluctuations in the chiral spin state does 
not cause any instability since the gauge fluctuations are 
suppressed and become massive due to the Chern-Simons 
term. Due to the broken time reversal and parity sym- 
metries, the chiral spin state does not flt within our clas- 
sification scheme. 

Spin liquids can also be constructed using the slave- 
fermion/tr-model approach. J4|, |42) Some gapped spin liq- 
uids were constructed using this approach. Gsl M Those 
states turn out to be Z2 spin liquids. But they are not 
symmetric spin liquids since the 90° rotation symmetry 



is broken. Thus they do not fit within our classifica- 
tion scheme. Later, a .Z2-gapped symmetric spin liquid 
was constructed using the 5/7(2) slave-boson approach 
(or the SU{2) projective construction) . Hll The PSG for 
such a state is Z2Aa;x0z. Recently, another Z2 state was 
constructed using slave-boson approach. ||63[ Q It is a 
Z2-linear symmetric spin liquid. Its PSG is given by 
Z2A0013. New Z2 spin liquids were also obtained re- 
cently using the slave- fermion/ a- model approach. Q It 
appears that most of those states break certain symme- 
tries and are not symmetric spin liquids. We would like 
to mention that Z2 spin liquids have a nice property that 
they are stable at low energies and can appear as the 
ground states of spin systems. 

Many spin liquids were also obtained in quantum 
dimer model, |t47|, [4§|j 4£,j50| ^ and in various numer- 
ical approach. |32[ ]5^, 54, BM It is hard to compare those 
states with the spin liquids constructed here. This is be- 
cause either the spectrum of spin-1 excitations was not 
calculated or the model has a very different symmetry 
than the model discussed here. We need to generalized 
our classification to models with different symmetries so 
that we can have a direct comparison with those interest- 
ing results and with the non-symmetric spin liquids ob- 
tained in the slave-fermion/cr-model approach. In quan- 
tum dimer model and in numerical approach, we usually 
know the explicit form of ground state wave function. 
However, at moment, we do not know how to obtain PSG 
from ground wave function. Thus, knowing the explicit 
ground state wave function does not help us to obtain 
PSG. We see that it is important to understand the re- 
lation between the ground state wave function and PSG 
so that we can understand quantum order in the states 
obtained in numerical calculations. 



X. SUMMARY OF THE MAIN RESULTS 

In the following we will list the main results obtained 
in this paper. The summary also serves as a guide of the 
whole paper. 

(1) A concept of quantum order is introduced. The 
quantum order describes the orders in zero-temperature 
quantum states. The opposite of quantum order - classi- 
cal order describes the orders in finite-temperature clas- 
sical states. Mathematically, the quantum order charac- 
terizes universality classes of complex ground state wave- 
functions. It is richer then the classical order that char- 
acterizes the universality classes of positive distribution 
functions. Quantum orders cannot be completely de- 
scribed by symmetries and order parameters. Landau's 
theory of orders and phase transitions does not apply to 
quantum orders. (See section I A) 

(2) Projective symmetry group is introduced to de- 
scribe different quantum orders. It is argued that PSG 
is a universal property of a quantum phase. PSG ex- 
tends the symmetry group description of classical orders 
and can distinguishes different quantum orders with the 



31 



same symmetries. (See section EVA and VIII E) 

(3) As an application of the PSG description of quan- 
tum phases, we propose the following principle that gov- 
ern the continuous phase transition between quantum 
phases. Let PSGi and PSG2 be the PSG's of the two 
quantum phases on the two sides of a transition, and 
PSGcr be the PSG that describes the quantum criti- 
cal state. Then PSGi C PSGcr and PSG2 C PSGcr- 
We note that the two quantum phases may have the 
same symmetry and continuous quantum phase transi- 
tions are possible between quantum phases with same 
symmetry. pM The continuous transitions between dif- 
ferent mean-field symmetric spin liquids are discussed 
in section M and appendix B which demonstrate the 
above principle. However, for continuous transitions be- 
tween mean-field states, we have an additional condition 

PSGi = PSGcr or PSG2 = PSGcr- 

(4) With the help of PSG, we find that, within the 
SU{2) mean-field slave-boson approach, there are 4 sym- 
metric SU{2) spin liquids and infinite many symmetric 
U{1) spin liquids. There are at least 103 and at most 196 
symmetric Z2 spin liquids. Those symmetric spin liquids 
have translation, rotation, parity and the time reversal 
symmetries. Although the classifications are done for the 
mean-field states, they apply to real physical spin liquids 
if the corresponding mean-field states turn out to be sta- 
ble against fiuctuations. (See section tV and appendix 
A) 

(5) The stability of mean-field spin liquid states are dis- 
cussed in detail. We find many gapless mean-field spin 
liquids to be stable against quantum fluctuations. They 
can be stable even in the presence of long range gauge in- 
teractions. In that case the mean-field spin liquid states 
become algebraic spin liquids where the gaples s exc ita- 
tions interact down to zero energy. (See section VIII) 

(6) The existence of algebraic spin liquids is a striking 
phenomenon since there is no spontaneous broken sym- 
metry to protect the gapless excitations. There should be 
a "principle" that prevents the interacting gapless exci- 
tations from opening an energy gap and makes the alge- 
braic spin liquids stable. We propose that quantum order 
is such a principle. To support our idea, we showed that 
just like the symmetry group of a classical state deter- 
mines the gapless Nambu-Goldstone modes, the PSG of 
a quantum state determines the structure of gapless exci- 
tations. The gauge group of the low energy gauge fiuctu- 
ations is given by the IGG, a subgroup of the PSG. The 
PSG also protects massless Dirac fermions from gaining 
a mass due to radiative corrections. We see that the sta- 
bilities of algebraic spin liquids and Fermi spin liquids 
are protected by their PSG's. The existence of gapless 
excitations (the gauge bosons and gapless fermions) with- 
out symmetry breaking is a truly remarkable feature of 
quantum ordered states. The gapless gauge and fermion 
excitations are originated from the quantum orders, just 
like the phonons are originat ed from transla t ion sym me- 
try breaking. (See sections |VIII C| , |VIIID| , |VIIIE| and 
discussions below Eq. (E^) 



(7) Many Z2 spin liquids are constructed. Their low 
energy excitations are described by free fermions. Some 
Z2 spin liquids have gapless excitations and others have 
finite energy gap. For those gapless Z2 spin liquids some 
have Fermi surface while others have only Fermi points. 
The spinon dispersion near the Fermi points can be lin- 
ear E oc \k\ (which gives us Z2-linear spin liquids) or 
quadratic E oc k (which gives us Z2-quadratic spin liq- 
uids). In particular, we find there can be many Z2-linear 
spin liquids with different quantum orders. All those dif- 
ferent Z2-linear spin liquids have nodal spinon excita- 
tions. (See section |l|, 0, and appendix B) 

(8) Many C/(l) spin liquids are constructed. Some C/(l) 
spin liquids have gapless excitations near isolated Fermi 
point with a linear dispersion. Those U{1) linear states 
can be stable against quantum fluctuations. Due to long 
range U{1) gauge fluctuations, the gapless excitations in- 
teract at low energies. The [/(l)-linear spin liquids can 
be concrete realizations of algebraic spin liquids, pj, p% 
(See section |l|, and appendix B) 

(9) Spin liquids with the same symmetry and different 
quantum orders can have continuous phase transitions 
between them. Those phase transitions are very simi- 
lar to the continuous topological phase transitions be- 



tween quantum Hall states.] 66 



We find that, 
at mean-field level, the UlCnOln spin liquid in Eq. (132 



(the staggered flux phase) can continuously change into 
8 different symmetric Z2 spin liquids. The SU2AnO spin 
liquid in Eq. ( po| ) (the uniform RVB state) can contin- 
uously change into 12 symmetric U{1) spin liquids and 
52 symmetric Z2 spin liquids. The SU2BnO spin liquid 
in Eq. (|3l|) (the 7r-flux phase) can continuously change 
into 12 symmetric U{1) spin liquids and 58 symmetric 
Z2 spin liquids. (See appendix B) 

(10) We show that spectrum of spin-1 excitations 
(ie the two-spinon spectrum), which can be probed in 
neutron scattering experiments, can be used to measure 
quantum orders. The gapless points of the spin-1 exci- 
tations in the UlCnOln (the staggered-flux) state are al- 
ways at fc = (tt, it), (0, 0), (tt, 0) and (0, n). In the pseudo- 
gap metallic phase of underdopcd high Tc superconduc- 
tors, the observed splitting of the neutron scattering peak 
(^, tt) ^ (7r±<5, tt), (^, 7r±<5) @, ^, ^ 0,jg, fl g |7| 
or (tt, tt) -^ {tt + 6,TT ~ S),{7r — S,7r + 6) uS, ITm at low 
energies indicates a transition of the UlCnOln state into 
a state with a different quantum order, if we can indeed 
identify the scattering peak as the gapless node. Non of 
the 8 symmetric Z2 spin liquids in the neighborhood of 
the the UlCnOln state can explain the splitting pattern. 
Thus we might need to construct a new low energy state 
to explain the splitting. This illustrates that detailed 
neutron scattering experiments are powerful tools in de- 
tecting quantum orders and st udyi ng transitions between 



quantum orders. (See section VII 



(11) The mean-field phase diagram Fig. ^ for a Ji- 
J2 spin system is calculated. (Only translation sym- 
metric states are considered.) We find four mean-field 
ground states as we change J2/J1'- the 7r-flux state (the 



32 



SU2AnO state), the chiral spin st ate (an S'C/(2)-gapped 
state), the [/(l)-hnear state Eq. (123) which breaks 90° 
rotat ion symmetry, and the SU{2) x SU{2)-\mear state 
Eq. (122). We also find several locally stable mean -field 
states: the [/(l)-gapped state UlCnOOa; in Eq. (104) and 
two Z2-linear states Z2Azzl3 in Eq. (|6|) and Z2A0013 
in Eq. (pq). Those spin liquids have a better chance 
to appear in underdoped high Tc superconductors. The 
Z2A0013 Z2-linear state has a spinon dispersion very sim- 
ilar to electron dispersion observed in underdoped sam- 
ples. The spinon dispersion in the Z2Azzl3 Z2-linear 
state may also be consistent with electron dispersion in 
underdoped samples. We note that the two-spinon spec- 
trum for the two Z2-linear states have some qualitative 
differences (see Fig. \\v^ and Fig. [lO| b and note the posi- 
tions of the nodes) . Thus we can use neutron scattering 
to distinguish the two states. (See section VI.) 

Next we list some remarks/comments that may clar- 
ify certain confusing points and help to avoid possible 
misunderstanding. 

(A) Gauge structure is simply a redundant labeling 
of quantum states. The "gauge symmetries" (referring 
different labels of same physical state give rise to the 
same result) are not symmetries and can never be broken. 
(See the discussion below Eq. (fsl)) 

(B) The gauge structures referred in this paper (such 
as in Z2, U{1), or SU{2) spin fiquids) are "low energy" 
gauge structures. They are different from the "high 
energy" gauge structure that appear in Z2, U{\), and 
SU{2) slave-boson approaches. The "low energy" gauge 
structures are properties of the quantum orders in the 
ground state of a spin system. The "high energy" gauge 
structure is a particular way of writing down the Hamil- 
tonian of spin systems. The two kinds of gauge structures 
have nothing to do with each other. (See discussions at 
the end of section I C and at the end of section IV A ) 

(C) There are (at least) two different interpretations of 
spin-charge separation. The first interpretation (pseudo 
spin-charge separation) simply means that the low energy 
excitations cannot be described by electron-like quasipar- 
ticles. The second interpretation (true spin-charge sepa- 
ration) means the existence of free spin-1/2 neutral quasi- 
particles and spin-0 charged quasiparticles. In this paper 
both interpretations are used. The algebraic spin liquids 
have a pseudo spin-charge separation. The Z2 and chi- 
ral spin liquids have a true spin-charge separation. (See 
section I C) 

(D) Although in this paper we stress that quantum 
orders can be characterized by the PSG's, we need to 
point out that the PSG's do not completely characterize 
quantum orders. Two different quantum orders may be 
characterized by the sam e PS G. As an example, we have 
seen that the ansatz Eq. (104) can be a C/(l)-linear state 
or a f/(l)-gapped state depending on the values of pa- 
rameters in the ansatz. Both states are described by the 
same PSG UlCnOOx. Thus the PSG can not distinguish 
the different quantum orders carried by the t/(l)-linear 
state and the C/(l)-gapped state. 



(E) The unstable spin liquids can be important in 
understanding the finite temperature states in high Tc 
superconductors. The pseudo-gap metallic state in un- 
derdoped samples is likely to be described by the un- 
stable UlCnOln algebraic spin liquid (the staggered 
flux state) which contains a long range U{1) gauge 
interaction. I33I M (See discussions at the end of section 
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(F) Although we have been concentrated on the char- 
acterization of stable quantum states, quantum order and 
the PSG characterization can also be used to describe 
the internal order of quantum critical states. Here we 
define "quantum critical states" as states that appear 
at the continuous phase transition points between two 
states with different symmetries or between two states 
with different quantum orders (but the same symme- 
try). We would like to point out that "quantum criti- 
cal states" thus defined are more general than "quantum 
critical points" . "Quantum critical points" , by defini- 
tion, are the continuous phase transition points between 
two states with different symmetries. The distinction 
is important. "Quantum critical points" are associated 
with broken symmetries and order parameters. Thus the 
low energy excitations at "quantum critical points" come 
from the strong fluctuations of order parameters. While 
"quantum critical states" may not be related to broken 
symmetries and order parameters. In that case it is im- 
possible to relate the gapless fluctuations in a "quantum 
critical state" to fluctuations of an order parameter. The 
unstable spin liquids mentioned in (E) can be more gen- 
eral quantum critical states. Since some finite tempera- 
ture phases in high T^ superconductors may be described 
by quantum critical states or stable algebraic spin liquids, 
their characterization through quantum order and PSG's 
is useful for describing those finite temperature phases. 

(G) In this paper, we only studied quantum orders 
and topological orders at zero temperature. However, we 
would like point out that topological orders and quan- 
tum orders may also apply to finite temperature systems. 
Quantum effect can be important even at finite tempera- 
tures. In Ref. Q, a dimension index (DI) is introduced 
to characterize the robustness of the ground state degen- 
eracy of a topologically ordered state. We find that if 
DI< 1 topological orders cannot exist at finite tempera- 
ture. However, if DI> 1, topological order can exist at 
finite temperatures and one expect a flnite-temperature 
phase transition without any change of symmetry. Topo- 
logical orders in FQH states have DI^l, and they cannot 
exist at finite temperatures. The topological order in 3D 
superconductors has DI=2. Such a topological order can 
exist at finite temperatures, and we have a continuous 
flnite-temperature superconductor-metal transition that 
do not change any symmetry. 

Although we mainly discussed quantum orders in 2D 
spin systems, the concept of quantum order is not limited 
to 2D spin systems. The concept applies to any quantum 
systems in any dimensions. Actually, a superconductor is 
the simplest example of a state with non trivial quantum 
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order if the dynamical electromagnetic fluctuations are 
included. A superconductor breaks no symmetries and 
cannot be characterized by order parameters. An s-wave 
and a d-wave superconductors, having the same symme- 
try, are distinguished only by their different quantum 
orders. The gapless excitations in a d-wave supercon- 
ductor are not produced by broken symmetries, but by 
quantum orders. We see that a superconductor has many 
properties characteristic of quantum ordered states, and 
it is a quantum ordered state. The quantum orders in 
the superconducting states can also be characterized us- 
ing PSG's. The IGG Q — Z2 ii the superconducting 
state is caused by electron-pair condensation, and the 
IGG Q = Z4 ii the superconducting state is caused by 
four-elcctron-cluster condensation. The different quan- 
tum orders in an s-wave and a d-wave superconductors 
can be distinguished by their different PSG's. The ansatz 
of the s-wave superconductor is invariant under the 90° 
rotation, while the ansatz of the d-wave superconduc- 
tor is invariant under the 90° rotation followed by gauge 
transformations Cj -^ ie^/^Cj. 

It would be interesting to study quantum orders in 
3D systems. In particular, it is interesting to find out 
the quantum order that describes the physical vacuum 
that we all live in. The existence of light - a massless 
excitation - without any sign of spontaneous symmetry 
breaking suggests that our vacuum contains a non-trivial 
quantum order that protect the massless photons. Thus 
quantum order provides an origin of light, pq] 

I would like to thank P. A. Lee, J. Moore, E. Frad- 
kin, W. Rantner, T. Senthil for many helpful discussions. 
This work is supported by NSF Grant No. DMR-97- 
I4I98 and by NSF-MRSEC Grant No. DMR-98-08941. 



APPENDIX A: CLASSIFICATION OF 
PROJECTIVE SYMMETRY GROUPS 



on elements GxT^ and GyTy of the PSG. Here Q is the 
IGG. This condition allows us to determine that there 
are only two different extensions (given by Eq. ( pi] ) and 
Eq. (16^)) for the translation group generated by T^ and 
Ty, if ^ = Z2. 

However, a bigger symmetry groups can have many 
more extensions. In the following we are going to con- 
sider PSG's for the symmetry group generated by two 
translations T^.y, three parity transformations Px,y,xy, 
and the time reversal transformation T. Since transla- 
tions and the time reversal transformation commute we 
have, 

{GxTx)-\GtT)-^GxTxGtT e g 

{GyTy)-\GTT)-^GyTyGTT G Q (A3) 

which reduces to the following two conditions on Gx,y(i) 
and Grii) 

Gx\i)G^Hi)Gxii)GT{i - x) e g 
G-\i)G^\i)Gy{i)GT{i-y) e G (A4) 

Since T-^p-^TPx = 1, T-^P-^TPy = 1, and 
T^^P^J^TPxy ~ 1, one can also show that 

GT\Px{i))Gp'ji)GT{i)Gp^{i) eg 
G^\Py{i))Gpl{i)GT(i)Gp^{i) eg 
GT\Pxy{i))G-pl{i)GT(i)Gp^^{i) eg (A5) 

From the relation between the translations and the par- 
ity transformations, TxP'^TxPx = T'^p-^TyPx = 

T T^^^T P 7^ — 1 p— 17^ p T^^l T^^^T P 

y y ^y^V ~ -'-X -'^y ^xJ^y — ^y J^xy ^xi^xy — 

Tx^PxyTyPxy = 1, wc find that 

{GxTx){G p^Px)^ GxTxGp^Px eg 

{GyTy)-\Gp^Px)-^GyTyGp^Px eg (A6) 



1. General conditions on projective symmetry 
groups 

The transformations in a symmetry group satisfy vari- 
ous algebraic relations so that they form a group. Those 
algebraic relations leads to conditions on the elements 
of the PSG. Solving those conditions for a given sym- 
metry group and a given IGG allows us to find pos- 
sible extensions of the symmetry group, or in another 
word, to find possible PS G's a ssociated with the sym- 
In section IV A , we have seen that the 
1 between translations in x- and 



metry group, 
relation TxTyP^ 



'T- 



y-directions leads to condition 

^x-^ x^y-^ y\^x-^ X ) V^y-^y) — 

GxTxGyTyT^ Gx Ty Gy e g; 



Gx{i)Gy{i - x)G-\i - y)Gy{i)-^ e g 



(Al) 



(A2) 



or 



{GyTy){Gp^Py)-'GyTyGp^Py eg 

{GxTxy\Gp^Py)-^GxTxGp^Py eg (A7) 

(GyTy) {Gp^yPxy) GxTxGp^^Pxy eg 
{GxTx)'\Gp^^Pxy)-'GyTyGp^^Pxy eg (AS) 

Gx{Px{i))Gp'ji + x)Gxii + x)Gp^ii) eg 

G-y\Px{i))G-pl{i)Gy{i)Gp^{i - y) eg (A9) 

Gy{Py{i))G-pl{i + y)Gy{i + y)Gp^ (i) eg 

G-\Py{i))Gplii)Gx{i)Gp^{i - X) eg (AlO) 

G-\Pxy(i))Gpl{x)Gx(i)Gp^^{i - x) eg 
G-\Pxy{i))G-pl(i)Gy(i)Gp^^{i - y) eg (All) 
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We also have P^yP.cP^yP-^ = PyP^P-^P-^ = 1. Thus 

Gp^yPxyGp^PxGp^yPxy{Gp^Py) EG 

Gp^PyGp^P,{Gp^Py)-\Gp^Px)-^ eg (A12) 

which impUes 
Gp.Ai)GpAP.y{i))Gp^^{P:.yP4i))Gpl{i)eg 

Gp^(i)Gp^Py(i))Gpl{PAi))Gpl{i) eg (A13) 
The fact T^ — I leads to condition 

G|(i) e g (A14) 

and P^ = f2 = p2^ = 1 leads to 

GpAi)GpAPxii)) &g 

Gp^{i)Gp^{Py(i))eg 

Gp.^i)Gp^^{P,y{i))eg (A15) 

The above conditions completely determine the PSG's. 
The solutions of the above equations for g = Z2, U{1), 
and SU{2) allow us to obtain PSG's for Z2, U{1), and 
SU{2) spin liquids. However, we would like to point out 
that the above conditions define the so called algebraic 
PSG's, which are som ewhat different from the invariant 
PSG defined in section IV A. More precisely, an algebraic 
PSG is defined for a given IGG and a given symmetry 
group SG. It is a group equipped with a projection P and 
satisfies the following conditions 

IGG cPSG, P{PSG) = SG (A16) 

P{gu)^P{u), ioT any u e PSG and g e IGG 

It is clear that an invariant PSG is always an algebraic 
PSG. However, some algebraic PSG's are not invariant 
PSG's. This is because a generic ansatz Uij that are in- 
variant under an algebraic PSG may be invariant under 
a larger invariant PSG. If we limit ourselves to spin liq- 
uids constructed using Uij, then an algebraic PSG char- 
acterizes a mean-field spin liquid only when it is also an 
invariant PSG at the same time. 

We would like to remark that the definition of invari- 



ant PSG can be generalized. In section IV A, the in- 
variant PSG is defined as a collection of transformations 
that leave an ansatz Uij invariant. More generally, a 
spin liquid is not only characterized by the two-point 

correlation {Uij)ai3 = (i^aiwa^/ but also by many-point 

correlations such as {Uijmn)ai3j\ = V^ai'4^(ij'^\mi^\ 
We may define the generalized invariant PSG as a collec- 
tion of transformations that leave many-point correlation 
invariant. It would be very interesting to see if the gen- 
eralized invariant PSG coincide with the algebraic PSG. 



2. Classification of Z2 projective symmetry groups 

We have seen that there are only two types of Z2 spin 
liquids which have only the translation symmetry. How- 
ever, spin liquids with more symmetries can have more 



types. In this section, we are going to construct all (al- 
gebraic) PSG's associated with the symmetry group gen- 
erated by Tx^yi Px,y,xy, and T for the case g ~ Z2. This 
allows us to obtain a classification of mean-field symmet- 
ric Z2 spin liquids. 

We start with Z2 spin liquids with only translation 
symmetry. First let us add the time reversal symmetry. 
An arbitrary ansatz has the time reversal symmetry if it 
satisfies 



GrTiuij) = 

T{uij) = 



j-ij 



j-ij 



(A17) 



For Z2 spin liquid, the condition Eq. (A4) becomes 

G-\i)G^\i)Gx{i)GT{i -x)^ 77,tT" 
G-\i)G^\i)Gy{i)GT{i - y) - Vytr" (A18) 

where rjxt.yt — ±1- For Z2 spin liquids, Gx.y c>c r" and 
the above four conditions (labeled by rjxt.yt = ±1 ) on 
Gt can be simplified 

G^\i)GT{i-x) = r,xtT° 
G^\i)GT{i-y)=VytT° (A19) 

This leads to four types of Gt labeled by rixt,yt — ±1 

GT{i) ^Vy^V^gr (A20) 

where gx satisfies g^ = ±r°. We see that gr has two 



gauge inequivalent choices gr = t 



Thus the sym- 



metry group generated by Tx.y and T has 2x4x2 = 16 
different extensions (or 16 different PSG's) ii g = Z2- 
There can be (at most) 16 different mean-field Z2 spin 
liquids which have only translation and the time reversal 
symmetries. 

Next let us add three types of parity symmetries. An 
arbitrary ansatz has the parity symmetries if it satisfies 



Gp,Px{uij) =Uij 

Px{i) ={-ix,iy) 

Gp^Py{Uij) ^Uij 

Py{Uij) =Up^(i)^p^(^j) 

Pyii) ={ix,-iy) 

Gp^yPxyyUij) ='>J-ij 

Pxyiuij) EE-Up^^(i)^p_^^(j-) 



(A21) 



(A22) 



(A23) 



For Z2 spin hquids, Eq. (|A9|), Eq. ( [ATc| ), and Eq. ( [AIi| ) 
reduce to 

Gx{Px{i))Gpl{i + x)Gx{i + x)Gp^ (i) =?7,p,T° 

G^\Px{i))Gpl{i)Gy(i)GpAi - y) ^VypxT" (A24) 
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Gy{Pyii))Gpl{i + y)Gy{i + y)Gp^{i) =77,p,r° 

G-'{Pyii))Gplii)G4i)Gp^{i - X) ^rj^pyr" (A25) 

Gy\P,y{i))Gpl{i)GAi)Gp^^{i - x) ^v,p^yT° 
G-\P.y{i))Gpl{i)Gy{i)Gp^^{i - y) ^Vvp^yr" (A26) 

WilGrc ijxpx^xpy^xpxy -'--^ clIlQ 'lypx,ypy,yr}xv Itl. 

We will consider the two cases Eq. ( plj ) and Eq. ( p^ ) 
separately. Fir st w e assu me G xji) — Gyj i) = r". In 
this case, Eq. ( [A2^ ), Eq. ( |A25D , and Eq. ( |A26D can be 
simplified 

Gpl{i + i;)G'p,, (i) =nxpxT° 
Gp'ji)GpAi-y)=f1ypxT'' (A27) 

Gpl{i + y)Gp^{i) =VypyT° 

Gpl (i)Gp, (i - &) =??.pyr° (A28) 

Gpl {i)Gp^, {i - y) ^nypxyr°. (A29) 

GpAi) =VxpxVy'px9P. 

GPy{i) =VxpyVy'py9Py 
Gp^y (i) =VxpxyVypxy9P.y (A30) 

where gp — ±t°, gf, = ±t°, and gf, = ±r°. 77's and 
51's i n the above equation are not independent. From 
Eq. ( [aT^ ), we find that 



We find 



Ixpxy lypxy Ixpxy 'lypxyiJ Pj^ 

- ± 77*^ rf^ 77*" if 

Ixpxy Ivpxv'xpxv'V 



xpxy lypxy Ixpxy lypxy 



(A31) 



which requires T]xpxy = Vypxy = Vpxy From Eq. ( |A13| ) we 
see that 



Ixpxy lypxy Ixpx lypx Ixpxy lypxy Ixpy lypy 
,-1 _ _i_,0 



9P.y9P.9P.y9p^ = ±r 
9Py9P.9pl9pl =±^° 



We find 



and 



VxpyVxpxyVypxVypxy — ^ 

VypyVypxyVxpx^xpxy ^-L 



9P^y9P^9P.y9pl ==±t" 
9Py9P.9pl9pl - ±^° 



(A32) 



(A33) 



(A34) 



Fro m Eg. (A33) we find r]xpx — Vypy ^-nd rjxpy 
Eq. (|A30[ ) becomes 



Vypx 



Gp^ (i) ^rfxpxV%y9P^ 

Gpy [i] =VxpyVxpx9Py 

Gp.yii) ^v;xy9P.y (A35) 

Now the three rjxpx, Ixpy, and "qpxy are independent. 
We note that gauge transformation Wi = rjli^xVwy with 
flwx,wy = ± does not change the form of Gx,y in Eq. ( |6l| ) 
and Eq. (|6j). Thus we can use such gauge transforma- 
tion to further simplify Gp^ ^^^^. We find that the gauge 
transformation Wi = (— )*^ changes rjpxy to —r^pxy Thus 
we can always set 7]pxy = 1- In the following we will 
choose the gauge in which ripxy = 1. 
We also find that Eq. (A5) requires 



Vxt = Vyt 



(A36) 



and 



9t^9p^9t9p, 



9T^9p^9T9Py =±T° 



9T^9pl9T9P:,y 



:±r" 



Thus we only have two types of Grii) 
Grii) =Vt9T 



(A37) 



(A38) 



labeled by 774 = ±1. 

In the following, we will list a ll the gau ge in equivalent 
solutions for g's from Eq. (|A34|) and Eq. (|A37|) . Most of 
them are obtained by setting g's to be one of r^, /i = 
0,1,2,3. 



9Pxy —'' 
9Pxy —'1 

9Pxy =i 
9Pxy =i 
9Pxy =i 

9Pxy —^ 
gpxy =^ 



9P. =T 
9P^ =*T^ 

9P. =r° 
9P. ^iT^ 
9P. =iT^ 

9P. =r° 
9P. =«T^ 



9Py 
9Py 

9Py 
9Py 
9Py 

9Py 
9Py 



9T =T 
9T =T^ 

9t =t^ 
9T =-r' 
9T =t' 

9t ^ir'- 
9T =iT- 



gpxy =T gp^ =^T gPy =it gr ^it 



9T =ir 



gpxy 


—IT 


9P. 


—T 


9Py 


9Pxy 


=iT' 


9P. 


=ZT^ 


9Py 


gpxy 


^ir' 


9P^ 


=ZT^ 


9Py 


gPxy 


=7ri 


9P^ 


=rO 


9Py 


9Pxy 


=iT' 


9P. 


=ZT^ 


9Py 


gpxy 


^ir' 


9P^ 


=ZT^ 


9Py 


gpxy 


=iT' 


9P. 


=it' 


9Py 



gr =iT- 
9T =iT- 



gr =«T 
gr =iT- 
9T ^ir- 
gr =iT- 



(A39) 
(A40) 

(A41) 
(A42) 
(A43) 

(A44) 
(A45) 
(A46) 

(A47) 
(A48) 
(A49) 

(A50) 
(A51) 
(A52) 
(A53) 
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gpxy 


=iri2 


9P. 


gpxy 


^it''' 


9P. 


where 







9Py 
9Py 



0. 



9t =«t^ 



(A54) 
(A55) 



(A56) 



V2 ' V2 ' 

The above 17 solutions, wh en co mbined wi th 8 choices of 
Vt, Vxpx, and rj^py (see Eq. ( |A38| ) and Eq. ( |A35| )), give us 



136 different PSG's for the case G^{i) = Gy{i) = r" 
For Z2 spin hquid with G^{i) ^ (-)'''t", Gy{i) = t° 



Eq. ( |A2^ ), Eq. ( [A25| ), and Eq. ( |A2^ ) can be simphfied 



Gp^{i + x)Gp^ {i) ^r]xpxT° 

Gpl{'>-)Gp, {i - y) =VypxT° 

Gp^ii + y)Gpy{i) =Vypyr° 
Gp^{i)Gp^{i-x) =T]xpyT° 

i-y'Gplii)Gp^^ii - x) ^rj,p,yT° 
i-y''Gpl{i)Gp^^{i - y) =r]yp,yT°. 

The above equations can be solved and we get 

Gp^ (i) ^VxpxVl'pxgp. 

GPy (i) ^VxpyVllygPy 
Gp^,{i) ={-)"""' VxpxyVll, 



(A57) 



(A58) 



(A59) 



xpxy lypxy 



9P. 



(A60) 



Eq. ( A15 ) and Eq. ( [A13[ ) still leads to rjxpxy — Vypxy and 
Eq. (A32). Thus 



GP^ (i) =VxpxVxpy9P^ 
Gpy{i) =VxpyVx'px9Py 

GP.M <-)'-''' 9P.y 



(A61) 



Eq. (A34) and Eq. (A37) are still valid here, which lead 



to the same choices for g's. For the case {—Y^Gxii) = 
Gy{i)_= T°, the 17 choices of .g's in Eq. ( [A39| ) to 
Eq. ( |A55 ) , when combined with 8 choices of rjt , Vxvx , and 
Vxpy agai n giv e us 136 different PSG's through Eq. ( |A38| ) 



and Eq. (|A6l|). 

Now we would like to consider which of the translation 
symmetric ansatz in Eq. (p3) or Eq. (p4h have the parity 
and the time reversal symmetries. We note that three 
parity symmetries also imply the 90° rotation symmetry. 

After two parity transformations Px and Pj,, we find 
Um in Eq. (|6^) or Eq. (p3) satisfies 



u_ 



VxpyVxvx9Py9P.Um.9p^ 9 P^ 



Ixpy ixpxb 

After the time reversal transformation, we have 



(A62) 
(A63) 



Thus Um = u'^T^ must satisfy 

7;" — n if r?"* r)"^ — 1 or r?™ — 1 
"■m ~"J ^^ 'Ixpy'lxpx — r 01 ;/j — 1 

uLt^ =vTpyr]Tpxgpygp.UmT^gplgpl 



(A64) 



in order to have the parity and the time reversal symme- 
tries. We see that uij — Q ii qt = t'^ and rjt — 1 and 
""L = 0; ' = 1: 2, 3, if rjt = VxpxVxpy and gr = ±gPygp^. 
There are 2 x 6 x 4 = 48 PSG's with gr = t° and 
rjt = 1. There are 2 x 1 x 2 = 4 PSG's with gr ^ t°, 
gr = ±gpygp^, Vt = VxpxVxpy, and r]t = 1. There are 
2 x 6 X 2 = 24 PSG's with gr = ±gPygp^, rjt = rjxpxVxpy, 
and rjt = —1. Since the ansatz that are invariant under 
the above PSG's have u^ = 0, those ansatz are actually 
invariant under larger PSG's with IGG equal or larger 
than SU{2). Thus there are at most 272 - 48 - 4 - 24 = 
196 different mean-field Z2 spin liquids that can be con- 
structed form Uij. 

For ansatz of type Eq. m3) the parity symmetries also 
require that 



up^(m-) =v7plv7p''ygp^umgp^ 

UPyim) ^V^pyV^lgp^UrngPy 
Up^yim) ^VTxy9plyUmgP^y 



(A65) 



For ansatz of type Eq. (|6J) the parity symmetries require 
that 

Up^(m) ^r]7px'n'^py9plum.9P, 
-1„ 



UPy(m) =VxpyVxp'x9PyUm9Py 
"Pxy(m) 



I ) Ipxygp^y ""rngP^. 



(A66) 



For each choice of g's and Ty's, Eq. (A64), Eq. ( A6E ), and 
Eq. (A66) allow us to construct Z2 symmetric ansatz Uij. 



3. Classification of U{1) projective symmetry 
groups 

In this section we will use PSG to classify quantum 
orders in U{1) spin liquids by finding the PSG with IGG 
Q = U{1). First we note that elements in the U{1) IGG 
must have a form e^^'"^'^ where Vi is a site dependent 
vector. We can always choose a gauge such that Vi all 
point to the same direction, say, the r^ direction. We 
will call this gauge canonical gauge. We also find that 
l^il must be independent of i in order for Uij to be non- 
zero and invariant under the IGG. Thus, in the canonical 
gauge, IGG has a form 



Q^{ 



e'''^ \ee [0,2^)} 



and the ansatz Uij has a form 



U^T^'+uljT^ 



(A67) 



(A68) 
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We see that the flux through any loops is in the r^ direc- 
tion. Due to the translation symmetry of the ansatz, the 
absolute value of the flux must be translation invariant, 
but the sign may change as we translate the loops. Thus, 
the loop operator have a form 



Pc, = irT'Pc^^oirT' 



(A69) 



where n^ = 0, 1 and Ci is loop with base point i. Here 
the two loops Ci and Cj are related by a translation and 
have the same shape. Now we can choose a different 
gauge by making a gauge transformation Wi = (ir^)"'. 
In the new gauge we have 



^ZJ 



-^vl^"i-"3■„?.^o 



<^r^) 



^r^' + (zr^)"'+"^«^r 



Pc, ^Pc^.o 



g^{e'(-y'O-"\0e[O,27r)} 



(A70) 



Since the loop operators are uniform, we will call the new 
gauge uniform gauge. 

Let us first work in the uniform gauge. From the trans- 
lation invariance of Pd 

Pc^ =G^{i)Pc,_M-\i) = G^{i)Pc^G-\i) 

Pc^ ^Gy{i)Pc,_,G-\i) = Gy{i)Pc,G-\i) (A71) 

we find that G^.y have a form 



Gyii) =g3{ey{i)) 



(A72) 



Now we switch to the canonical gauge. We note that 
a gauge transformation that keep an ansatz to have the 
form in the canonical gauge Eq. (A68) must have a one 
the following two forms 



Wi=g3{9{i)) 
Wi^iT^gsiOii)) 



(A73) 
(A74) 



if we require that Uij ^ 0. (More precisely, we require 
that any two points on the lattice can be connected by 
several non-zero Ujj's. We will call such an ansatz con- 
nected.) Thus for spin liquids with connected Uij, Gx.y 
must take one of the above two forms in the canonical 
gauge, since Gx.y are two special gauge transformations. 



From Eq. (A72), we find that Gx,y have a form 



Gy{i)^HTY^gsi0y{imT'r-^ 

Gx{i)M-^rTg3{Uimr')"'-' 



(A75) 



in the canonical gauge. Thus ni can only be one of the 
following four choices: rii — 0, m — {1 — (— )*)/2, Ui = 
(1 - {-Y'')/2, and m = {1- {-y^)/2. In these four 
cases, Gx,y take one of the above two forms and Uij can 
be connected. 

Let us consider those cases in turn. We will work in 
the canonical gauge. When rij — 0, Gx.y have a form 

G.(i)=53(e.W), Gy{i) = g^{ey{i)) 



Since the gauge transformation Wi — gsiOi) keep an 
ansatz and its PSG in the canonical gauge, we can use 
such kind of gauge transformation to simplify Gx. y by 
setting 0y{i) = and Oxiiy = 0,ix) = 0. Now Eq. (^) 
takes a form 



Gx{i)Gxii-y)-'^g3i^) (A77) 

for a constant ip. This allows us to obtain 

Gx{i) = g3{iy^ + Ox), Gy{i) = gsiOy). (A78) 

The translation symmetric ansatz has a form 

UiA+m = ipmgsi-myixip + (j>) (A79) 

where pm > 0. The above ansatz describes particle hop- 
ping in uniform "magnetic field" with e^'^'^ flux per pla- 
quette. In this case ip/ir should be a rational number 
(ys/TT ~ p/q (between and 1) so that the ansatz can be 
put on a finite lattice. Thus (/j/tt should be viewed as a 
discrete label and different rational numbers between 
and 1 will gives rise to different type of spin liquids. 
When ni = (1 + (— )'')/2, Gx,y have a form 

Gx{i) - g3{ex{i))iT\ Gyii) = g3i9y{i))iT' (A80) 

Again we can use the gauge transformation Wi = gz[9i) 
to simplify G x.y by setting 9y{i) — and 6x{iy = Q,ix) = 
0. Now Eq. (K3) takes a form 

Gx{i)T^Gx{i - y)-Vi -ff3(^) (A81) 

or 

9x{i) + 9x{i - y) ^ V (A82) 

for a constant (p. This allows us to obtain 

Gx{i) = gsiH'^Ar + 9x)ir\ Gy{i) = g3{9y)iT\ 

(A83) 

where (ptr G [0,7r). A gauge transformation Wi = 
5'3(— (— )'''0*r/2) change the above to 

Gxii) = g3{9x)iT\ Gyii) = g3i9y)iT\ (A84) 

The translation symmetric ansatz has a form 

Ui,i+m. = iPmgz ((-)*0m) (A85) 

When ni = (1 + (—)'=" )/2, Gx,y have a form 
Gx[i) = gz[9x{i))iT\ Gyii) = gzi9yii)) 



After using a gauge transformation Wi 



(A86) 

g3i9i) to sim- 
plify Gx,y by setting 9yii) — 0, Eq. ( |A2| ) takes a form 

Gxii)Gxii - y)-^ ^g3ip) (A87) 

or 

(A76) 9xii) - 9xii - y) = ^ (A88) 
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for a constant Lp. This allows us to obtain 

Gx{i) = 9z{iyV + dx)iT\ Gy{i) = gsiOy). (A89) 

where (/p e [0, tt). A gauge transformation Wi = 
g'^[—iyif/2) change the above to 

G.(i) = g^{e^)iT\ Gy{i) = g^iOy). (A90) 

The translation symmetric ansatz has a form 

Ui,i+m=ipmg3{{-y"4>m) (A91) 



To summarize, Eq. ( |A7S| ), Eq. ( |A84| ) and Eq. ( |A9q ) are 
the most gene ral tra nslati on PSG's t hat a llow non-zero 
Uij. Eq. ( JA79D , Eq. ( |A85D , and Eq. ( |A9l| ) are the most 
general translation symmetric mean-field ansatz for U{1) 
spin liquids. 

Next, we would like to include more sy mmet ries. We 
first consider the translation PSG in Eq. (A7S). When 
(y9 = the translation PSG has a form 



Gx(i)=g3(e.), G,(i)=.93(e,). 



(A92) 



The corresponding spin liquids will be called type UlA 
spin liquids. When (p = tt the translation PSG has a 
form 



G,(i) = (-)'«g3(^.), Gy{i)^g^{ey). 



(A93) 



and the corresponding spin liquids will be called type 
UIB spin liquids. For other value of Lp, we will call the 
corresponding spin liquids type Ul™ spin liquids, where 
m/n = ip l-K m od 1. The translation PSG's Eq. ( [A84| ) 
and Eq. ( A90 ) will correspond to type UlC and type 
UID spin liquids respectively. 

Let us first consider the type UlA spin liquids. To add 
the time reversal symmetry, w e no te that, for type UlA 
spin liquid, the condition Eq. (|A4|) becomes 



G^i(i)GT(i-y)ef/(l) (A94) 

This leads to two types of Gt 

Gt = 53(i • Vi + Ot), igaii ■ ^t + St^ (A95) 

Since T^ = 1 and G^ G U{1), the above becomes 

Gt - Vl-tf^l\93{0t), ^gsii ■ Vt + dtV (A96) 

To add three types of parity sy mme tries, we note that, 
for t ype UlA spin hquids, Eq. (|A9|), Eq. ( |A10D , and 
Eq. ( All| ) reduce to 



Gpl{i + x)GpSi)eU{l) 
Gp'ji)GpJi-y)eU{l) 

Gplii + y)Gp^{i)eU{l) 
GpHi)Gp^ii-x)eUil) 



(A97) 



(A98) 



Gpl{i)Gp^^{i-y)eU{l). (A99) 

We find that Gp p p can have the fohowinff forms 

Gp^ =53(i • Vpx + ^pt), m{i ■ ^px + ^px)t^ 
Gpy =93(1 ■ tPpy + 9py), igsii ■ fpy + epy)T^ 

Gp^y =93{i ■ <-Ppxy + 9pxy), igsii ■ <Ppxy + Spxy)T^ 

(AlOO) 

Note that the gauge transformation Wi — gsii ■ 0) does 
not change Gx.y Thus, we can use it to simplify Gp^^p 
and get 

Gp^ ^gsiiy'Ppx + 6px): igsi'ix'Ppx + Opx)T^ 

Gpy ^gsiixVpy + Opy), igaiiyVpy + (^py)T^ (AlOl) 

Gp^y =g3ii ■ ^pxy + Qpxy), igsii ' iPpxy + 9pxy)T^ 



From the condition Eq. (A15), we find that 
Gp^ ^Vllx93i.epx), iVxpx93{0pxy 

GPy =Vxpy93{dpy), iVyhygsiOpyy (A102) 

Gp^y =53 ((ia; — iy)'Ppxy + (^pxy), 433 ((*£c + iy)(Ppxy + Spxy)T 



G p p p^„.T should also satisfy Eq. ( A13| ) and 
Eq. (A5). We find Gp_^ can be obtained from Gp^,p 
through Eq. ( |A13| ) and Gt from Gp^^p^^p^^ through 
Eq. (|A5|) . This leads to the following 24 sets of solutions 



Gp. 

Gp.y 



--vlixg^iSp 



Gi 



lypxi)-3\'^pxh '-'Py — VypxgsyOpy 

=g3{0pxy), g3{0pxy)iT 



Gt =493{0t)\n,=-i, vig3{0tW 



(A103) 



and 



Gp. =Vxpxg3{0px)iT^ , Gpy = rfypxg3{0py)iT^ 

Gp.y =g3{Qpxy), g3{Spxy)iT 



Gt =4g3{et)W- 



Tj;g3[0t)iT^ 



(A104) 



Whe n co mbined with the type UlA translation PSG 
Eq. ( A92 ) , the above 24 sets of solutions give us 24 dif- 
ferent PSG's. A labeling scheme of the above PSG's is 
given below Eq. (|70|). 

Now let us consider the form of ansatz that is invari- 
ant under the above PSG's. The translation symmetry 
requires that 



ulr' 



"L^^ 



(A105) 



The 180° rotation symmetry requires that, for Gp^ 

'riypxg3{Qpx),Gp^ = nl"pxg3{Qpy), 



Vypx^-rn 



77"* u'' 

lynx ^r, 



lypx^-rn — 'lypx"'m, (AlUoj 

and for Gp^ = Vxpxg3{0px)iT^,Gp^ = Vxpxg3{dpy)iT^, 



Ur 



Vx 



^U- 



(A107) 
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The time reversal symmetry requires that, for Gt 



<-r^r 



ioT Gt = vi93{0t)iT^ 



ur^ = -vruLr° 



vrnLr\ 



(A108) 



(A109) 



The above properties of Uij can be obtained after reahz- 
ing that itij ca n be o btained by combining Uij with the 
ansatz in Eq. (A114). We see that the mapping has a 
meaning of adding 7r-flux to each p laquette. 

Using the mapping Eq. ( All2| ) and the results for 
the type UlA symmetric spin liquids, we find that the 
type UIB symmetric spin liquids are also classified by 24 
PSG's. Gp^^p p^ T of those PSG's are given by 



We find the following 8 sets of ansatz that give 
rise to U{1) symmetric spin liquids: U1A00[0, 1]1 and 
U1A11[0,1]1, 



'^i.i-\-m '^~^rn^ ' 

UlAO?i[0, l]x and UlAa;l[0, l]x, 

Ui,i+m. =u'^T° + ul^T 
..0 



(Alio) 



0, if m = even 
0, if m = odd. 



(Alll) 



Other 16 PSG's lead to SU{2) spin liquids and can be 
dropped. 

To obtain the PSG's for type UIB symmetric spin 
liquids, we would like to first prove a general theorem. 
Given a PSG generated by Gx,y,T and Gp^^p^^p^^, the 
following generators 

G^i) H-y^Gxii), Gy{i) =Gy{i), 

GpJi)^GpJi), Gp^ii) ^Gp^ii), 

Gp^^W-l-^'^Gp^.W, GT{i)=GT{i). (A112) 

generate a new PSG. The new PSG has the same IGG 
and is an extension of the same symmetry group as the 
original PSG. If an ansatz Uij is described by a PSG 
{Gx,y,Ti Gp^,p^,p^^), a new ansatz described by PSG 
iGx,y,T, Gp^^p^^p^J can be constructed 



<') 



Uv 



'Uij, for {jx - ix)ijy - iy) 



Uij =0, for [jx ~ ix){iy - iy) = odd 



: even 
(A113) 



The new ansatz Uij has the same symmetry and the same 
gauge structure as Uij. 

To obtain the above result, we note that the following 
ansatz 



— (\(jy-iy)i!^-r'^ 



Uij =(-) 



for {jx - ix){jy -iy) ^ even 



Uij =0, for {jx - ix){jy -iy) = odd 



(AIM) 



has all the translation, parity, and the time reversal sym- 
metries and has an SU{2) invariant gauge group. The 
PSG of the ansatz has a subgroup 



Gx{i) ={-r^r\ Gy{i) 

GpAi)^r\ Gp^ii) 

G p. Ji) ={-)'-''' r' 



Gp^ ^Vyi,x93{0px), Gp^ ^Vypx93iOpy), 
{-T'-'^Gp^^ ^g^iOpxy), 9z{0pxy)iT^ 

Gt =7i93{et)W=-i, vi9M)^r\ (A116) 



and 



Gp^ =Vxpx93{&px)ir^ , Gp^ = Vx'px93{Spy)iT^ , 



i-y-'^Gp^^ ^dsiOpxy), 93{0pxy)iT\ 

Gt =vi93{0t)\vt=-u vt93{0t)tT\ 



(A117) 



A labeling scheme of the above PSG's is given below 
Eq. (0). 

Next we consider the form of ansatz that is invariant 
under the above PSG's. The translation symmetry re- 
quires that 



Wi,i+m — (, ) 



'!'7/„_ = 



(A118) 



The 180° rotation symmetry requires that, for Gp^ 

Vvpxg3{0px),Gp^ ^Vl'px93{0py), 



r?"* u 



mxmj,„t 
m 



%px{- }""""" u'^ (A119) 

and for Gp^ = ■n%x93{0px)iT^ ,Gp^ = Vxpxg3{dpy)iT\ 

um = vTpx^-m = ^^. (-)"-"" "L (A120) 



The time reversal symmetry requires that, for Gt = 
r]ig3{0t)\r,,=-i 



-{-)^Ur 



for Gt = 493{0t)iT^ 



Urr. = -Vru'mr° 



rirnLr\ 



(A121) 



(A122) 



We find the following 8 sets of ansatz that give 
rise to U{\) symmetric spin liquids: U1B00[0, 1]1 and 
U1B11[0, 1]1, 



:.z-\-m 



-U^ T^ 



^m =0: if ^^x = odd and niy — odd. (A123) 
UlBOn[0, l]x and UlBa;l[0, l]a;, 



=w° r" + utr^ 



(A115) 



Ui^i-\-m 

u^ =0, if m = even 

-0, if m,x = odd or niy 



m 

,3 



odd . (A124) 
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Other 16 PSG's lead to SU{2) spin liquids and can be 
dropped. 

Next we consider the type UlC spin liquids. To add 
the t ime reversal symmetry, we note that the condition 
Eq. (A4) becomes 



T^G:r\'i'yGT{i-x)e U{1) 
T^G^\iyGT{i-y)eU{l) 

This leads to the following Gt 

GT=53((-)V(i)), g,ii-rH^))^T\ 
where 4>{i) satisfies 

(j>{i + x) ~ (j){i) —{—yifi mod 27r 
(j){i + y) - (/)(i) =(-)V2 mod 27r 



(A125) 

(A126) 

(A127) 



The solution exist only for two cases where tpi — (f2 
mod tt: 



[i + x)- (j){i) = - i-y29t mod 2n 
(i + y)- (j){i) = - (-)*26lt mod 27r 



(A128) 



and 



4>{i + x)~ (j){i) = - i-y{20t + n) mod 2tt 

(/.(i + y) - 0(i) = - (-)^20t mod 27r (A129) 

The two solutions are given by 

0(i) ^^t + i-yot 

cbii)^^t + i-yOt+i,T: (A130) 

Thus Gt can take the following four forms 

GT=ff3((-)Vt + ^*), (A131) 

g3{ixTr + hyft + Ot), 

gz{{-yvt + OtyT\ 

g3{ix7r+{-yipt + et)iT^. 
1 and G^ G U{1), the above becomes 



Since T^ 



Gt ^VxtVvf93{0t), 

gsii-yipt + otyr^ 

g3{ixTT+{-y(pt+et)iT^. 



(A132) 



*xt,yt 



±1. 



where r]^ 

To add three types of parity symmetries, we note that, 
for t ype UlC spin liquids, Eq. ( |A9| ), Eq. ( |A10[ ), and 
Eq. ( All| ) reduce to 



T^Gpl{i + xyGpJi)eU{l), 
T'Gp^iiyGpAi-y)eU{i), 

T'Gpl{i + yyGp^{i)eU{l), 
T'GpHiyGp^{i-x)eUil), 



(A133) 



(A134) 



T'GpliiyGp^^ii-x)eU{l), 
T'GpliiyGp^^ii-y)eUil). 



(A135) 



After a calculation similar to that for Gt, we find that 
Gp^.p ,p^ can have the following forms 



Gp, =ff3((-)Vpx +Op,j:), 

gsiixT^ + {-y^px + Opx), 

gsiixTT + {-y-^px + dpx)iT^; 



Gpy ^gaii-y^py+Opy), 

53(ix7r+ (-)Vpy +Spy), 
g3{{-yVpy+0py)iT^, 
53(«x7r + (-)Vpy + Gpy)iT'^\ 



Gf 



=53((-)Vpa;i/ + (^pxy), 
gsi^xT^ + {-y^pxy + Opxy), 
gaii-y^pxy +dpxy)iT^, 
gaiixT^ + {-y<Ppxy + Opxy)iT'^\ 



From the condition Eq. ( A15 ), we find that 
Gp, ^rixpxVy'pxgai^px), 

gsii-Yfpx + dpx)iT^, 

gs.{ixi^ + {-y^px + 6px)iT^\ 
Gpy --vilyVilygaiOpx), 

53(ix7r + i-y^Ppy + Opy)iT^; 



G, 



=53 i 



.), 



g3{ixT^+{~y^+dpxy): 

gzii-y^pxy + 



^^T'■. 



(A136) 



(A137) 



(A138) 



(A139) 



(A140) 



(A141) 



The first condition in Eq. (A13) requires Gpx,Py should 
have the same number of -PT The second condition in 



Eq. (AlS) further requires that ipp^ = (ppy mod 7r/2. We 
note that Gx.y in Eq. (A84) are invariant under gauge 
transformation Wi = g3[[—y<f))- Using such a gauge 
transformation, we can set (ppx = and ippy = mod 
7r/2. This leads to 



GP^ =Vl'pxVy'pxg3{0px), 

Vxpx93{0px)iT'^; 

Gpy ^Vl'py'fly'pygsi&px), 

vllxVlpgsidpyW; 



(A142) 



(A143) 
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Gp^y —g^iOpxy), 

93{ixTr+{-yj + 9pxy), 

gsii-T'^pxy + Opxy)iT^; 



(A144) 



W e find Gp^^ can be determined from Gp^.p^ tlirougli 
Eq. (|A13[) and Gt from Gp^pp^^ tlirougli Eq. (^). 
Thus from Eq. (A13) and Eq. (A[), we find the following 
60 sets of solutions for Gp p p t- 



<■ X-,-' y^->- xy J 



Gp^ —Vl"pxVypx93i(^px), Gp^ 

Jpxy), 



Vypx Vxpx93\"py)j 

Gt ^vt93(0t)U=-i, v'pxy93(0tW\ (A145) 



^Pxy Ipxy93vtpxy A 



xpxilo\"pxJ''i , GPy — 'nxpx^pxy93\Ppy}1T , 

pxy'7 + "pxy): 

Gt ^vi93{0t)U=-i, v'pxvVi93{0t)^r\ (A146) 



Gp^ ^Vxpx93{0px)iT^ , 
Gp^y =r]'ply93in^ 



Gp^ =VxpxVy'px93{0px) 
Gp^y ^93{0pxy)iT^ , 

Gt ^493{0t)W=-i- 



Gf 



Vypx Vxpx93\"py)j 



(A147) 



Gf 



=VxpxVy'px93{0px), Gpy = VypxVl'pxgsiOpy), 



Gp^y ^93{vlxy^ + Opxy)iT'^, 

Gt =Vpxv493{9t)iT\ 



(A148) 



Gp^ =lxpx93i^px)iT^ , Gpy = rfypxVlxv93{0py)iT^ , 

Gp^y ^93ivlxyJ + dpxy)iT\ 

Gt ^vi93{0t)W=^i, vWply93{0t)iT\ (A149) 



Whe n co mbined with the type UlC translation PSG 
Eq. (A84), the above 60 sets of solutions give us 60 dif- 
ferent type UlC PSG's. A labeling scheme of the above 
PSG's is given below Eq. ^). 

Now let us consider the form of ansatz that is invari- 
ant under the above type UlC PSG's. The translation 
symmetry requires that 



UiA+m = u°^T° + (-)*uj 



(A150) 



For PSG's UlC[00,nn][0,n, l]n, U1C11[0, 1]", and 
UlCa;l[n, x]n, the 180° rotation symmetry requires that, 



"L = i-r^i.- 



The time reversal symmetry requires that, 



yO ^ „c„^"^ 3 



i~ru. 



m„,3 



The ansatz have a form 



Tn^3 3 



u^ =0, if m = even. 

which describe SU{2) spin liquids. 
For PSG's 



(A151) 



UlC[n0,0n][0,n, l]n, UlCll[n,x]ri, 

UlCa;l[0,l]n, (A152) 

the 180° rotation symmetry requires that 



3 =u3 

m' "'m "'m,' 



{-rul^. u: 



The time reversal symmetry requires that 



^"m = -i-rw. 



m„.Q ,,3 



-(-ru: 



m' 



The ansatz have a form 



"^2,1+7 



--ulr' + i-rul 



u^ =0, if m = even. 



For PSG's 

UlC[00,nn][0,l]l, U1C11[0, 1]1, 

the 180° rotation symmetry requires that, 



(A153) 



(A154) 



m^ 3 



^mi m 



{~ru. 



The time reversal symmetry requires that. 



^rm "m 



The ansatz have a form 



?Ti^,3 ^3 



=0, if m = odd. 



(A155) 



The ansatz gives rise to U{\) x U{1) spin liquids since 
Uij only connect points within two different sublattices. 
For PSG's 



UlC[nO,On][0,l]l, UlCa;l[0, 1]1, 

the 180° rotation symmetry requires that 



(A156) 



"m V J "m' ""m "r 



The time reversal symmetry requires that 



qfi — —7;*^ 11'^ = 7/'^ 

"m "•mJ "m "m- 



The ansatz have a form 

For PSG's 

UlC[00,nn]la;, U1C11[0, l]^;, 



(A157) 
(A158) 
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the 180° rotation symmetry requires that, 



i-ru 



rrt^,3 



The time reversal symmetry requires that, 



"i;. = -(-rc, uL^i-ru: 






The ansatz have a form 



uL, —Q, '\i m ~ odd. 



(A159) 



The ansatz gives rise to U{\) x U{1) spin hquids since 
Uij only connect points within two different sublattices. 
For PSG's 



UlC[?iO,On]lx, UlCxl[0, l]a;, 

the 180° rotation symmetry requires that 

"m V / "ml "m "m ■ 

The time reversal symmetry requires that 



y!L^-{-rui„ ui,^i-rui,. 



(A160) 



The ansatz have a form 



m^,3 ^3 



""m, — O7 ii m — even 
u^ =0, if m = odd. 

For PSG's 

UlC[00,nri][ri,a;]l, UlCa;l[n, x]l, 

the 180° rotation symmetry requires that, 



(A161) 
(A162) 



771^,3 

771" 



The time reversal symmetry requires that 



The ansatz have a form 



UiA-i 



<~ruir^ 



uL =0, 



^n^ 3 3 

if TOj: = odd or tti^ = odd. 



(A163) 



The ansatz gives rise to (C/(l))* spin liquids since Uij 
only connect points within four different sublattices. 
For PSG's 



UlC[?iO,On][n,a;]l, UlCll[n, a;]l, 

the 180° rotation symmetry requires that. 



(A164) 



\ I '-'"im "■■ 



3 
m 



The time reversal symmetry requires that, 



"■m ~ \ I "■rrn 



(-) 



771 :c^, 3 



The ansatz have a form 

^rra =0, if nix — even or niy = odd 

ul^ =0, if nij, = odd. (A165) 

For PSG's 

UlC[00,nn]a;a;, \]lCxl[n,x\x, (AI66) 

the 180° rotation symmetry requires that. 



(-)-u 



m.^,3 
m ■ 



The time reversal symmetry requires that, 



ui. = -{-r^< 



The ansatz have a form 

U^ =0, if TTlj, = 



(-)'"«u 



7^71,^,3 

m' 



odd or rrij, = odd. (A167) 



The ansatz gives rise to (t/(l))'* spin liquids since Uij 
only connect points within four different sublattices. 
For PSG's 



UlC[nO,On]ra, UlCll[n,x]x, 

the 180° rotation symmetry requires that. 

The time reversal symmetry requires that. 



(A168) 






The ansatz have a form 



(-)"»u: 



777 1,^,3 

77^" 



^rn =0; ^^ "^a: = ^dd Or niy — 



=0, if TO,, 



odd. 



(A169) 



The type UID spin liquids always break the parity gen- 
erated by Pxy and Eq. ( All ) cannot be satisfied. Thus 
there is no type UID symmetric spin liquid. 

Last, let us consider the type Ul™ spin liquids. Instead 
of finding a classification of Ul™ spin liquids, here, we 
will just consider the following example: 



^i,i+x \^ •) ^i.i+y 



I "'^'^ ■ \ 3 

X9z\ '>'x)t 

n 



(A170) 



One can check that the above ansatz describes a sym- 
metric spin liquid. Its PSG is given by 

77777- 
Gx =g-i{-—iy + ^x), Gy = 53(^y), (A171) 
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The form of Gx,y tells us that Eq. (A17C) indeed de- 



scribes a UIJJ' spin liquid. Using the labeling scheme for 
the U1[A,B,C] PSG's, we can label the above PSG by 
W^xxxn. From the above example, we see that there 
are infinite different spin liquids of type Ul™ , at least 
one for each rational number m/n between and 1. 

In summary, we find 8 type UlA, 8 type UIB, and 30 
type UlC symmetric C/(l) spin liquids. But there is an 
infinite number of type Ul™ spin liquids. 



4. Classification of SU{2) projective symmetry 
groups 

In this section we will use PSG to classify quantum 
orders in mean-field symmetric SU{2) spin liquids. We 
need to find the extensions of the symmetry transforma- 
tions when IGG G = SU{2). First, we assume that, for 
a SU{2) spin liquid, we can always choose a gauge such 
that Uij has a form 



j-ij 






(A172) 



We will call this gauge canonical gauge. In the canonical 
gauge, IGG has a form Q — SU{2). Here we will only 
consider spin liquids described by non-zero Uij. In this 
case the gauge transformations that keep Uij to have the 
form in the canonical gauge are given by 

Wi = r^{i)g (A173) 

where r]{i) = ±1 for each i and g £ SU{2). The gauge 
transformations G^^y 
take the above form: 



transformations G^^y associated with the translation also 



Gxii) ^'nx{i)gx 

Gyii) ^Vyii)9y 



(A174) 



Note that gauge transformation Wi = ri{i)T^ still keep 
Uij in the canonical gauge. So we can use such gauge 
transformation to simphfy Gxii) and Gy{i) (see Eq. (|49|)) 
and get 

Gxii) =gx 

Gy{i) ^iiy{i)gy (A175) 

with riy{ix — 0, ij,) — 1. Now Eq. ( |A^ ) takes a form 

r]y{i-x)ny{i)&SU{2) (A176) 

We find that there are only two different PSG's for trans- 
lation symmetric ansatz 

Gx{i) ^gx Gyii) =gy (A177) 

Gx{i) ^gx Gy{i) H-y-gy (A178) 

The two PSG's lead to the following two translation sym- 
metric ansatz 



^i.i+m — ^771^ 



^i^i-^rm — t, ) 



m^iy 



(A179) 
(A180) 



Next we will consider the case Gx{i) = gx and Gy{i) = 
gy and add more sy mme tries. First let us add the three 
parities Px.,y,xy Eq. ( |A9| ), Eq. ( |A10D , and Eq. ( [All| ) can 
be simplified 

G-pl{i + x)GpJi)eg 
Gp^{i)Gp^{i-y)eg (A181) 

Gpl{i + y)Gp^{i)eg 
Gplii)GpJi-x)eg (A182) 

Gpl{i)Gp^^{i-x)eg 
Gp'ji)Gp^Ji-y)eg. (A183) 

Gp^ (i) =VxpxVy'px9P. 
GPy (i) =VxpyVy'py9Py 
Gp^y{i)=Vx'pxyVypxy9P.y (A184) 



We find 



where gp^ G 5C/(2), gp^ e 5[/(2), and gp^^ G 5(7(2). 
77's in the above equation are not independent. From 
Eq. (|A15D, we find that 



,9k. e g 



(A185) 



which requires rjxpxy 
see that 



Ixpxy lypxy Ixpxy lypxy^ 

rjypxy = itpxy From Eq. ( |A13|) we 



Ixpxy lypxy Ixpx typx Ixpxy lypxy Ixpy lypy 

9P^y9P^9P^y9pl e Q (A186) 



We find 



^xpyf]xpxyf]ypx^ypxy — -L 
^ypyfjypxy'^xpx'^xpxy — ^ 



(A187) 



When combined with rixnxv = Vvn xm we see that rjxp; 



riypy and rjxpy = ijypx- Eq. (|A184[) becomes 



GpAi) 
Gpyii) 

Gp.yii) 



~ Ixpx 'Ixpy yPx 
~ Ixpy Ixpx yPy 
~ IpxyclPxy 



(A188) 



Now the three rjxpxy V xr>y : and rjpxy are independent. 

Similarly from Eq. (A4) and Eq. (A5) we find there are 
only two types of Gt'- 



Grii) =vl9T 



(A189) 



labeled by ryt = ±1. We note that rjt = I implies Uij = 
for the SU{2) spin liquids. Thus we can only choose 



Gxii) = (-)V 



(A190) 
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Eq. ( |A177l ), Eq. ( |A188| ) and Eq. (|A190D give us PSG's 
for symmetric SU(2) spin liquids. They are labeled by 
Vxpx,xpy,pxy Again we can use the gauge transformation 
Wi — (— )*^ to set rjpxy = 1. Thus there are only four 
PSG's labeled by rixpx,xpy 

Now we consider the case Gxji) — g x and Gy{i) — 
{-Y'^gy Eq. (^), Eq. ( |A10D , and Eq. ( |AllD have the 
form 



Gp;(i + y)GpJi)Gg 
G-pl{i)Gp^{i~x)eg 

i-y^Gpi{i)Gp^^{i-x)eg 

i-y-Gpl(i)Gp^^{i-y)eg. 



(A191) 



(A192) 



(A193) 



We find 



Gp^i) 



Ixpx lypx yPx 
~VxpyVypy9Py 



ixpyiypy 

Gp^y{i) ={-)"''"' vllxyn 



pxy lypxy 



gp.y 



(A194) 



where gp^ e 5C/(2), gp^ e 5C/(2), and gp^^ e 5[/(2). 77's 
in the above equation are not independent. We find that 

Vxpxy Vypxy ^ Vpxy^ Vxpx Vypy aUQ Ijxpy Vypx- 

After setting r]j,x v — 1 through gauge transformation 
Wi = (-)*-, Eq. ( ^194| ) becomes 



GP^ (i) =VxpxVx'pygP^ 

GPy (i) ^VxpyVx'pxQPy 

Gp.Ai)H-y-''^9p.y 



(A195) 



and the two rj^px and rjxpy are independent. The gauge 
transformation associated with the time reversal trans- 



formation is stiU gi ven b y Eq. ( |A190[ ). Eq. ( [A178[ ), 
Eq. ( A195 ) and Eq. ( A190 ) give us PSG's for symmetric 
SU{2) spin hquids. 

We see there are total of 2 x 2^ = 8 different SU{2) 
PSG's. A labeling scheme of those 8 SU{2) is given below 
Eq. (|l|). Those S'C/(2) PSG's are algebraic PSG's. In the 
following we will see which of them lead to symmetric 
SU{2) spin hquids. 

First we co nsider which of the translation symmetric 
ansatz in Eq. (A17£) have the parity and the time rever- 
sal symmetries. After two parity transformations Px and 



Py, we find Ur, 



in Eq. (A179) satisfies 



After the time reversal transformation, we have 

-Urrv - i-)"^gTUm9T' (A198) 

Thus the SU2AnO and SU2A0n symmetric ansatz have 
a form 



Urr,. =M™T° 



=0, if m = 



(A199) 



The other two PSG's SU2A[00, nn] leads to vanishing Uij 
and should be dropped. 



For the translation symmetric ansatz in Eq. ( AI8O) the 
180° rotation symmetry requires that 



-i-y 



'^ ['^oopy'nxpx ) ^' 



m 



The time reversal transformation requires that 



i-ru. 



m 



(A200) 



(A201) 



Thus the SU2BnO and SU2B0n symmetric ansatz have a 
fornr 



„,0 



=0, if m — even 



(A202) 



The other two PSG's SU2B[00, nn] leads to vanishing Uij 
and are dropped. 

We see that only 4 of the 8 SU{2) PSG's leads to sym- 
metric SU{2) ansatz. Thus there are only 4 SU{2) sym- 
metric spin liquids at mea n-field level. 

For ansatz of type Eq. (A17£) the parity symmetries 
also require that 



.nix „m„ 



'^P^(m) — VxpxVxpy'^rn 
Up^{m) ^VxpyVxpx^rn 
Up^y{m) ='n'^yUm 



(A203) 



For ansatz of type Eq. (AI8C) the parity symmetries re- 
quire that 



'^P^(m) —VxpxVxpy'^rn 



„m^ „m„ 



Up^y(m) -(-) 



m^m-u ^m 



(A204) 



For each choice of t/'s Eq. (|A199|), Eq. (|A202|), 
Eq. (A203), and Eq. (A204) allow us to construct ansatz 



Uij 



for symmetric SU{2) spin liquids 



APPENDIX B: SYMMETRIC PERTURBATIONS 
AROUND SYMMETRIC SPIN LIQUIDS 



u-m = ivxpyVxpx)'^ gPygp^Umg p^ g p^ 



(A196) 



Construction of symmetric perturbations 



'^m ~ [VxpyVxpx) 






(A197) 



Let us consider a perturbation Suij around a general 
mean- field ansatz Uij . We would like to find all the sym- 
metric perturbations that do not change the symmetries 
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of the original ansatz. Let PSGq be the PSG of the 
ansatz Uij. Clearly, the PSG of the perturbed ansatz 



Uii 



5u. 



the two ansatz 



ij, PSGi, is a subgroup of PSGq. If we require 



i-tj 



and 



Hj 



'ij 



to have the same sym- 



metry, then the two PSGq and PSGi must satisfy 

PSGo/IGGo = PSGi/IGGi (Bl) 

where /GGo,i are the IGG of PSGq^i. We see that the 
low energy gauge group of the perturbed ansatz is equal 
or less than the low energy gauge group of the original 
ansatz. 

In the next a few subsections, we will use the following 
steps to find symmetric perturbations. We first choose 
IGGi to be IGGq or a subgroup of IGGq. Second, 
we find all the gauge incquivalent subgroup of PSGq: 
PSGi G PS Go, which has IGGq as its IGG and satisfies 
Eq. (Bl). Last, we find all the ansatz that are invariant 
under PSGi. 



Symmetric perturbations around the Z2-Iinear 
state Z2A003Z 



First let us apply the above approach to find all the 
symmetric ansatz near the the Z2-linear state Eq. (39). 
Here by symmetric ansatz we mean the ansatz with 
the translation, the time reversal and the three parity 
symmetries. The IGG that leaves the Z2-linear ansatz 
Eq. (||) invariant isQ = Z2. The PSG of ansatz Eq. (||) 
is given by 






Gy{i) =rO 

GpAi)=r" 
Grii) =(-)V° 



(B2) 



This is one of the PSG of the Z2 spin liquids labeled by 
Z2A0032;. We note that the IGG for any ansatz Uij is 
at least Z2. Thus the above Z2 PSG is already minimal, 
ie non of its subgroup can be regarded as the PSG of 
some symmetric ansatz. Therefore, all the symmetric 
perturbations around the Z2-li near state are invariant 
under the above PSG Eq. (|b|). From Eq. (|A64[ ) and 
Eq. ( A65| ), we find that the most general symmetric spin 
liquid with PSG Z2AtItIt^t^ have the form 



'^i,i+m 


=uLt-'\i-- 


=1,2,3 


1,2 
^P.yim) 


-~ul;i 




^P.yim) 


="L 




1,2,3 
""PAm) 


— "■m 




1,2,3 
^Pyim) 


=ul^2,3 





=0, for m = even 



(B3) 



3. Symmetric perturbations around the 

t/(l)-linear state UlCnOln 

The above analysis can also be used to obtain all the 
symmetric perturbations around the f/(l)-linear ansatz 

in Eq. (p2|). The invariant gauge group is Q = {e^^^ }. 

The ansatz is invariant under translations by x and y 
followed by gauge transformation ir*^ and zr^", where 
ir^ = i{cos9T^ + sm0T^). The ansatz also has the time 
reversal and the three parity symmetries. The PSG of 
the ansatz is generated by 



G,(i)=iA, 
GpAi)M'y^93{0p.) 



Gyii) ^ir'- 
Gpy{i) = {-y^g3{epy) 
Grii) =(-)'53(^t) (B4) 



where ga{0) = e"-"^ and 0x,y,px,py,pxv^ can take any 
values. We see that the ansatz Eq. ( p2| ) is labeled by 
UlCnOln. 

First let us consider the symmetric perturbations that 
do not break the U{1) gauge structure. Since the per- 
turbed ansatz are required to invariant under the same 
IGG and have the same symmetry as the original U{1)- 
linear ansatz Eq. (p2), the perturbations must be invari- 
ant under the original PSG Eq. (B4). The translation 
symmetry require the perturbations to have a form 



Su^l^r^ + i-ySut^r' 



(B5) 



The 180° rotation symmetry PxPy requires that 
and the time reversal symmetry requires 



Thus the symmetric ansatz with PSG Eq. (B4) are given 

by 



i„,3 ^3 



"^rn —^1 ^'^^ '"^ ~ even 



7/0 


=uL 


"Px„(m) 


= -'uL 


^P.(m) 


H-r^u'^j 


"■P„(m) 


={-r'<^ 



(B6) 



The above represent most general ansatz around the 
t/(l)-linear state that do not break any symmetries and 
do not change the quantum order in the state. 

Next we consider the symmetric perturbations that 
break the C/(l) gauge structure down to a Z2 gauge struc- 
ture. The IGG becomes Q = Z2 for the perturbed ansatz. 
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We first need to find subgroups of Eq. (B4) which have 
the reduced IGG and the same symmetries. The elements 
in new PSG have the following form 



Gp^ii) = ± (-)*-g3(V) GPyii) = ± i-y^gsiOpy) 

- '■' ' -■-'- GTW=±(-r53(^T) 

(B7) 






where 0x,y,px,py,pxy,T each takes a fixed value. The ± 
signs are independent from each other and come from 
the Z2 IGG. Different choices of 9x,y,px,py,pxy,T give us 
different subgroups which lead to different classes of Z2 
symmetric perturbations. 



To obtain the consistent choices of 9 



X ,y ,px ,py ^pxy 



T, we 



note that Gx.v, Gp ^,p„,p^ and Gt r nust s atisfy Eg. (A2), 
Eq. (I A^) Eq. (|A9|)7 E q . 1^10|) , Eq. ([aTTI) , Eq. @, and 
Eq. (A13|), with Q = {±r'^}. Those equations reduce to 



T ^T "T ^T 



= ±t" 



A53-i(0T)T«-53(0T)=±rO 



^g3\0pxy^g3{0px)^±T° 
yyg3{9px)^±r^ 



^"53 



^pxj 



(B8) 



(B9) 



(BIO) 



When Gx = ir^, and Gy — ir^, we fin d th e follo wing 8 
solutions for Eq. 



Eq. ( |B9|), Eg . (|Biq ), Eq. (|B11[ ) 



Eq. ( pT^ ), Eq. ( pT^ and Eq. JBU]), with G = {±7^ 



Gx{i)^iT'^, Gy{i)^iT^, 
{-y^GpAi)^{'y^Gp^{i)^T°, tr' 
Gp(i)=zr\iT^ i-yGT{i)^i 



(B18) 



We can make a gauge transformation Wi = {ir^Y (see 
Eq. ([4^)) to change the above to 

Gx{i)=T°, Gy{i)=T°, 

Gp^(i) = Gp^{i) = T\ t{-rr^ 
Gp^^d) =*t\ z(-)V2 GT(i) = (-)V°, ^T^■ 

Q ={±t"} (B19) 

We can use a gauge transformation Wi = (— )*^ to change 
Gp^^{i) — i{—yT^ to Gp^^{i) ~ iT^ without affecting 
other G's. Now we see that Gp^ (i) — {{—Yt'^ and 
Gp^^{i) = ir^ give rise to gauge equivalent PSG's. Thus 
we only have 4 different PSG's 

Gx{i)^T", Gy{i)^T", 

Gp^{i)=Gp^{i) = T\ ^{~Yr^ 



Gp^^{i)=iT^ GT(i) = {-yT\i 
Q ={±rO} 



(B20) 



'gz\(^pyy'9AQp 



±r^' 



A53"'(^P.)T'^53(^py)=±r° 



^Ox^Opxy^Oy^Op^y — -^ ^0 



93\0T)g3\epx)g3{0T)g3{Opx) = ± r° 



93\(^T)93\epy)gsieT)93{0py) 

53"'(^t) A-"53(^r)T^-'' = ± r' 



.0 



r''-''ff3(^p.)r^-''53~'(^py)=±r° 

93{6py)93{0px)g3^{0py)g3^{0px) = ± T° 

Since F^^ ^ p2 ^ p2^ ^ 7^2 ^ 1^ we also have 



ff3(M=±T° 



53 \"py 



)=±r° 



We can choose a gauge to make 6'^: = 0. Eq. 
two solutions 



Gx ^it\ 
Gx ^ir^, 



Gy =iT^ 
Gy =iT^ 



(Bll) 



(B12) 



When Gx = ir^, and Gy = ir"^, we find t he foll owin g 
4 solutions for Eqjp|), Eq. ( |B9|), E q. ( pTo| ), Eq. ( |Bll[ ), 
Eq. ( pT^ ), Eq. (pl3|), and Eq. (|B14D, with G = {±V^ 



G,(i) =iT\ Gy{i) = zt2, 



(B21) 



We can make a gauge transformation Wi = {ir^y (ir^)*« 
(see Eq. (ES)) to change the above to 



(B13) 



(B14) 



(B15) 
q) has 

(B16) 
(B17) 



i^3 



Gxii) =i-p'r", Gyii) 
GpSi) = Gp^{i)=T\ i{-) 
Gp.„(i)-(-r-^''*ri GT(i) = (-)VO, zr^; 

S ={±T°} (B22) 



We note the 4 PSG's in Eg . ( |B22| ) can be obtained from 
the 4 PSG's in Eq. (B20) through the transformation 
Eq. ( |AT13| ). 

We find that all the symmetric spin liquids around the 
t/(l)-linear state Eq. ( ^ that break the U{\) gauge 
structure to a Z2 gauge structure can be divided into 



eight classes. Using PSG's in Eq. (B20), we find transla- 



tion symmetry requires the ansatz to have a form 



- "^ -^ (B23) 
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The 180° rotation (generated by PxPy) symmetry re- 
quires that 



Symmetric perturbations around the 

5f/(2)-gapIess state SU2AnO 



Um =U- 



which imphes u^ = 0. The time reversal T symmetry 
requires that 



or 



r'^Um.T^ 



The four ansatz for PSG's in Eq. (B20) are given by 

Z2A[T«T0,r3T3]TiT0: 



_ I I 



(B24) 



l\^'^ —Q, for m = even 



a.iidZ2 A[tItI,tItI]t^tI 



UiA+m =ulnT^ + uI^t"^ 



Using PSG's in Eq. (B22), we find translation symme- 
try requires the ansatz to have a form 



"2,2+™, — I, ) ^rn' 

the 180° rotation symmetry requires that 



(B26) 



{-y^-^yum = (-) 



^ ('_^(4x+mx)m!,yt 



In this subsection, we would like to consider the sym- 
metric perturbations around the 5'C/(2)-gapless ansatz 
Eq. (|3^), which describes a SU2AnO spin liquid. The 
invariant gauge group is G = SU{2). The PSG of the 
ansatz is generated by 



Gx{i) =gx, 
GpAi)=(.-y^9p. 
Gp^y (i) =gpxy 



Gy{i) =gy 

Gpy{i) ={-)"' 9pv 
Grii) H-Ygr (B29) 



where gx,y,px,py,pxy,T G SU{2). Thus the 5t/(2)-gapless 
state is labeled by SU2At°t^. 

First let us consider the symmetric perturbations that 
do not break the SU{2) gauge structure. To have the 
SU{2) gauge structure, the perturbations must be in- 
variant under the gauge transformations in G and satisfy 



Su. 



XJ 



To have the symmetries, the perturbations 



(B25) must be invariant under PSG in Eq. (B29). The transla- 



tion symmetry require the perturbations to have a form 



5uiA+m = Su^T° 



(B30) 



The 180° rotation symmetry PxPy and the time reversal 
symmetry T require that 



-5ul, =Su%,{-y 







,1,2,3 
^m 



=0, for rrix — even or my = even 
^m" —0) f'^'' ^^x = odd and niy — odd 
The time reversal T symmetry requires that 

Um — V ) ^m 



T^UmT^ 



Z2B[TlTl,T^T^yT°: 



The four ansatz for PSG's in Eq. (B22) are given by 



'^i,i-\-m (, } 






(B27) 



,1,2,3 



andZ2B[T0r0,r3T3]TV3: 

Ui^i+mM~y^"'A^Lr'+ulT') (B28) 

"^m ^0: foi' "^2: = ^dd and niy = odd 

The eight different Z2 spin liquids have different quan- 
tum orders. They can transform into each other via the 
C/(l)-linear spin liquids, without any change of symme- 
tries. Those transitions are continuous transitions with- 
out broken symmetries. 



Thus the symmetric ansatz with PSG Eq. (B29) are given 
by 



U- -^ -U° T° 
^z,i+m — ^m 



"m — Oj f<3^ "^ — even 



7/° -7/" 



— ('_V"^7/0 

P^im) —\ I "m 

m 



(-)"^< 

\irrr.) =(-)"''"^ 



(B31) 



The above represent most general ansatz around the 
S'?7(2)-gapless state that do not break any symmetries 
and do not change the quantum order in the state. It de- 
scribes the most general S'C/(2)-gapless state with quan- 
tum order SU2Ar°r°. 

Next we consider the symmetric perturbations that 
break the SU{2) gauge structure down to a U{1) gauge 
structure. The invariant gauge group becomes G = U{1) 
for the perturbed ansatz. We first need to find sub- 
groups of Eq. (B29) by choosing a fixed value for each 
gx.y.px.py.pxy.T- We choosc gxy.px.py.pxy T in such a wa y 
that E£^ (|A2|), Eq. (g|), Eq. (g9|), Eq. ( |ATo| ), Eq. (^), 
Eq. ( |A5| ), and Eq. ( A13 ) can be satisfied when we limit G 
to a [/(I) subgroup. (Those equations are always satis- 
fied when G — SU{2)). Since the original invariant gauge 
group is formed by constant gauge transformations, its 
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f/(l) subgroup is also formed by constant gauge trans- 
formatfons. We can choose a gauge such that the C/(l) 
invariant gauge group is given hy G — {53(^)1^ € [0, Stt)}. 
To obtain the consistent choices of gx,y,px,py,pxy,T, we 
note that Gx,y, Gp^,p„,p^ and Gt must satisfy Eg (A2), 
Eq. (|A^) Eq. ©/Eq. JKl^ , Eq. (|aTT1) , Eq. @, and 
Eq. |\13D , with g = {g3{0)\e e [0, 2tt)}. Those equations 
reduce to 



J=^9y9x ^9y ' e U{1) 



(B32) 



5.~'5r 5x5T e[/(l) 5;'5t 5y5T ef/(l) (B33) 



dxgpxQxgpx et/(i) 5y ^9px9y9px et/(i) (B34) 



9y9pv9y9pv e^(l) 5x ^9py9xgpy &U{1) (B35) 



5y ^9pxy9x9pxy et/(l) 5^ ^gpxy9y9pxy &U {\) (B36) 



5t 9px9T9px eC/(l) ffj, 9py9T9py et/(l) 

9T^9pxy9T9pxy ^U (1) 



(B37) 



9pxy9px9pxy9py <^U{1) gpygpxQpyOpx ^U {I) (B38) 

Since P^^ ^ p2 ^ p2^ = 72 = 1, we also have 



2 
ypxy 



gly eL/(i), 



sLy ec/(i), 



apt, 
^2 



g^ eL/(l). (B39) 



Solving the above equations, we find 16 different PSG's 
with C/(l) invariant gauge group. The following are their 
generators and their labels. 

UIAt^t? [r°, r^][T" , rl] ( which is gauge equivalent to 
UlATOr|[r",Ti][rO,r|]): 



Gp.y{i)=93{0pxy), iA^" 
UlArirl [rO,ri][TO,rl]: 

Gp.„(^)=33(6'p:.y), jA^" 

UlCT£Tf^[ri^,ri][r!l,ri]: 

(-)^-Gp^(i)=g3(^p.) 
Gp.„(i)=ff3(^p..y), iA^" 



G,W=53(e,) (B40) 

-)■'=' Gp„(i)=.g3(ep,) 



Gy{i)=g3{ey) (B41) 

-)■'«Gp„(^)=^A« 
-)^GT(i)-,g3(^T), it'^ 



Gy(i)^iT'^y (B42) 

-)■'=' Gp,(i)=.93(0p,) 
-)^GT(i)=ff3(^T), it'^ 



UlCrlri[r!^,r|][T'^,rl]: 

G,(2)=iA, Gj,(i)=zA (B43) 

{-^Gp^ii) - zr"- (-)■'" Gp,(i) = *A» 

Gp.Ji) = g3(Vy)> »^'''^" (-)'GT(i) = g3(^T), */^ 

The 180° rotation symmetry PxPy and the time rever- 
sal symmetry T require that 



^ — i^—i — m ^2,i+mv / ^- 

— Ui^+m —grUi^i+rngx (~)" 



We find that all the symmetric spin liquids around the 
5L/(2)-gapless state Eq. @ that break the SU{2) gauge 
structure to a U{\) gauge structure can be divided into 12 
classes. They are given by (using abbreviated notation) 
UlCnO[0,lln: 



3 _3 



= ul,T^' + {-yu^^r 



u^ —Q, for m — even 



(B44) 



G,(^)=i/^ Gy{i)=iT''y, 
{-Y'-Gp^d) =g3(V), {-r^Gp^ii) = gsiOpy), 

Gp^M =g3i0pxy), it'--- i-yGrii) - gsiOry, 



UlCnO[0,ll 



Ui,i+m = ul,T° + {-yul^T^ 

^771 ~0i foi' "^ = even 
n't, —0, for m ~ odd 



(B45) 



G,(^)=^A, GJ,(^)=^/^ 
(-)*^GpJi) =53(V), i-r^Gp^ii) = gal^p,), 
Gp.Ji) =53(0p.y), *A- (-)'GT(i) = »r^-; 

UlCxl[0, l]n: 



3 _3 



W^,i+rr^ = W^l^r" + (-)*U^r' 



,0,3 _ 



0, for m = 



(B46) 



G,(^)=^T«^ G,(^)=^r»^ 
(-)*^Gp,(i) =*tV, (-)'»Gp,(i) = *r^™, 

Gp.„(i) =53(^p.y), »A- (-)'GT(i) - gsiOry 

UlCxl[0, l]x: 



,0 



= ^$„r" + (-)^uLr3 



,3 



w^ =0, for m = odd 

Gx{i)^iT''% Gy{i)=tT'y, 
i~y- Gp^ (i) =*r«- , i-yy Gp^ (i) = * A 

Gp.y {i) ^gziepxy), ^A- {-yGrii) = it' 



(B47) 



' UT . 
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UlAri,0[0, l]a- (which is gauge equivalent to U1A00[0, 1]1) 



Ui.i+rr,. = u™r° + uLr^ 



(B48) 



^m =0' fo^' ''^ ~ even 

"m =0' foi' "^ = odd 

G.(i) -33(0.), G,(i)=.g3(0,), 

(-)^^Gp,(i) =.93(V), H'^GPyii) = 93{0py), 

Gp.Ai) =g3iep.y), *tV« i-yGrii) = zr"-; 
UlAxl[0, l]a;: 

?/ ■ ■ I = 7/'^ T*^ + 7/'^ T'^ 

M^ =0, for m = odd (B49) 

G,{i) ^g^ie,), Gyii) ^ gsiOy), 

(-)'"= Gp,. (i) ^ir'^'"' , (-)'" Gp„ (i) = ir«™ , 

Gp.,(i) =g3{ep.y), *A- (-)'GT(i) = zr^^; 

The 12 different C/(l) spin hquids have different quantum 
orders. 

We can use a gauge trans form ation Wi = {ir^Y to 
make ansatz Eq. ( B44 ) to Eq. ( B47 ) translation invariant. 
Weget UlC[7iO,a;l][0,l]?i: 



'^m ~^- fo'" '"^ ~ even 
G,(i) =tO, G^(i) = r°, 
GpAi)^Gp^{i)^T\ ZT^ 
Gp^^{i)=T\iT' GT{i)^i-rr"; 

andUlC[nO,a;l][0,l]x: 

uL =0, for m = even 



G,(z) =tO, Gy{i) = r°, 
Gp.(i)=Gp„(i)=T", zri 

Gp(i)^T°,lT' GT{i) = 



(B50) 



(B51) 



^r^; 



Now we consider the symmetric perturbations that 
break the SU{2) gauge structure down to a Z2 gauge 
structure. The invariant gauge group becomes Q ~ Z2 
for the perturbed ansatz. We ch oose a fixed value for 
each .ar.i;.nx.m;.nxii^ in Eg. (|B29|) such that Eq. (A2), 



iX.V .VX.W.VX' 



Eq. (^), Eq. (^), Eq. (IATOD; Eq. (|All|), Eq. (^), 
and Eq. ( Al3| ) can be satisfied when we limit Q to a, Z2 
subgroup. Since the original invariant gauge group is 
formed by constant gauge transformations, its Z2 sub- 
group is also formed by constant gauge transformations, 
which is given by ^ = {ztr*^}. 

To obtain the consistent choices of gx,y,px,py,pxy,T, we 
note that G^^y, Gp^^p^^p^^ and Gt must satisfy Eq. (|A^), 



Eq. 
Eq. 



( J AJ), Eq. (|A9|), Eq. ([aTo|), Eq. ([aTiI), Eq. (^), and 



(A13), with Q = {±t"}. Those equations reduce to 



JxQyg^ ^gy ^ = ±T° 



(B52) 



gx^gT^9x9T = ± r" gy^gr^gygr = ± r° (B53) 

9xgpxgx9px = ± T° gy^gpx9y9px = ± t° (B54) 

9ygpygy9py = ± r° g'^g'ygxgpy = ± r" (B55) 

9y^9pxygxgpxy = ± T° g^^ g^^ygygpxy = ± r° (B56) 



5t ^gpxgrgpx = ± r" 5j, ^gpygrgpy = ± t° 



'1„-1 



5t gpxygrgpxy = ± r 



(B57) 



gpxygpxgpxygpy = ± r" gpygpxgpygp^ = ± r" (B58) 

Since P^ ^ p2 ^ p2 ^ 7^2 ^ ;^^ ^g g^jg^ j^j^^g 



'2 =±rO g 



2 -±r°. 



«2 - -I- T-0 „2 _ , 

gpxy — ±T grp — zLT . 



(B59) 



From Eq. ( |B58| ), we see that gpxPpy = ±gpygpx- 9px,py 
has the following 5 gauge inequivalent choices 



gpx =T"°, ir^ 
Qpx — ^'^ 



9py — ^ 7 ^^ 

5py =2r 



(B60) 
(B61) 



Also, according to Eq. ( |B58| ), gp. = ±gpxygpygpxy This 
requires gp. = gpy if [gpx,gpy\ = 0. Similarly, we also 
have Qx = gy if [gxTgy] = 0. Many solutions can be 
obtained by simply choosing each of gx,y,T and gpx,py,pxy 
to have one of the four values: (t°,?t-'^'2'^). We find the 
following 65 solutions 



Gx —Gy — T , 



i-Y^Gp^ ={-Y'Gp^ = r", ZT^ Gp^^ =t\ zr^; 

(B62) 



Gx — Gy — T , 



-y - ' , (-)*Gt =r°, 

{-Y'-Gp^ H-Y'Gp^ = zri Gp,.„ =zt2; (B63) 



Gx — G,, — T , 



(-)^^Gp, =7ri 
'P 



Gp.., =zri2; 



Gx —Gj, — T 



{-Y^Gp^ ^IT^ 



{-YGt -^T^ 



(B64) 



^Y-Gp^ H-Y^Gp^ = T°, zt3 Gp,„ -rO, zr^; 



(B65) 
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(^x —try — ''' J 



{-y''Gp^={-r'Gp^=T\ I 



(B66) 



„-^i.2,3 r^ i-A 



Gx ^it\ 



Gy —IT , 

i-yGT^T°,zr^ 
i-Y^Gp^ =^r2; (B77) 



Gx ^Gy = T°, {-YGt =^r^ 

(B67) 



Gx ^iT\ 



Gy ^It\ 

{-YGt =t°, ir^ 



i-y^Gp^ =r"; (B78) 



Gx ~Gy — T , 

i-y^Gp^ =^r^ 



{-YGt ^t\ 
{-y«Gp^ ^ir' 



G 



p =IT 



Using a gauge transformation Wj — (iT^)*(±)'==, wc 
can change Eq. ( B69 ) - Eq. ( B75 ) to standard Z2 form 
Eq. (^) or Eq. (|62D. We choose the ± to remove the 
(B68) (-)* factor in Gp^^. We get 



Gx Gy — IT" 



(-)'^Gp, H-y-Gp^ = r", tr' Gp^^ =r°, zr^; 



(B69) 



Gx — Ct^ — IT 



i-y-Gp^ H-y^Gp^ ^ r0,zri'2,3 Gp.„ ^zr^; 






(B70) 



Gj; — Gy — IT , 



(-)'Gt =r«, 



(-r^Gp, =(-^Gp^=*^' Gp,„=t", 



„-^l,2,3. 



(B71) 



Gx =Gy = rO, {-YGt =rO, 



Gp =Gp =r°, zt3 



Gp =r°, ir^; 



(B79) 



Gx =Gy = T , {—y^Gr =T , 

Gp. ^Gp^ = T°'\ i-yiT^'^ Gp^^ ^it'; (B80) 

Gx =Gy = T , {—YGt —t , 

Gp. =Gp^ ^ i-Yzr' Gp^^ =r°, zr^^^'^. 



G, -G, - r°, (-)^Gt =ir3, 

Gp. =Gp„ = T°, zT^ Gp,„ =T°, ir^; (B82) 



Ga; — Gj, — IT 



'Gt ^it\ 



Gx —Gy — T , 



{~yGT =it\ 



i-y-'Gp^ ={-y^Gp^ = t\ ir^ Gp^^ =rO, zr^; Gp^ =Gp^ = t"'^ (-)^t1-2 Gp^^ =it^- (B83) 



(B72) 



Gx ^Gy = zr3, {-YGt =zr3, 

(-)''"'Gp, H-y^Gp^ = r°, zt1^2,3 g^^^ ^^^1. 



Ga; — Gj, — T , 

Gp. =Gp^ = (-)^ri Gp^^ =r' 



(-)'Gt =zt3, 

•°, zri'2^3. (B84) 



(B73) 



Gt =ir\ 



Gx =Gy = zt3, (-)^Gt -zr3, 



^0 ,vl>2,3. 



(B85) 



(-)^^Gp, ={~y^Gp^ = zri Gp,„ =r", zr 



.0 ,vl.2,3. 



(B74) 



G, =Gy = it\ 



{-YGt ^IT^ 



'X — "-fy 

{-y^Gp^ ={-yyGp^^T\ zt1'2.3 Gp^^ =t°,zt1'2.3. 



Ga; — Gy — T , 

/^ /^ „0.3 ,' \i„- 1.2 /-< _0 „■ 1 

Gp. =Gp^ = r , (-) zr Gp^^ =t , zt 

Using a gaug e tra nsform ation Wi = (zT^)*="(zr^)*'', we 
can change Eq. (B7£ ) - Eq. (B78) to more standard forms 
(seeEq. (|l|) and Eq. ^): 

G, =(-r«rO, ^ — ° 



Gp.„ ={~y^'y 



^ *" *'' IT 



Gy=r' 

Gt =(-)V0, zr^ 



(B75) (-)'^Gp^ =ir2, (-)^''Gp„ =zt1; (B86) 



Gx — ZT , Gj, — ZT , 

Gp.„=Zt12, (_)iGj,=T",ZT3, 

(-)*^Gp, =zt2, i-y-Gp^ =zt1; (B76) 



G. =(-)*''t", Gy =tO, 

Gp.„ =(-)^^'nTi2, Gt =(-)V0, zt3, 

(-)'''Gp. =zt\ (-)'=^Gp, =zt2; (B87) 
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G. -(-)*«tO, 


Gy =t\ 


Gp.„ H-r-''^T^^ 


Gt =(-)V0, iT^ 


Gp^ -r°, 





(B88) 



Now we can list all PSG's that describe the Z2 spin 
liquids in the neighborhood of the S'?7(2)-gapless state 
using the notation Z2A... or Z2B.... We find that the 
65 PSG's obtained before lead to 58 gauge inequivalent 
PSG's. 

In the following, we will list aU the 58 PSG's. We 
will also construct ansatz for those PSG's. First let us 
consider PSG of form Z2A.... For those PSG the ansatz 
can be written as 



^x,z-\-m 



(B89) 



the 180° rotation and the time reversal symmetries re- 
quire that 



^Um{-y 



Wcfind 



^i.i+m ^^m^ ' ^m'^ ' ^m^ 



,0,1.2 



=0, for m = even 



For PSG 



Z2^rlriT0'i'2.3^3 



(B94) 



(B95) 



In the following we consider the constraint imposed by 
the 180° rotation symmetry and the time reversal sym- 
metry. 
For PSG 



Z2At^t1t\° 



(B90) 



(here we have used the notation [a, b] [c, d] to represent 
four combinations ac, ad, be, bd) , the 180° rotation sym- 
metry generated by PxPy requires that 

u =u (-)'^ = u'' 

^ — m ^m\ } m 

The time reversal symmetry T requires that 

^m '^'m\ ) 

The above two equations require that Uij ex t'^ , which 
describe SU{2) spin liquids. 
For PSG 



the 180° rotation and the time reversal symmetries re- 
quire that 



=Um(" 



'i-y 



The ansatz has a form 



Ui.i+m — Wm, 


r" + ul^T^ + ul^r 


uL-0, 


for m ~ even 


u]^ =0, 


for m = odd 


For PSG's 




Z2AtIt^_ 


[t\1 



Z2At1t1tO'1'2.3^3 
Z2A[tItI,tItI]t^^^-\1 



(B96) 



(B97) 



the 180° rotation and the time reversal symmetries re- 
quire that 



Z2At3t3t0^1'2.3^o 
Z2A[r°T°,T^T^]T°'i'V! (B91) 

the 180° rotation symmetry and the time reversal sym- 
metry require that 



U—Yn ^m ^ 



The above two equations give us 



m 
m 



For PSG 



^In'^ =0, for m = even 



Z2At^tIt^^t\ 



(B92) 



(B93) 



-m ^t^m 



^m, ^ ^m,^ 



m 



The ansatz has a form 



u^ =0, for m = even 
ul^ =0, for m = odd 



(B98) 



For PSG Z2Arl r^T^V^, the 180° rotation and the 



time reversal symmetries require that 



m _ „t 



Wc find Uij ex T^ . The spin liquid constructed from u 
is a SU(2) spin liquid. 



^3 
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For PSG's 



For PSG 



Z2AtItIt°^^^^-^tI 



(B99) 



the 180° rotation and the time reversal symmetries re- 
quire that 

The ansatz has a form 

(BlOO) 
There are six PSG's of form Z2B.... The first two are 



11 12 



Z2B[TWl,TlTiy\l. (BlOl) 

The 180° rotation symmetry requires that 



0,1,2 _ 



"^m — Oj ^'^^ ''^x = odd or ray — odd 
The time reversal symmetry requires that 

We find 

u^^'2 =0, for m = even (B102) 

For PSG 

Z2BtItIt^t°, (B103) 

the 180° rotation symmetry requires that 



(-)'^™''^m = (-) 



^ f-_^(«x+mx)mj,yt 



^TTi ~0, for nix — even or m^ = even 



1,2.3 



utL~'" —0, for rria; = odd and niy = odd 



The time reversal symmetry requires that 
The above two equations give us 



uif '^ =0, for m — even 



(B104) 



Z2B[TlTl,TlTiy\l 



(B105) 



the 180° rotation symmetry requires that 



or 



"^m ' —0, for nix — even and niy = even 
u^ —0, for nix — odd or niy — odd 

The time reversal symmetry requires that 



3 3 



UiA+m =(-)'^'"« (uI^t' + ul,T^) (B106) 

u]^ =0, for nix = even and niy = even 
For PSG 



Z2BtItIt\1 



the 180° rotation symmetry requires that 



(B107) 



(_)».™.y^ = (-) 



_ ^_^(ix+nix)myyt 



or 



u^ —0, for mx = even or niy = euen 

,1,2,3 



=0, for 



= odd and 



= odd 



The time reversal symmetry requires that 



^^m '^ ^m^ 



The ansatz has a form 

Ui,i+m =(-)'-"" (Mj„ri + ul^T^) (B108) 

Uj^ =0, for nix — odd and ?7ij, — odd 

We find that there are 52 different Z2 spin liquids in 
the neighborhood of the S'[/(2)-gapless state Eq. (pCJ). 
Those Z2 spin liquids can be constructed through Uij . 



5. Symmetric perturbations around the 

S'f/(2)-linear state SU2BnO 

In this subsection, we would like to consider the sym- 
mctr ic p erturbations around the SU {2) -linear ansatz 
Eq. (|3l|), which describes a SU2BnO spin liquid. The 
invariant gauge group is G — SU{2). The PSG of the 
ansatz is generated by 



Gx{i) H-y^g., 

GpAi)M~y^9px 
Gp.y{i)H-y^''9pxy 



Gy{i) ^gy 

Gpy{i)M~y'9py 
GtW H-ygr (B109) 
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where gx,y,px,py,pxy,T G SU{2). 

First let us consider the symmetric perturbations that 
do not break the SU{2) gauge structure. To have the 
SU{2) gauge structure and the symmetries, the pertur- 
bations must be invariant under PSG in Eq. ( B109| ). The 
translation symmetry require the perturbations to have 
a form 



:.t-\-m 



{-y-'^-sui^ 



(BllO) 



The 180° rotation symmetry PxPy and the time reversal 
symmetry T require that 

.0 jr„.0 / \m 



-5u^l^ ^su^-y 



Thus the symmetric ansatz with PSG Eq. (B109) are 
given by 

Ui,i+rr^ -(-)*^"'' li^r° (Bill) 

"^m =0) foi' "^ = even 

The above represent most general ansatz around the 
5'f/(2)-linear state that do not break any symmetries and 
do not change the quantum order in the state. 

To obtain other symmetric perturbations around the 
5?7( 2)-line ar state SU2BnO, we can use the mapping 
Eq. (|A112|) and Eq. (|A113| ). We first note that the PSG's 
that describe the spin liquids around the 5t/(2)-linear 
state can be obtained from those around the SU{2y 
gapless state SU2AnO. This is because mappin g de - 
scribed by Eq. ( |A112| ) ma ps the SU2AnO PSG Eq. ( P29| ) 
to the SU2BnO PSG Eq. ( |B109| ). Using this results, we 
can obtain all the PSG's that describe the spin liquids in 
the neighborhood of the 5'C/(2)-hnear state. The SU2AnO 
PSG for the 5f7(2)-gaplcss state have 16 dif feren t sub- 
groups with IGG=C/(1) (see Eq. (picl) - E q. (p43|)) an d 
58 subgroups with IGG=Z2 (see Eq. (|b62D - Eq. (|b88D). 
Therefore, the SU2BnO PSG of the 5[/(2)-linear state 
also have 16 different subgroups wi th IG G—Uj l) an d 58 
subgroups with IGG=Z2 (sec Eq. ( ^62| ) - Eq. (|B8^)). 

The 16 s ubgr oups with IGG=C/(1) can be obtained 
from Eq. (B4C) - Eq. ( p34S ) through the mapping 
Eq. ( [ATi^ ). They arc UIBT^tO [r^, ri][rO , rl] ( which 
is gauge equivalent to UlBT^^r^fr^, r-'^][T'' , r]:]): 

Gx{i) = i-y-g3{0x), Gy{i) = g^{dy) (B112) 

{-y-Gp^ii) = gsiOpx) i-yGrii) = gsiOr), ir'^ 
hy^Gp^ii) = g,{9py) (-)*-*" Gp,„(i) - g3{9pxy), *A- 

UlBriT|[rO,ri][T£,Ti]: 

Gx{i) = i-y-g3{9x), Gy{i) = g3{0y) (B113) 

hy^Gp^ii) = *T«- i-yGrii) = asiOr), it'^ 
i^y^Gp^ii) = zA" i-y^-'^Gp^^ii) = 53(Vy), *^'^^" 



UlCr^ri^TO [tO,t|], UlCrOrVirl and UlCrOTOrV°: 

Gxii) = z(-)'» A, Gyii) = tr'y (B114) 

i~y-GpAi) = gsiOpx) i-yGrii) - 53(^t), ^r'- 
i~y^Gp^ii) = gsi0py) i-)'-'yGp^^ii) = gsi0pxyy »A- 

UlCTJr|[rO,rl][rO,Ti]: 

Gxii) = i(-)^''T«^ Gyii) = iT^'y (B115) 

(-)'==GpJi) = zr"- i-)'Grii) ^ gsier), ^r'^ 
i-)''Gp^ii) = zA» (-)^^^«Gp,^(i) = g,i9pxy), *A- 

The labels for the last two equations are obtained by 
making a gauge transformation to put Gx,v in a more 
standard form (see Eq. (pT2^) and Eq. (pl2^). 



After obtaining the PSG's, we can construct the ansatz 
which are invariant under those PSG's. We note that for 
the above PSG's the time reversal symmetry T requires 
that 



— Ui^i+m —gTUi.i+rngr (~)"* 

gr ==T°, iT^ 
and the 180° rotation symmetry PxPy requires that 



-I, — I— m f^i,i+m 



(-)"^=uT. 

V / —x—m,—x 



(B116) 



When Gxii) = i(— )*»t^^, Gyii) = ir^^ , Uij has a f orm 
Ui^+rr^ = (-)'-"« (uj^r" -f i-ful^T^). Eq. (^11^ ) rc- 
duces to 

^rn =0; ^^ "^a: = even and nriy — even 
^m —0' if ^^x — odd and niy — odd 

When Gxii) = (-)'" 53(^2;), Gyii) = gaidy), Uji h as a 
form Ui,i+m = (-)'=="" (u^tO + u^t3). Eq. ( [Bll6| ) re- 
duces to 

"rn =0) if "^2: = even and my — even 
u^ =0, if nix = odd or niy = odd 

Using the above results, we find that all the symmetric 
spin liquids around the 5'C/(2)-linear state Eq. (^Tl) that 
break the S'C/(2) gauge structure to a C/(l) gauge struc- 
ture can be divided into 12 classes. They are given by 
UlCnO[n, l]n: 

Ui.i+rr. = i-)'^"'H^^Lr" + i-yul,T^) 

(B117) 

u^^ =0, if m = even 

Gxii)^ii-y-T'-, Gy^)^iT'y, 

i-y^'GpSi) =53(V), i-y^Gp^ii) = g^ie„), 

i-y^^'^Gp^^d) ^g^iOpxyy *A- i-yGrii) = 53(^t); 
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UlCnO[n,x]l: 



i-y^"'^iul,r' + i-rulr') 



(B118) 



u^ —0, if m = even 



^m — Oj if ^^x = odd or iriy = odd 

i-y^Gp^i) =53(V)- (-)^''Gp„(i) = 53(^p,), 
(-)"■'« Gp^Ji) =53(Vy), *A- (-)'GT(i) = it'-; 

VlCn[n,x]n: 

(B119) 
""tA =0; if m = even 



UlCll[n,x]x: 



(-)'^"''("^^° + (-)'"Lr^) 



(B120) 



""rn, ^Oj if m = even 



u'L. —0, if nij: = odd or TTiy = odd 






G,(i) =z(-)'«r«% G,(i) = zr««, 
(-)^-GpJz) =iT«-, (-)'«Gp„(i) = ir"-, 
(-)'"'" Gp,„(i) =53(Vy), *^'''^'' (-)'GT(i) - *T^-; 

U1B00[0,1]1: 

Ui^i+rr. = (-)'^"« (zi^r° + ul^r') (B121) 

u^ =0, if m = even 
^m =0j if '^a: = ^dd Or ?7iy — odd 
G,(i) =(-)^«.g3(0.), G,(i) = .g3(0,), 
(-)^^GpJi) =53(V), (-)^''Gp„(i) = 53(0p,), 
(-)'^^«Gp,„(i) =53(Vy), *^'''^" (-)'GT(i) - *r''-; 

UlBa;l[0, l]a;: 

Ui^i+rr. = (-)^=='"'' (^CrO + ul,T^) (B122) 

u^ =0, if m = even 
'^m — Oj if "^K = odd or ttIj, = odd 
G,{i)H~y'93i0.), Gy{i)^gsieyy 
(-)'^GpJi) =ZT^-, (-)'«Gp„(i) = zr"-, 
(")'^^«Gp,„« ^93iOp.y), it'--- i-yGT{i) = zr" 



— ,V"T . 



can change Eq. ( |B114| ) and Eq. ( |B115D to 



G, =iA, Gy ^ir'y, 
Gp. -(-)^==53(ep.), Gp„ = (-)^«.g3(ep,), 
GP., =(-)'==^"53(((-)''' - i-y-)l + 0pry), 

Gt =(-)'53(eT), (-)'^53(^T)^rl (3^23) 

Gp. =(-)^r^-, Gp^^iA" 
GP., =(-)■'==*« g3(((-)^'' - H'-)l+Or,.yy 

i-y-'^93i{{-y^ + i-y'')j+ep.yyr\ 



Gt ={-y93iOT), {-y-gzieryr^ 
We note that 



{-y-'^93{{{-y^-[-y-)-^+0p.y) 



(B124) 



and 



<-y-93{{-yj+e'^,y) 



(-r^^ssia-r^ + i-D^+v.) 



(B125) 



'-93{{-yj + e',,y) 



(B126) 



Thus the above two sets of PSG's are labeled by 
UlCrO tO [tO , tI][t^_ , tI] and VICtItI [t^ , ri] [t^ , rl] 
respectively (see Eq. I JMj ) - Eq. (fz^)). We also note 
that the PSG UlCr^r^rlr" is gauge equivalent to 
UlCT°T°r^r° . Using the gauge transformation Wi = 



.g3(-(-)V/8) we can change Eq. (§123|) and Eq. ( |B124| ) 
to 



Gj; = ir''" 



Using the gauge transformation Wi — g3((~)*"7i'/4) we 



G, =zt" 
Gp. -(-)'-g3(V) Gp„ = (-)^''g3(ep,) 

Gp.„ =i-y-'^93a{~y^ - (-)'-) J + Vy)>53(^p.,)«r\ 

Gt =(-)'<?3(0t), i-y«g3{{-y7r/A)^T''- (B127) 



G, =ZT«% Gj,=iA 
Gp,. =*.g3(V + i-y^y Gp^ = igs{9,y - H'^y 
Gp., =(-)^-'''g3(((-)^« - i-y-)^ + Op:..y),g3{0p.yW, 
Gt H-y93{0Ty {-y^93{{-y^/^)ir'- (B128) 

Then we can use the gauge transformation Wi = (jt^)* 
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to change Eq. (|B127|) and Eq. (|B128D to 



Gp., ^(-y^^ygsiii-y^ - {-yy)^ + {-y9p,y), 

93{{-y0p^y)iT\ 
gt ={-y93{{-yoT), {-y^gzii-yeTyr^'' (B129) 



"m =0: fo^' "^ ~ even (B133) 

Gx ^g3iH%), Gy = gsiH'ey) 
Gp. =(-)■'^93((-)'ep.)^r^ 

Gp, ^{-yyg3{H'e,yyr' 

Gp., ={-y^'^93{{hy^ - {-y^)^ + {-ye,.y), 



93{{-)%.vyr^\ 

Gt ={-y93{{~yeT) 



G.=g3{i-yo,), Gy ^ gsii-yey), 

Gp. ={~y-93{H%.yT'' Gp^ = i-y^g3iH%yyr'^ 
Gp., ^{-y^^yg,i{{-y^ i-y^)^ + (-)^v.)> 
93ii-yOp^yyT\ 

Gt ={-y93{{-yeT), {-y^93{{-yOTyr'' (BISO) 

We can use a gauge transformation Wi = 



[IT 



i\i. 



'g3(- (-)^7''/8). q3((-) '"7r/4) to simplify the ansatz 
Eq. ( |B117| ) - Eq. (|B12C| ). After the gauge transfor- 



mation, the IGG is given by {53((— )*0)}- The ansatz 
has a form Ui^i^rn = '"m'''^ + '^m'''^ fo^' '^ ~ '''^'^ ^^"i 
Ui,i+m = ""Vr" + ^^'''^ for m = even. We find the 



m 



ansatz Eq. (B117) - Eq. (B12C) become 



Um —0, for m = even (B131) 

G. =53((-r^.), G, = .g3((-r0,) 
Gp. =.93((-r^p.) Gp^ ^ gsii-yepy) 

Gp., M-y^'^93{{{-y^ {-y^)l + (-)%.y), 

93{[~yepxyyT^, 
Gt H-y93{{-yoT) 



'^m T^Oi foi' ''^x — even and my = odd 
'^m 7^0' for TTia; = odd and TOj^ = ewe?! 
'^m 7^0' fo^' "^a: = even and tti^ = ewe?! (B132) 
G. =.g3((-re.), Gy^9:i{{-y9y) 
Gp. ^93{{-)%x) Gp^ = 93{{-yOpy) 

Gp., M-y^'^93{{hy^ - {-y^)l + H%xy), 
93i{-yop^yyT^^ 

Gt =(-)'" 53((-)'0T)*ri 



UiA+m =u\nT^ + U^T^ + U^r^ 

w^ 7^0, for TOj; = even and niy — odd 

"^m 7^0, for nix — odd and niy — even 

""m, 7^0: fo^' "^2; = even and niy — even (B134) 

Gx^g3{{-)%), Gy^gsii-yey) 

Gp^ H-y^93{{-y0pxyr\ 

GP, H-y^93{i-)%yyr' 

Gp., H-y'''"93{{{-y^ {-y^)l + {-yo,xy), 

93{{-)%xyyT^\ 

Gt H~y'93{{~yoTyr^ 



In Eq. O132D , Eq. ( |B133| ) and Eq. ( |B134| ) we have made 



additional gauge transformation (r 



12 ^12 



)~>(7 



'■)■ 



Using the mapping Eq. ( A112 ), we can obtain all the 
PSG's for the Z2 symmetric spin liquids near the S'C/(2)- 
linear state from the 58 Z2 PSG's obtained in the last 
subsection for the 5't/(2) -gaple ss state. We note that, 
under the mapping Eq. ( A112 ), a Z2 PSG labeled by 
Z2Aabcd will be mapped into a PSG labeled by Z2Babcd 
and a Z2 PSG labeled by Z2Ba6cd will be mapped into 
a PSG labeled by Z2Aabcd. In the following, we will list 
all the 58 Z2 PSG's for the spin liquids near the SU{2) 
linear spin liquid. We will also construct ansatz for those 
PSG's. First let us consider PSG's of form Z2B.... For 
those PSG's the ansatz can be written as 






(B135) 



In the following we consider the constraint imposed by 
the 180° rotation symmetry and the time reversal sym- 
metry. 
For PSG 

Z2B[t1tI,t^tI][t°,t^]t1 
Z2Bt1tIt^t°_ (B136) 

the 180° rotation symmetry generated by PxPy requires 
that 
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For PSG 



^m ' =0, for nix = odd or niy = odd 
The time reversal symmetry T requires that 

The above two equations give us 



„,0 



For PSG 



0, for m ~ even 

odd or m^ 



(B137) 



odd 






(B138) 



the 180° rotation symmetry requires that 

or 

u^ =0, for nix = even or niy — even 

The time reversal symmetry requires that 

^m ^m V ) 

The above two equations give us 



For PSG 



the 180° rotation symmetry requires that 



(B139) 
(B140) 



,04,2 



=0, for nix = even and rUy = even 



The time reversal symmetry requires that 

We find 

u^i.2 ^Q^ f^j, ^ ^ g^g^ (B141) 



^0^0 _3_3 1^0,1.3^3 
Z2i?TiT>°'1^2,3^3 



(B142) 



the 180° rotation symmetry requires that 



yxnij,„_ ^ ^_^»Ti('_^(ix+nix)mj,yt 



(-)*-"«u^ = i-n-y 



or 



^rn =0i f'^'' '^a; = ewen and TOj, = even 
u^ ' =0, for nix = odd or niy = odd 

The time reversal symmetry requires that 

-Um ^T^UmT^i-)"" 

The ansatz has a form 

Ui,i+rr^ =(-)*^"« (u^t" + n^T^ + U^T^) (B143) 

M^ =0, for m = even 



=0, for nix — odd or niy — odd 



For PSG's 



Z2BtIt^t^tI 
Z2BrlTlr"^i'2,3^3 

Z2B[TV+,r^T3]T"^i'V3 (B144) 



the 180° rotation symmetry requires that 



^'"a7/.__ — ('_V«x+mx)myyt 



(_)*.™«^^ - (-)( 



or 



u^ =0, for nix = ewen or my = even 
u^ ' =0, for nix = odd and niy — odd 

The time reversal symmetry requires that 

The ansatz has a form 

Ui,i+m =(-)*^"« (uj^ri + ul,r^ + ul,T^) (B145) 

■"m, =0i f'^'' '^ = even 

u^ =0, for 777,1; = odd or ttIj, = odd 

For PSG 

Z2BTlr^riVi, (B146) 

the 180° rotation symmetry requires that 



or 



nix = even and 777,^ — even 
ul„ —0, for nix = orfd or rriy — odd 



M^i'2^0, for 



-'m 
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The time reversal symmetry requires that 
We get 

Ui,i+rr. =(-)^^"« {utr" + uI,t' + li^T^) (B147) 



The first two are 



'^rn =0; f'^'' '^ = even 

^m ^Oi for TTix — even or niy = even 

For PSG's 



Z2BtItIt°^^^^'^tI 
the 180° rotation symmetry requires that 



(B148) 



(-)'^™''^m = (-) 



^ ^_^(«x+m^)mj,yt 



u^ =0, for nix — even or TOj, ~ 
u^ ■ —0, for nix = odd and m^ = 

The time reversal symmetry requires that 



odd 



Um — T Urn^ 



The ansatz has a form 



Ui,i+rr^ =(~)'='™'' (w^tI + U^T^) 



(B149) 



,1,2 



=0, for TOj; = odd and tti^, 



odd 



There are six PSG's of form 7/2K... whose ansatz have 
a form 



^x.z+m 



(B150) 



Z2A[rlr2,r|T!]riV!!. 



Their ansatz have a form Eq. (B94). For PSG 



(B151) 



(B152) 



Z2AT"rOT-V° 



the ansatz have a form Eq. (p93). For 



Z2A[TlTl,TlTiy^Tl (B153) 

the 180° rotation and the time reversal symmetries re- 
quire that 









u 


-m — 


3 3/ 


.)- 










- 


^m — 


3 3 






which 


gives 


us 




















Uj 


,i+m 


=uln 


t1 + 


uLr' 










ull^ 


= 0, 


for 


m - 


: even 



For PSG 



Z2AtUIt\1 



(B154) 



(B155) 



the ansatz has a form Eq. ( B100| ). 

In summary, we find that there are 12 classes of per- 
turbations around the SU (2) -linear spin liquid that break 
the SU{2) gauge structure down to a C/(l) gauge struc- 
ture, and there are 58 classes of perturbations that break 
the SU{2) gauge structure down to a Z2 gauge structure. 
The resulting C/(l) and Z2 spin hquids can be constructed 
through Uij . 
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